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The monochromatic x-ray diffraction patterns for liquid nitrogen, nitric oxide and nitrous 
oxide have been obtained and analyzed by the Fourier integral method to obtain the atomic or 
electronic distributions. The diffraction patterns show one main peak and only faint peaks or 
plateaus at larger angles. The distribution curves show diatomic aggregates in liquid nitrogen 
and nitric oxide comparable to those found in the gaseous states. A linear molecule for nitrous 
oxide is consistent with the present data while a triangular molecule is not. 


INTRODUCTION 


HE diffraction patterns produced by x-rays 
when scattered by liquids are of interest 
both in a qualitative way and as regards the in- 
formation they can give concerning the ‘“‘struc- 
ture” of the liquid. The x-ray diffraction patterns 
produced by monatomic liquids may be analyzed 
in a systematic way to give the radial distribu- 
tion of atoms about a given atom, using the 
theory of Zernicke and Prins,! and Debye and 
Menke.? The diffraction pattern produced by a 
polyatomic liquid is not so directly interpreted 
in terms of the distribution of atoms, but does 
give, by a direct analysis, the distribution of 
electrons in the liquid. A successful analysis by 
either of these theories depends first on’a careful 
determination of the diffraction pattern out to 
large angles for that liquid, the relative intensity 
of the scattered x-rays, as well as the position of 
the diffraction peaks, being of importance. 
* Now at University of Arkansas, Fayetteville, Arkansas. 
1 F. Zernicke and J. Prins, Zeits. f. Physik 41, 184 (1927); 
Zeits. f. Physik 56, 617 (1929). 


* P. Debye and H. Menke, Physik. Zeits. 31, 797 (1930); 
Ergeb. d. Tech. Rontgenk. (1931), Vol. 2. 


The present work was undertaken both to seek 
to improve on earlier work on the diffraction 
pattern produced by liquid nitrogen and for the 
first time actually to carry out analyses of this 
pattern and those of liquid nitric oxide and ni- 
trous oxide so as to obtain direct information 
concerning the atomic and electronic distribu- 
tions in these liquids. 

The diffraction pattern produced by liquid 
nitrogen was one of the first to be obtained, pic- 
tures having been taken by Keesom and de 
Smedt? about 1923. Using molybdenum radiation 
filtered through zirconium and an evacuated 
camera of 2.75-cm radius, they found two peaks 
in the pattern for liquid nitrogen, the main 
peak at sin @/A=0.139 and a second peak at 
about 0.208. 


THEORY AND METHOD OF EVALUATING 
THE FOURIER INTEGRAL 
The equation used for liquid nitrogen has been 
developed previously‘ in the literature and may 
3 W. H. Keesom and J. de Smedt, Proc. Amst. Akad. Sci. 


26, 112 (1923). 
4B. E. Warren, J. App. Phys. 8, 645 (1937). 
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be written 
4rR%p(R) = 4 Rips +2R/x [ sti(s) sin rsds, 
0 
where 


s=4r sin 0/X, 4(s) =(I— Nf?) /Nf?, 


I is the experimentally determined scattered in- 
tensity which has been corrected for absorption 
in the sample, the partial polarization of the 
scattered x-rays and the incoherent scattering 
background, Nf? is the scattering which would be 
expected if the scattering by the different atoms 
did not interfere, and po is the average density of 
atoms (atoms per cubic angstrom). The function 
4rR’p(R) gives the distribution of atoms about 
any given atom in the liquid. 

The equation used for the polyatomic liquids 
nitric oxide and nitrous oxide has been given in 
connection with the analysis of certain com- 
pounds and glasses.5 The development of this 
equation is based on the assumption that the 
dependence of the atomic scattering factors upon 
sin 6/X for the different atoms in the liquids is 
such that they may be taken proportional to an 
average structure factor f, per electron with Z 
(or K, an effective number of electrons per atom) 
as the proportionality factor. It should be pointed 
out that this assumption is particularly good for 
compounds of nitrogen and oxygen, which ele- 
ments are adjacent in the periodic table. The 
equation is 


¥v Kn4aRYom(R) =44R%pp Eu Zn 
42R/x f si(s) ein Reds, 
0 


where i(s) =(I—N >omfn?m)/Nf2, I is the fully 
corrected experimentally determined scattering 
intensity, V >> uw fm? is the scattering which would 
be expected if the scattering from the different 
atoms did not interfere, fe= ou fm/d-m Zm, and 
po is the average density of the electrons (elec- 
trons per cubic angstrom). >> a signifies summa- 
tion over the unit of composition in the liquid. 
The graph of }\u Kn4aR?pn(R) vs. R thus ob- 
tained is the result of a superposition of the 
electronic distribution curves about each atom 


5B. E. Warren, H. Krutter, and O. Morningstar, J. Am. 
Ceramic Soc. 19, 202 (1936). 
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of the composition, each curve multiplied by the 
effective number of electrons in the origin atom 
for that curve. Thus this curve is more difficult 
to interpret unambiguously than one obtained 
for liquids composed of one kind of atom. 

The experimentally determined and fully cor- 
rected intensity curve was placed on an absolute 
basis by choosing a value of N so that J= Nf? 
(or, for the generalized method, J=N > w fm?) 
out at some large angle where it is assumed that 
all interference effects have ceased. The Fourier 
integral appearing in these two equations was 
evaluated in two ways, one using graphical in- 
tegration and the other using the method of 
trigonometric interpolation which is described 
elsewhere. ® 

The appearance of discrete peaks in the dis- 
tribution curves is of significance. The area under 
the first peak of the 47 R’p(R) curve gives directly 
the number of nearest neighbor atoms at this 
position for liquids of one kind of atom. Two 
methods suggest themselves for interpreting dis- 
crete peaks in the 47R? }> Kmpn(R) curve in 
terms of atomic distributions. The first is to 
calculate from the measured area of a peak the 
number of atoms of type j surrounding one atom 
of type 7 using the formula »=area/2n;K;K, 
where 7; is the number of atoms of type 7 in the 
smallest possible unit of composition. The other is 
the calculation from an assumed plausible struc- 
ture of the area of the expected peak using the re- 
lation that the area is equal to yw Km>-.K; 
where >>, K, serves to add up all the electrons due 
to atoms at a distance R’ (inside a shell with R’ as 
average radius) from atom m and > yw K~» serves 
to move the origin of the shell over all atoms of 
the composition. This last method was of special 
value in interpreting the results for nitrous oxide. 


EXPERIMENTAL 


The apparatus consisted primarily of the x-ray 
source, the system for recording the scattered 
x-ray intensity, and the cooling system for main- 
taining the liquid at a known témperature and 
pressure. Auxiliary apparatus also had to be pro- 
vided for obtaining and purifying the substances 
to be studied. The source of radiation was a 
molybdenum target tube with the Ka radiation 


* G. C. Danielson and C. Lanczos, J. Frank. Inst. 233, 
365, 435 (1942). 
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TABLE I, 
System Pressure Temp. Cooling agent 
Ne Closed, photographic, 3.3 atmos. 89°K Liquid air 
counter 
Nez Closed, photographic 3.3 “ 64°K Liquid air and 
pump 
NO Closed, photographic 2 ” 125°K Liquid methane 
(natural gas) 
N20 Closed, photographic 3 ” 196°K Dry ice and ether 
N:O Closed, photographic 3 . 173°K Dry ice and ether 
and pump 














being reflected from rocksalt. For the photo- 
graphic work a General Electric CA-2 tube was 


used and with the Geiger-Mueller counter sys-: 


tem,’ used to obtain valuable supplementary 
data on liquid nitrogen, a demountable grid- 
controlled molybdenum target tube*® was used. 

Most of the work was done photographically 
with a specially designed® camera. A glass Dewar 
and sample holder 8 mm in diameter were sealed 
to the camera and the camera was evacuated. 
The x-rays monochromatized by the rocksalt 
crystal were collimated by two slits and the width 
of the incident beam at the sample was about 1.5 
mm. The photographic film to record the inten- 
sity of the diffracted x-rays was placed around 
the outside of the cylindrical evacuated camera, 
the x-rays passing through a thin (0.004’’) alu- 
minum window to reach the film. The radius of 
the camera was 9.53 cm. The “‘camera’’ used with 
Geiger counter detection for liquid nitrogen con- 
sisted of a metal Dewar with a sample chamber 
having flat beryllium windows.'® 

The cooling system used for the liquids re- 
ported on here used a closed sample system, as 
contrasted with the open system used on liquid 
oxygen.!! The closed system has the great ad- 
vantage that it doesn’t require as large quantities 
of the substance to be studied as is true for the 
open system. Furthermore, the pressure on the 
gas being studied could be increased to produce 
liquefaction in some cases where the cooling 
agent used would not liquefy the gas at atmos- 
phere pressure. The cooling agents used were 
liquid air, liquefied natural gas (high percentage 


7A. Eisenstein and N. S. Gingrich, Rev. Sci. Inst. 12, 
582 (1941). 
8 A. Eisenstein, Rev. Sci. Inst. 13, 208 (1942). 
(19en) Eisenstein and N. S. Gingrich, Phys. Rev. 58, 307 
>. 
(ean Eisenstein and N. S. Gingrich, Phys. Rev. 62, 261 
). 
1 P. C. Sharrah and N. S. Gingrich, J. Chem. Phys. 10, 
504 (1942). 


51, 3194 (1929). 





methane) and dry ice. The liquid air and dry ice 
were used as the cooling agents for nitrogen and 
nitrous oxide, respectively. Nitric oxide was the 
most inconvenient of the gases to liquefy since 
it is a solid at liquid-air temperature and since 
its critical temperature is below the temperature 
of dry ice. Liquid methane has a boiling point 
which made it of value as a cooling agent and the 
natural gas available was a rather high percent 
methane so that it was used. A brass liquid-air 
container was supported in the reservoir of the 
glass Dewar to serve as a condenser for the na- 
tural gas which would fall to the bottom of the 
Dewar to boil and maintain a fairly constant 
temperature at the sample. (There was a gradual 
rise of this temperature due to fractional dis- 
tillation as the liquefied natural gas boiled away.) 

Table I summarizes the temperatures and cool- 
ing agents used for the various liquids studied. 
The temperatures listed for nitrogen and nitrous 
oxide are probably correct to within two degrees. 
The temperature for nitric oxide is only an aver- 
age, the rise during an exposure being from 8° 
to 10°. A pump was used to reduce the pressure 
and hence lower the temperature in certain cases, 
as shown in Table I. 

The nitrogen used was of two qualities, one the 
usual commercial product, and the other some- 
what more pure (99.8%).* Commercial nitrogen 
was passed through sulfuric acid and phosphorous 
pentoxide in being admitted to the sample holder 
for the photographic method. The purer nitrogen 
was passed through a dry ice trap in admitting 
it to the sample holder for use with the Geiger 
counter system. 

Nitric oxide was obtained from the reaction of 
sulfuric acid dropped into a dilute solution of 
sodium (or potassium) nitrite and potassium 
iodide.’ * The gas resulting from this reaction 
was condensed and frozen out in a liquid-air trap. 
But even with great care a large percent of im- 
purities appeared with N.O; and N,.O, being 
easily detected. The first purification step sought 
to eliminate most of the N,O; and N.O, by per- 
mitting a large portion of the liquid in the trap 


*Ohio Chemical Company, 1177 Marquette Street, 
N. E., Cleveland, Ohio. 

122 W. A. Noyes, J. Am. Chem. Soc. 47, 2170 (1925). 

13H. L. Johnston and W. F. Giauque, J. Am. Chem. Soc. 
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Fic. 1. (a) The x-ray diffraction patterns of liquid nitrogen at 89°K and at 64°K. (b) The atomic distribution 
curves for liquid nitrogen at 89°K and at 64°K. 


to boil off and be condensed in a second liquid-air 
trap, the remaining liquid in the first trap being 
rich in these impurities which then were dis- 
charged from the system through a second tube 
with stopcock. The remaining amount of these 
two impurities was then further greatly reduced 
by boiling the liquid in one trap and freezing it 
out in the other liquid-air trap, thus driving it 
repeatedly back and forth through a potassium 
hydroxide wash and a dry ice trap placed be- 
tween them. The purified nitric oxide was finally 
frozen out and held in reserve in one of the liquid- 
air traps from which it could be introduced di- 
rectly into the closed sample holder where it 
would be maintained in the liquid state by the 
liquefied natural gas used as a cooling agent there. 

Nitrous oxide (N2O) of 99%* purity was passed 
through an ice and salt trap in which the tem- 
perature was around —10°C. This would tend to 
eliminate any traces of moisture, N2O3 and N2O, 
which might be present. 

Five pictures were taken of liquid nitrogen at 


* See asterisk footnote, p. 437. 


89°K and four at 64°K. Three pictures were 
taken of liquid nitric oxide in the neighborhood 
of 125°K. Pictures were taken of liquid nitrous 
oxide at both 196°K and 175°K, but the differ- 
ence between the resulting patterns was so slight 
that the three pictures taken at the lower tem- 
perature were the only ones used in the analysis. 
A large number of patterns were determined for 
each substance because of the broad peaks or 
plateaus beyond the main peak, which were 
typical of all the liquids studied. These relatively 
small peaks are of importance in the analysis and 
thus had to be determined with considerable 
accuracy. 


DISCUSSION OF RESULTS 


Figures 1(a), 2(a), and 3(a) show the experi- 
mentally determined diffraction patterns fully 
corrected for absorption and polarization.“ For 
the patterns obtained photographically, the ab- 
sorption correction used was that for a cylin- 
drical sample completely bathed in a parallel 


4S. Katzoff, J. Chem. Phys. 2, 841 (1934). 
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Fic. 2. (a) The x-ray diffraction pattern of liquid nitric oxide at 125°K. (b) The electronic distribution 
curve for liquid nitric oxide at 125°K. 
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Fic. 3. (a) The x-ray diffraction pattern of liquid nitrous oxide at 173°K. (b) The electronic distribution 
curve for liquid nitrous oxidesat 173°K. 
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TABLE II. Diffraction pattern peak positions in sin 0/2. 








First Second Third 
0.139 0.208 

Nitrogen 89°K .144 .26 0.42 
Nitrogen 64°K 148 .26 42 
Nitric oxide 125°K 15 26 44 
Nitrous oxide 173°K 155 .30 50 


Substance 





Nitrogen* 








* See footnote 3, p. 435. 


beam of x-rays, as discussed by Blake.'® How- 
ever, it has been shown!* that for sodium, with 
cylindrical samples having diameters of from one 
to five times the width of the beam, there is little 
difference in the ratio of the intensity of scatter- 
ing at the main peak to that at angles as large as 
were used here. The absorption correction ap- 
plied to the patterns obtained for nitrogen using 
the Geiger-Mueller counter was for a flat sample. 
These curves are averages from several experi- 
mental patterns. Table II gives the position, in 
terms of sin@/A, of the intensity peaks (or 
plateaus). Values for nitrogen from the previous 
values of Keesom and de Smedt® are given for 
comparison. The weak plateau appearing at 
about sin 6/A=0.42 in Fig. 1(a) was not reported 
by them. The position of the main peak and the 
existence of their second peak are open to ques- 
tion. It is believed that with the larger camera 
here used the present peak positions are likely 
to be more nearly correct. 

Figures 1(b), 2(b), and 3(b) show the atomic 
and electronic distribution curves obtained from 
analyses of the intensity patterns. These curves 
are averages, more than one complete Fourier 
analysis having been carried through for each 
substance. In all the analyses carried out for a 
given substance, the final distribution curves 
were fairly reproducible up to nearly 6A, but 
beyond that distance the various analyses differed 
noticeably. The shape of the discrete peak ob- 
tained for all the liquids analyzed depended 
somewhat on the fitting, and hence its area was 
altered to some extent. A feature of all the dis- 
tribution curves which persisted to a greater or 
lesser extent, depending on the fitting, but which 
never completely disappeared, is the rise im- 
mediately following the fall to zero of the distri- 
bution curve after the first peak. This rise leveled 


1 F, C. Blake, Rev. Mod. Phys. 5, 180 (1933). 
16 G. A. Mitchell, unpublished work. 
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off in the case of nitrogen and nitrous oxide but 
did not resolve itself into a discrete peak. 

The peak positions for the distribution curves 
are given in Table III. 

Resolving the more distant aggregations into 
two peaks (third and fourth) for nitric oxide and 
nitrous oxide is somewhat artificial, but conveni- 
ent for comparison with the graphs for nitrogen 
(and oxygen"). An interatomic distance of 1.095A 
has been reported elsewhere for the nitrogen 
molecule and 1.15A for the nitric oxide molecule."” 
Distances within the triatomic nitrous oxide 
molecule N — N —O have been reported as 1.15, 
1.23A, and 2.38A, respectively.'® Recent values 
for these same distances in nitrous oxide, based 
on electron diffraction experiments, are 1.12A, 
1.19A, and 2.31A, respectively.'® 

The area of the first peak of the distribution 
curve for liquid nitrogen at both 89° and 64°K 
is 1.03 atoms. The area of the first peak in the 
distribution curve for nitric oxide is 108.9 elec- 
trons” and for nitrous oxide 225.5 electrons.’ 

The area under the first peak for nitrogen is 
only slightly different from one, so that, within 
the experimental error, each atom in the liquid 
has one permanent nearest neighbor at the dis- 
tance of the first peak, which is comparable to 
the interatomic distance found by other methods 
for the diatomic nitrogen molecule. 

In order to interpret the meaning of the dis- 
crete peak for nitric oxide and nitrous oxide, the 
formula n=area/2n:K;K; may be used to calcu- 
late the number of nearest neighbors from the 
known area of the peak in electrons.? Applying 
this equation to nitric oxide and seeking to calcu- 
late the number of nitrogen atoms surrounding 
each oxygen atom, one has 2;=no=1. Further- 


TABLE III. Atomic concentration positions in 
angstroms. 








Substance Second Third Fourth 
Nitrogen 89°K “a 2.6 4.0 
Nitrogen 64°K e 2.6 4.0 
Nitric oxide 125°K a (2.6) 3.4 
Nitrous oxide 173°K ; SF | 3.4 











17 G, Herzberg, Molecular Spectra and Molecular Structure 
(Prentice-Hall, 1939) p. 491. 

18 W. G. Penney and G. B. B. M. Sutherland, Proc. Roy. 
. Soc. London, A156, 672 (1936). 

19 Verner Shomaker and Robert Spurr, J. Am. Chem. 
Soc. 64, 1184 (1942). 
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more, one assumes plausible values for the 
effective number of electrons per atom which 
may be taken as K;= Ko=8.3 and K;=Ky=6.7. 
Substitution gives n=(108.9)/2(1)(8.3) (6.7) = 
0.978, which is close to one which is the number 
of nitrogen atoms per oxygen atom, showing that, 
within the experimental error, nitric oxide shows 
a diatomic aggregation in the liquid state with 
the interatomic distance shown by the position 
of the first peak in the distribution curve. 

Now applying the same equation to nitrous 
oxide to seek to determine the number of nitrogen 
atoms around one oxygen atom at the distance 
of the first peak, one has n;=o=1, and using the 
approximate values for K;=Ko=8 and K;=7, 
one gets m=225.5/2(1)(8)(7) =2.05=2. This re- 
sult is two, within the experimental error, and 
it suggests a triatomic unit in the liquid state, 
but the fact that the scattering from nitrogen 
and oxygen is comparable demands some fur- 
ther calculations. With no information about the 
structure assumed in other work for N.O, one 
might reasonaly assume two structures; one an 
approximately equilateral triangular array of 
three atoms with the side of the triangle equal to 
the peak position on the distribution curve, and 
the other a linear array of three atoms in the 
molecule, with the short interatomic distances 
approximately the same and equal to the first 
peak position in the distribution. Assuming these 
two, the expected area of the first peak will be 
calculated according to the equation area= 
Xiu Km Sa Ka and compared with the observed. 

For the triangular array, using Ko=8 electrons 
and Ky=7 electrons, the area of the peak at 
R’=1.2A is Ko(Kynw +Kwn«@)) +Kynua (Ko+ Knew) 
+Ky 2) (Ko +Kwy)) = 8(14) +7 (15) +7 (15) = 322 
elect.2 225.5 elect.?, and since the observed area 


was 225.5 elect.?, the equilateral triangular array 
is rather definitely eliminated. 

The area for the linear array Na)—Ne—O, 
using the same equation is Ko(Kw@))+KN@ (Ko 
+KwNy)+Kw (Kn) =8(7) +7(8+7) + (7) (7) 
= 209 elect.2= 225.5 elect.?. Recalling that the 
peak areas are probably not known to better than 
10 percent this makes the observed 225.5 elect.* 
approach 209 elect.? almost within the experi- 
mental error. The observed distribution curve is 
thus found to be consistent with the existence of 
the linear molecule N — N —O in the liquid state 
comparable to the linear molecule found in the 
gaseous state.1* The present data alone, how- 
ever, could not decide between N—N—O and 
N-—O-—N, the latter requiring an area of the first 
peak of 224 elect.?, this being itself in favor of the 
latter. Now assuming the linear molecule to be 
correct, it is interesting to calculate the area of 
the peak which one would expect due to the 
large N—O distance (2.4A). Its area would be 
KwKo+KoKn=2(7)(8) =112 elect.?. This peak 
is evidently masked by the interpenetration 
within this distance of the atoms of adjacent 
molecules in thermal motion. 
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The Rankine magnetic balance has been used to determine the magnetic susceptibility of 
oxygen and nitric oxide at atmospheric pressure in a field of about 15 oersteds. The mass 
susceptibility of water at 20°C is assumed to be —0.7200X10-*. The mean of ten separate 
measurements on oxygen gives (106.3+0.2) X 10-* for the mass susceptibility, with an average 
deviation from the mean of 0.9 10~*. The value of xT is 0.997 +0.002; the average deviation 
from the mean being 0.008. The mean of six separate measurements on nitric oxide gives 
47.2 10~ for the mass susceptibility, and 0.0590 x 10~* for the volume susceptibility, with an 
average deviation from the mean of 0.6X 10~* and 0.0007 x 10~, respectively. 





INTRODUCTION 


N 1934 Rankine! described a magnetic balance 
suitable for the comparison of feebly magnetic 
materials in extremely weak fields. He later con- 
structed such a balance but it did not function 
with the high sensitivity anticipated, largely be- 
cause of ferromagnetic impurities in some of the 
structural materials.2 An improved form of the 
balance which functioned with the requisite 
sensitivity was first developed by Iskenderian,' 
who used it to determine the susceptibility of 
oxygen in a field of about 40 oersteds.* His 
measurements gave the value 0.979 for xm7, 
which he interpreted as being in good agreement 
with the theoretical value’ of 1.000,° considering 
that the gas used was tank oxygen which prob- 
ably contained some diamagnetic impurities. 
The present work was undertaken to obtain 
more precise measurements on oxygen and to 
obtain measurements on nitric oxide at low fields. 
Oxygen is of particular interest since the theo- 
retical value of xu7J can hardly be questioned, 
and is the value obtained by experiment, within 
1 percent or so, at high field strengths. If, there- 
fore, measurements relative to water were made 
at low fields and a close check with theory ob- 


* Now at Eastern New Mexico College, Portales, New 
Mexico. 

1A. O. Rankine, Proc. Phys. Soc. 46, 1 (1934). 

2 A. O. Rankine, Proc. Phys. Soc. 46, 391 (1934). 

* Haig P. Iskenderian, Phys. Rev. 51, 1092 (1937). 

4 Haig P. Iskenderian, Phys. Rev. 52, 1244 (1937). 

5J. H. Van Vleck, The Theory of Electric and Magnetic 
Susceptibilities (Oxford University Press, London, 1932), 


p. ' 

6 Calculated from Van Vleck’s formula, using recent 
values of the physical constants. (See Raymond T. Birge, 
Rev. Mod. Phys. 13, 233 (1941).) 


tained, it would remove any doubts as to a dis- 
crepancy between the susceptibility of water in 
high and low fields and would furnish additional 
evidence of the fitness of the Rankine balance for 
magnetic measurements on gases. 

Nitric oxide was chosen because the theoretical 
value of the susceptibility is known and no results 
on this gas at low fields have been reported. 
However, the mass susceptibility has been meas- 
ured at high fields by several investigators with 
results as follows: 46.68X10-* (Weiss and Pic- 
card’); 48.7 10-® (Bauer, Weiss, and Piccard’); 
48.8010-* (Sone®). Van Vleck!® has made a 
careful calculation of the volume susceptibility 
using data from the molecular spectrum of the 
molecule. The value® obtained at 20°C and 760- 
mm pressure is 0.0597 X 10-®. The corresponding 
volume susceptibilities from the data of Bauer, 
Weiss, and Piccard and of Sone are 0.0609 x 10-° 
and 0.0610X10~-°, respectively, or a little more 
than 2 percent above the theoretical value. 


THE RANKINE BALANCE 


The Rankine balance depends for its action 
upon the fact that when a test specimen is placed 
near a suspended magnet the magnet experiences 
a force arising from a distortion of its field by the 
specimen. This force is proportional to the sus- 
ceptibility of the specimen, and is balanced by 
the application of a small non-uniform field over 


7P. Weiss and A. Piccard, Comptes rendus 157, 916 
(1913). 

8E. Bauer, P. Weiss, and A. Piccard, Comptes rendus 
167, 484 (1918). 

® T. Sone, Sci. Rep. 11, 139 (1922). 

10]. H. Van Vleck, Phys. Rev. 31, 587 (1928). 
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the magnet, produced by current in a nearby 
wire. From the magnitudes of the currents neces- 
sary to balance the forces arising from various 
test samples their relative susceptibilities can be 
deduced. The procedure in determining the mag- 
netic susceptibility of a liquid or gas with this 
type of balance has been given elsewhere.® 
Figure 1 is a schematic diagram of the balance 
used in the present work. The magnet C is of an 
alloy called Nipermag and was obtained from 
Cinaudagraph Corporation. The weight is 0.5 g 
and the pole strength about 10 e:m.u. The beam 
Fisa very thin Duralumin strip about 6 cm long, 
to which is attached a counterweight G and a 
small plane electrometer mirror E for measuring 
the deflection. The torsion head H is of OFHC 
copper and is supported by a tripod made of 
Pyrex tubing. The base plate A is also of OFHC 
copper and is secured toa heavy hard wood board 
by means of four short OFHC copper legs. The 
lengths of the upper and lower quartz fibers are 
about 24 cm and 17 cm, respectively, and their 
torsion constants of the order of 10-* and 10-° 
dyne-cm per radian. The test cell B is of Pyrex 
and is approximately a cylinder about 4 cm in 
diameter and 1.5 cm long. The suspended system 
was enclosed in the heavy Pyrex cylinder JN, 
provided with a plane glass window and a plate 
glass top which could be sealed on with soft wax. 
Preliminary work indicated that radiation fall- 
ing on the instrument apparently affected the 
suspension. To guard against this, a polished 
copper cylinder M was fitted inside the Pyrex 
cylinder and the upper fiber enclosed in the 
heavy-walled copper tube K. To avoid trouble- 
some electrostatic effects, a thin copper box was 
placed over the test cell to shield it from the 
magnet. The cell-connecting tubes and filling 
arrangement, as well as the balance proper, were 
protected from temperature fluctuations by Celo- 
tex housings. Continuous circulation of air within 
the housings was necessary to avoid small tem- 
perature gradients which affected greatly the 
operation of the balance. The temperature of the 
room was kept constant to about 0.1°C. The mo- 
tion of the suspended system was controlled 
from an adjoining room by varying the current 
in a straight wire beneath the base of the instru- 
ment. This control current was measured by 


means of a standard resistance and a L& N type 
K potentiometer. 

As was expected, it was almost impossible to 
obtain any consistent readings during the day 
when the laboratory was in general use. Conse- 
quently, measurements were made late at night 
when conditions were best. During a run, which 
required about two hours, it was necessary to 
keep under strict control the motions of all 
laboratory apparatus containing iron. It was ob- 
served that the parking of cars in certain loca- 
tions outside the building produced a noticeable 
deflection of the beam. Even alternating currents, 
under certain conditions, were observed to affect 
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Fic. 1. Magnetic balance. A—OFHC copper base plate; 
B—test cell; C—magnet; D—quartz fibers; E—electrom- 
eter mirror; F—Duralumin beam; G—counterweight; 
K—copper tube; H—torsion head; L—Pyrex tripod; M— 
copper shield; N—Pyrex cylinder. 





444 


TaBLe I. The magnetic susceptibility of oxygen at 20°C 
and 760-mm pressure. 








Trial (49) 9, X10" (At) py X10" « X106 x X10 Xl 





105.7 
107.3 
107.9 
104.8 
105.9 
106.6 
106.5 
105.8 
104.3 
103.8 
104.0 
107.4 


0.1407 
0.1429 
0.1437 
0.1395 
0.1407 
0.1420 
0.1418 
0.1409 
0.1389 
0.1382 
0.1384 
0.1430 
0.1397 104.9 
0.1415 106.3 
0.0012 0.9 
0.0003 0.2 


9959 
9830 
9830 
9821 
9866 
9867 
9867 
9897 
9892 
9943 
9826 
9779 
9829 


1893 
1907 
1909 
1857 
1878 
1890 
1887 
1876 
1857 
1860 
1840 
1893 
1863 
Mean values 

Av. deviation from mean 
Probable error in mean 


—_ 
SO WD SIO Cr Go to 








* Deleted from the final calculations. 


the balance. However, only a very few trials were 
vitiated by such disturbances when measure- 
ments were made late at night. 


OXYGEN 


Oxygen was prepared by electrolysis of water 


acidulated with sulfuric acid. After leaving 
the generating apparatus the gas was bubbled 
through distilled water, then passed through a 
tube containing copper oxide heated to 400°- 
500°C. After passing through tubes containing 
calcium chloride and phosphorus pentoxide, the 
gas was collected in a 500-ml bulb, from which it 
could be led directly into the test cell by turning 
a stopcock. 

Table I shows the results of thirteen consecu- 
tive measurements on oxygen. In this table 
(At)o, is the difference between the control cur- 
rents, in amperes, with oxygen in the test cell and 
with the cell evacuated. (Az) w is the correspond- 
ing quantity for water. The last three columns in 
the table give, respectively, the volume suscepti- 
bility, the mass susceptibility, and the product 
of the molar susceptibility by the absolute tem- 
perature. Those trials marked with an asterisk 
have been deleted in calculating the mean values 
and probable error. There was reason to believe 
that the gas had become contaminated when 
these measurements were made. Before making 
the next run (No. 12) the entire system was 
pumped out and refilled. 


A. BURRIS AND C. D. HAUSE 


NITRIC OXIDE 


Nitric oxide was prepared by the reduction of 
nitrous acid by mercury." The mixture used was 
a 20 percent sulfuric acid solution containing 2 
percent by weight of sodium nitrite, with enough 
mercury to provide the necessary surface for a 
reasonable rate of evolution of the gas. After 
leaving the generating flask, the gas was bubbled 
through mercury to remove nitrogen dioxide and 
some of the higher oxides probably present, and 
to prevent too violent a reaction by maintaining 
some pressure on the reacting solution. The gas 
was dried by passing over phosphorus pentoxide. 

Analysis of the gas was carried out by absorb- 
ing the nitric oxide and measuring the volume of 
the residual gas. Two different absorbing solu- 
tions were used; the first a saturated solution of 
sodium sulphite,” made slightly alkaline with 
potassium hydroxide, and the second a saturated 
solution of potassium dichromate," mixed with 
one-fifth its volume of concentrated sulfuric acid. 
The average purity of the gas as indicated by this 
method of analysis was 99 percent nitric oxide by 
volume, with very little variation. 

The correction to the observed value of x was 
made assuming the additivity law and also that 
the residual gas was entirely nitrogen. This latter 
assumption is supported by the work of Moser," 
Weiss and Piccard,’ and Sone.® The value of x at 
20°C was calculated assuming Curie’s law. 

Table II shows the results of seven consecutive 
measurements on nitric oxide. The symbols, here 
applied to nitric oxide, have the same significance 


TABLE II. The magnetic susceptibility of nitric oxide at 
20°C and 760-mm pressure. 








Trial (41) No 107 (Ai) Ww X10 x X108 





0.0581 
0.0540 
0.0606 
0.0592 
0.0580 
0.0594 
0.0587 
0.0599 
0.0007 


46.5 
43.2 
48.5 
47.4 
46.4 
47.6 


756 9796 
699 9740 
789 9740 
757 9534 
742 9538 
763 9580 
748 9551 47.0 
Mean values 47.2 
Average deviation from mean 0.6 








* Deleted from the final calculations. 


11 F, Emich, Monatsh. f. Chemie 13, 73 (1892). 

2 E. Divers, J. Chem. Soc. 75, 82 (1899). 

13G. von Knorre, Die Chem. Ind. p. 531 (1902). 
14L. Moser, Zeits. f. anal. Chemie 50, 401 (1911). 
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as those of Table I. Trial No. 2 is definitely out 
of line with the other values and has been deleted 
in calculating the mean and average deviation 
from the mean. 


DISCUSSION OF RESULTS 


The performance of the Rankine balance in 
this experiment has been found thoroughly satis- 
factory, with the precision in a single determina- 
tion limited largely by accidental disturbances 
which could not be controlled. The data show 
that the effect of these disturbances is an entirely 
random one, with the maximum deviation of any 
one measurement from the mean of less than 2 
percent and an average deviation of less than 1 
percent in the case of oxygen. The corresponding 
figures for nitric oxide are only slightly higher. It 
is felt that considerable significance may, there- 
fore, be attached to the mean values. These 
values are, in the authors’ opinion, about as 
precise as can be obtained with the balance 
in its present form under ordinary laboratory 
conditions. 

Measurements on electrolytic oxygen give a 
mean value of 0.997 for x7, which is in excellent 
agreement with the theoretical value of 1.000. 
The agreement with the results obtained by other 
methods using high field strengths is also good. 
The average of the results of seven other investi- 
gations,® using different methods, is 0.990, with 
a maximum deviation from the mean of about 3 
percent. 

In view of the results on oxygen, it is difficult 
to reconcile, on the basis of experimental error, 
the 3 percent difference between the value ob- 
tained here and that obtained by Sone for nitric 
oxide. The highest single value obtained in the 


%E. C. Stoner, Magnetism and Matter (Methuen and 
Company, London, 1934), p. 342. 


present work is a little more than 0.5 percent 
below that given by him. 

An important factor which enters into meas- 
urements on nitric oxide is the question of the 
purity of the gas. Most measurements have been 
made on the more or less impure gas and the 
observed value of the susceptibility corrected for 
the impurities. In the present analysis the two 
different absorbing solutions gave results in 
agreement to 0.2 percent. Since it is usually 
agreed that absorption methods are reliable, espe- 
cially when used on a gas of high purity, it is 
probable that the 99 percent purity stated is very 
close to the correct value. It therefore seems im- 
probable that a more accurate analysis would 
change the susceptibility by more than a fraction 
of 1 percent. The gas used by Sone was 95 percent 
nitric oxide by volume. He states: ‘“‘. . . thus the 
correction for the susceptibility due to nitrogen 
was about 3 percent.”’ According to the authors’ 
calculations, the correction would be nearer 5 
than 3 percent. However, this would result in 
an even greater discrepancy between the two 
values. 

It will also be noted that the volume suscepti- 
bility obtained here is about 1.2 percent below 
the theoretical value. Van Vleck® points out in 
his paper that x in his equation x«=hAv/3kT, 
might be about 5 percent larger, due to the fact 
that the value of B is different in the two- 
component normal states of the molecule. He 
states that by a crude estimate this would affect 
the susceptibility by about 1 percent. It is ap- 
parent from his final equation that an increase in 
x lowers the susceptibility. Hence, it is possible 
that his value 0.0597 X 10-® may be as much as 1 
percent too high. This would bring the value 
0.0590 X 10-* into excellent agreement with the 
theoretical’ value, while it would appear that 
0.0610 X 10-°, obtained by Sone, is too high. 
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A simplified method for the determination of the depolarization factor of Raman lines is 
reported which requires only one exposure and is, therefore, independent of fluctuation of arc 
current and any other factors influencing the intensity of Raman lines. Two Polaroid films are 
placed at the front of the slit of the spectrograph, one above the other. The electric vectors of 
the radiation are at right angles to one another. Between the Polaroid films and the slit are 
introduced two half-wave mica plates (one above the other) which orient the electric vectors of 
the light so that they are both vibrating in the same direction upon striking the prism faces. 
Thus differential reflection at these faces is eliminated. The method was checked with carbon 


tetrachloride. 


INTRODUCTION 


INCE the discovery of the Raman effect, 
workers in this field have been interested in 

the determination of the factor of depolarization 
for the Raman lines of the liquids they have 
studied. This information is necessary to a better 
understanding of the type of vibration or rotation 
involved. Cabannes and his co-workers! were 
among the first to devise an accurate method for 
this type of measurement. The diffuse radiation 
was broken up into its two components by means 
of a double-image prism and then depolarized by 
means of a suitably orientated crystal. The two 
components were then focused one above the 
other on the slit of the spectrograph. Trumpy’ 
and Rao*® also made use of the double-image 
prism. Reitz‘ carried out a very intensive study 
of the method. Both Cabannes and Reitz made 
use of the calcite combination of Hanle.’ Gupta® 
developed a unique arrangement. A double-image 
prism was used in such a way that the entering 
light was vibrating at an angle of 45° from the 
principal plane. It would leave the crystal vi- 


*We gratefully acknowledge a fellowship in chemistry 
from the Allied Chemical and Dye Corporation of which 
ok Harold T. Baker was the recipient for the year 1942- 
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1J. Cabannes and co-workers, Ann. de Phys. [X] 19, 
229 (1933). 

2B. Trumpy, Zeits. f. Physik 88, 226 (1934). 

3A. V. Rao, Proc. Ind. Acad. Sci. 2A, 49 (1935). 

4A. W. Reitz, Zeits. f. physik. Chemie B33, 368 (1936). 

5 W. Hanle, Ann. d. Physik 11, 885 (1931). 

¢ J. Gupta, Ind. J. Phys. 10, 313 (1936). 


brating in the same way, but the separation of 
the two components is in the vertical plane or in 
the direction of the length of the slit. Thus the 
differential reflection of the prism was greatly 
reduced. Langseth et al.7 used a Wollaston prism 
behind the slit of the spectrograph in the col- 
limator tube. This arrangement resulted in a dis- 
placement of one of the components so that they 
were recorded in different positions on the plate. 
Bhagavantam® interposed a Nicol between the 
Raman tube and the slit of the spectrograph. 
Two spectrograms were taken on the same plate 
with the vibration direction of the Nicol parallel 
to the track of the light in the liquid and the 
other perpendicular, the conditions of illumina- 
tion being identical in both cases. Corrections 
were introduced for the differential reflections at 
the prism faces. Cleveland and Murray® used a 
fixed film of Polaroid between the slit and the 
Raman tube. The position of the arcs was 
changed by 90° in two exposures. The tube con- 
taining the filter solution acted as a lens causing 
the light from the arc to converge onto the Ra- 
man tube. Edsall and Wilson" also employed the 
tube of filter solution as a lens. They used plane- 
polarized light, obtained with a Polaroid film, 
as the primary source of light in the first exposure 
and unpolarized light in the second exposure. The 


7A. Langseth, J. U. Sorensen, and J. R. Nielsen, J. 
Chem. Phys. 2, 407 (1934). 
(1983) Bhagavantam, Ind. J. Phys. 5, 60 (1930); 7, 79 

9F. F. Cleveland and M. J. Murray, J. Chem. Phys. 7, 
396 (1939). 

10 J. T. Edsall and E. B. Wilson, Jr., J. Chem. Phys. 6, 
124 (1938). 
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intensity of the unpolarized light was reduced so 
as to obtain equal densities for lines having a 
depolarization factor of 6/7. Stitt and Yost! in- 
troduced a half-wave mica plate before the Nicol 
prism. In the first exposure the diffuse rays 
passed through the mica plate vibrating in the 
same direction as upon entering and in the 
second exposure, the mica plate was so orientated 
that the rays were rotated through 90° upon 
leaving the crystal. Some of the methods recorded 
above involve one exposure only, which, however, 
is of great length and they are not simple in exe- 
cution, while others having features of simplicity 
require two exposures. It was desirable to devise 
a method which would feature both of these 
desiderata. The use of a Polaroid film’® and a 
mica half-wave plate! are most convenient and 
both are used in the method described below. 
However, the particular arrangement developed 
here allows one to obtain results with one 
exposure only. 

Two separate experiments are obviously un- 
desirable because all the experimental conditions 
for the two runs must be kept constant. The light 
fluctuations, change in the liquid (decomposi- 
tion), temperature, humidity conditions, etc., 
cannot be easily controlled. With a single ex- 
posure, however, these conditions may vary 
without affecting the results. 


METHOD 


An attempt has been made to develop a 
method that would not demand any elaborate 
equipment and would, at the same time, satisfy 
most of the requirements necessary to give a 
fairly accurate quantitative measurement of the 
depolarization factor. For the determination of 
this factor, the following simple method was em- 
ployed. Two pieces of Polaroid were cut in such 
a way that one of them permitted light vibrating 
vertically and the other one allowed radiation 
vibrating horizontally to pass. The first piece was 
placed above the other one in front of the slit. 
Between each Polaroid film and the slit was 
placed a half-wave-length mica plate. One of 
these mica plates allows the light to pass through 
unchanged in its direction of vibration. The other 
one causes the light vector to be rotated through 


"F. Stitt and D. M. Yost, J. Chem. Phys. 5, 90 (1937). 
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90°. In the latter case, the mica half-wave plate 
was so orientated that it separated the polarized 
light passing through into its two components, 
one of which was retarded a half-wave-length 
with respect to the other. When these two waves 
leave the mica, they recombine in a wave that is 
plane polarized, but whose direction of vibration 
is at an angle of 90° from its original direction. 
The arrangement is shown diagrammatically in 
Fig. 1. 

The two components are measured simul- 
taneously and are found on the plate, the one 
above the other. The Polaroid films and mica 
plates are mounted in a special metal slide that 
can be pushed into place in the front of the slit 
of the spectrograph. This equipment is used only 
for depolarization measurements and is removed 
for regular Raman runs. 

With many compounds, the time of exposure 
for ordinary Raman work may take as long as 
twenty or more hours. When the ideal arrange- 
ment of Cabannes' and Reitz‘ is used for de- 
termination of the depolarization factor, the in- 
tensity of the light becomes many times less and 
so almost prohibitive lengths of exposure are 
necessary. To get around this difficulty, Placzek 
and Van Wijk,” Ornstein and Stoutenbeck," 
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Fic. 1. P=Polaroid 
films, M=mica_half- 
wave plates, S/=slit, 
Sp=spectrograph. Ar- 
rows indicate electric 
vector. Dotted line in- 
dicates height of slit. 
Thickness of mica plate 
is 0.048 mm. 
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es Placzek and W. Van Wijk, Zeits. f. Physik 67, 582 
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754 (1933). 
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Fic. 2. Raman apparatus. A = arcs, F = filters, 
R=Raman tube, S/=slit (carrying Polaroid films and mica 
plates), L=lens, G=grid, Sp=spectrograph. Side view 
=longitudinal cut through the axis of the Raman tube. 


and Nijveld and Gerding,’* have employed 
Wood’s arrangement in order to obtain greater 
light intensity. Use was made of liquids whose 
depolarization factor (p) was known by the ideal 
method. The values obtained by this method 
were plotted against those obtained with Wood’s 
arrangement. Since highly polarized lines have 
a depolarization factor of about 0.4 when this 
method is used, and lines that are somewhat 
less polarized give almost identical results, it is 
nearly impossible to differentiate very accurately 
between lines with p=0.05 and p=0.30. Lines 
that are near 6/7 in value can be determined 
accurately because they are almost completely 
depolarized, and so the orthogonality of the sys- 
tem is not such a crucial factor. 

If an arrangement were used somewhere be- 
tween the ideal and the conventional Raman 
set-up, then it should be possible to obtain values 
lower than 0.4 and thus be able to differentiate 
better between highly polarized and moderately 
polarized lines. The experimental arrangement 
here described is an attempt in this direction. 


RAMAN EQUIPMENT 


The general features of the apparatus used are shown in 
Fig. 2. The placement of the Raman tube, filter solutions, 


“W. J. Nijveld and H. Gerding, Rec. trav. chim. 59, 
1198 (1940). 
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and arcs is similar to that described by Cleveland and 
Murray.’ The two arcs, the two filter tubes, and the Raman 
tube lie in a horizontal plane—all being parallel to one 
another. It was found convenient to build a completely 
adjustable metal frame to hold the Raman tube, filter 
solutions, arcs, and grid. The exciting radiation from the 
arcs passes through a grid of blackened tin strips which 
keeps the light nearly perpendicular to the Raman tube. 
Two elliptical, polished metal pieces were placed behind 
the arcs in order to increase the intensity of the light. The 
metal framework and equipment were mounted on a plat- 
form, which moved on a track, making it easy to push the 
equipment up to the spectrograph. The outside diameter of 
the arcs and Raman tube was 15 mm and that of the filter 
tubes was 21 mm. The distance from the center of the 
Raman tube to the center of a filter tube was 33 mm and 
the distance from the center of the filter tube to the center 
of an arc was 35 mm. The front plate of the Raman tube 
and the slit of the spectrograph were about 17 to 18 cm 
apart, and in this region was placed a lens of focal length 
4 to 4.5 cm, which converged the scattered light onto the 
slit. The ratio of image size to object size was approxi- 
mately 1: 1. 

The support of the arcs, filters, Raman tube, and grids 
consisted of a rectangular metal frame. Two pieces of steel 
(1’"* 3’"X41") were each bent twice at right angles to 
form open end rectangles of sides 14, 13, and 14”, respec- 
tively. These two open rectangles were placed one above 
the other with a ;;’’ opening between them. A piece 
(235’’X3""X 13") was welded onto the open ends of these 
rectangles. The other two corners were made rigid by 
welding on two small rectangular strips (2 75" 3/3”). 
This unit gave a framework with slots on three sides per- 
mitting the sliding of the equipment. A sliding center piece 
was made by taking two pieces (1’’ X }’’ 17”) and bending 
each of them at right angles (both bends inward) at a dis- 
tance of 24” from each end. These pieces were placed one 
above the other and spaced ;¢” apart, and welded together 
with two strips (3X %4’’ x3”). This section can slide in 
the slots of the frame perpendicularly to the 14’’ sides and 
parallel to the 13” sides. 

Into another piece (3X 1" X47%’") two 7%” holes and 3” 
hole were drilled. The half-inch hole is located 4” from one 
end and the other smaller holes are needed to attach the 
piece to the middle of the solid 13” side of the frame. A 
metal rod (10” long and 3” in diameter) is fitted into the 
hole just described and soldered or welded in place. This 
unit is the support for the back half of the framework. The 
front end is arranged similarly, except that a suitably sized 
metal tube is attached to a metal plate. When the frame is 








Fic. 3. Raman 
assembly. y=tube 
and filter support; a, 
b=blackened grid 
pieces. 
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fixed onto the platform it can be rotated about the axis of 
the Raman tube. The 14” sides are parallel to the Raman 
tube axis. 

To support the Raman tube, arcs, and filter tubes, ad- 
justable metal supports move in the slots of the 13” front 
piece and the 13’’ movable piece. These supports are 
elongated Y’s as shown in Fig. 3. A bolt with a wing nut 
holds the support to the frame. It can be raised and 
lowered or moved sideways. Thus any desired spacing of 
the Raman tube, arcs, and filter tubes can be secured. The 
grid of metal strips which gives the exciting light a direction 
perpendicular to the Raman tube was also made adjustable, 
so that changes in dimensions could be obtained when re- 
assembling the grid. Two types of pieces as shown in Fig. 3 
are the main parts. Blackened washers serve as spacers. 
The ‘“‘a”’ metal pieces (33’"X 3” or 33’ X §”) form the light 
grid and are made of thin tin sheet. The “‘b”’ pieces (7’’ X 8”) 
are of heavier metal and thus give it more rigidity. They 
are made by taking two metal strips and soldering a small 
metal piece across each end. In order to assemble the grid, 
four long bolts (%’’ diameter) are passed through each of 
two “‘b” pieces. The “a’’ pieces are then placed perpen- 
dicularly onto the ‘‘b’”’ pieces, while the washers serve as 
spacers. Two ‘‘b’’ are necessary for every four ‘‘a’’ pieces. 
This arrangement is repeated until the desired width is 
obtained. Different sizes can be made to suit other purposes. 


SPECIAL ITEMS 


Two useful techniques in connection with Raman work 
which have been found very useful are ultrafiltration and 
the employment of a Tyndall beam in order to test for 
optical clearness. Solutions of salts cannot be purified by 
distillation and so filtration must be resorted to. Since 
ordinary filter paper will not remove much colloidal matter, 
ultra-filtration was employed to great advantage as de- 
scribed by Holmes.“ The Tyndall beam is secured by 
focusing the light from a carbon arc into a narrow beam. 
The Raman tube is then held in this beam and observed at 
right angles. If colloidal matter is present in only small 
amount, the beam can scarcely be seen as it passes through 
the liquid. Distilled water is a good reference standard. 


TABLE I. Depolarization factors (p) of carbon 
tetrachloride. 








Raman Accepted 
line values 
(em~!) Run I Run II Run III Average Ref. (17) 


218 0.80 0.83 0.82 0.82 0.857 

314 81 85 84 83 857 

459 07 .06 .09 .08 05 
762, 790 82 79 .80 .80 719 











®H. N. Holmes, Manual of Colloid Chemistry (John 
Wiley and Sons, New York, 1934) third edition, pp. 28-30. 
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Fic. 4. Pyrex mer- 
cury arc. 


After a number of trials, a very satisfactory design for 
a mercury arc was evolved and used for the polarization 
measurements. The glow discharge tubes containing neon 
and mercury,’ while found very satisfactory for the meas- 
urement of Raman lines, did not furnish enough radiation 
when their illumination was cut down by the grid. The 
exposure times became quite prohibitive in length. For this 
reason, a change was made to a very satisfactory design of 
a Pyrex mercury arc. When such an arc operates at 
moderately cool temperature, very little continuous back- 
ground is found on the photographic plates. Low tempera- 
ture operation is, therefore, an important fact to consider 
in their construction. Arcs having large heat capacity art 
able to dissipate energy more readily than those of smaller 
size. A large mercury reservoir for each of the two electrodes 
will accomplish this purpose. The arcs are built in the shape 
of the letter ‘‘H’’ as shown in Fig. 4. The cross bar (15 mm 
Pyrex tubing) serves as the source of illumination. When 
these arcs are cooled by a breeze from fans and from com- 
pressed air jets, very efficient light sources are obtained. 


RESULTS 


When the depolarization factor was measured 
with the described arrangement, set to give the 
best possible orthogonality and the least possible 
convergence error, the results obtained for carbon 
tetrachloride were in good agreement with the 
accepted values!’ as shown in Table I. 

It appears that the method is sufficiently ac- 
curate for all ordinary purposes of Raman spec- 
troscopy. The relative ease of operation and the 
speed of execution commend it to the investigator 
of Raman spectra. 


16 G, Glockler and C. E. Morrell, J. Chem. Phys. 4, 15 
(1936). 

17F. F. Cleveland and J. M. Murray, Res. Pub. Illinois 
Inst. Tech. 2, 1 (1941). 
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Raman frequencies, relative intensities, and depolarization factors are listed for 1-iodo-1- 
propyne, 1-bromo-1-propyne, 1-chloro-1-propyne, and 3,3-dimethyl-1-butyne. The relative 
intensities and depolarization factors were obtained by use of a Gaertner microdensitometer. 
Assignments of the frequencies to the different vibration types and calculated values of the 
heat capacities for the ideal gaseous state at 1 atmos. pressure are given for 1-iodo-1-propyne 
and 1-bromo-1-propyne. Examination of the results in the 2200-cm™ region leads to the sug- 
gestion that the resonance splitting of the triple bond fundamental for disubstituted acetylenes 
may frequently be due to the combination frequencies (2900—700) and (1375+700); and 
that the doubling of the 2100-cm™ fundamental for the monosubstituted acetylenes, 3-methyl-1- 
butyn-3-ol and 3,3-dimethyl-1-butyne, may be due to the second overtones of the highly po- 


O C 


larized, “‘breathing’’ frequencies near 700 cm™ of the groups C—C— and ct, respectively. 
| | 


INTRODUCTION 


ISUBSTITUTED acetylenes almost in- 
variably have two or more lines near 2200- 
cm, rather than the expected single line corre- 
sponding to the triple bond carbon vibration fre- 
quency. Badger! has suggested that the multi- 
plicity is caused by a Fermi resonance interaction 
of the acetylenic fundamental with an overtone. 
However, not a great deal of success has accom- 
panied attempts to point out which vibrations 
are involved in the interaction. Crawford? in his 
study of dimethylacetylene ascribed the doubling 
to the first overtone of the infra-red active funda- 
mental at 1126 cm. Murray and Cleveland’ by 
use of both Raman and infra-red data for methyl- 
phenylacetylene found it necessary to assume a 
rather large anharmonicity to account for the 
phenomenon. 

The splitting of the 2200-cm— fundamental is 
so general for disubstituted acetylenes that the 
fundamental or fundamentals involved must oc- 
cur consistently in a wide variety of these com- 


* Presented at the American Physical Society meeting, 
ro ae University, January 1943; Phys. Rev. 63, 220 
). 
1R. M. Badger, J. Chem. Phys. 5, 178 (1937). 
2B. L. Crawford, Jr., J. Chem. Phys. 7, 555 (1939). 
*M. J. Murray and Forrest F. Cleveland, J. Chem. Phys. 
8, 133 (1940). 
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pounds at approximately the same frequency. 
The carbon-hydrogen frequencies‘ near 1375 and 
2920 cm and probably also the carbon stretch- 
ing frequency near 700 cm satisfy this condition 
reasonably well. The (1375+700) combination 
frequency falls near the 2200-cm— fundamental 
but the difference tone (2920—700) comes even 
nearer to coincidence. A number of acetylenes 
for which spectra have been reported show a line 
near 700 cm~, but it is missing in others. In the 
authors’ own work, for example, this line is not 
found in those cases for which the spectra as a 
whole are easily obtainable with short exposures. 
The weakness of the line corresponding to the 
700-cm~! fundamental does not, however, neces- 
sarily preclude its interaction with the 2920-cm™' 
fundamental to give the resonance effect with the 
2200-cm™ acetylenic frequency. 

Since rather simple symmetrical-top molecules 
are possible in the disubstituted acetylenes and 
since all of the fundamentals for molecules of this 
symmetry are permitted in the Raman spectra 
by the selection rules, it was hoped that com- 
pounds of this type would aid in solving the 
problem. The molecules selected were methyl- 


‘It was suggested to the authors, see J. Am. Chem. Soc. 
62, 3185 (1940), that hydrogen atoms on adjacent carbons 
might in some way be responsible for the doubling. 
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iodoacetylene (1-iodo-1-propyne) and methyl- 

bromoacetylene (1-bromo-1-propyne). The spec- 

trum of t-butylacetylene (3,3-dimethyl-1-butyne) 

was also obtained since it had been observed® 

that a monosubstituted acetylene, dimethyl- 

ethynyl carbinol (3-methyl-1-butyn-3-ol), with 
O 


the group — adjacent to the triple bond 


c 


gave a double line, 2089 and 2112 cm-. It was 
expected that ¢-butylacetylene would show a 
similar doubling, which indeed it did, although 
the intensity order of the lines differs in the two 
compounds. 


EXPERIMENTAL 


The spectrograms used in obtaining the frequencies and 
the polarization spectrograms were made with the tech- 
nique and apparatus described in previous papers from this 
laboratory.6 The spectrograms were obtained with the 
samples in the liquid state and excitation was by Hg 4358A. 
The intensities were determined by use of a Gaertner micro- 
densitometer. This instrument was calibrated for the Agfa 
Superplenachrome Press film by replacing the Raman tube 
by a white cardboard diffuser and varying the distance of 
an auto-lamp from the cardboard in order to get known 
relative intensities of illumination. One-hour exposures 
were made for each position of the lamp. A calibration 
curve was then made by plotting the densities, at a wave- 
length corresponding to Ayv~1400 cm™, against the relative 
intensities. 

The inverse square law method, described above, was 
subsequently checked by a method involving the use of 
crossed Nicol prisms. Finally, the calibration curve was 
checked, with satisfactory results, by obtaining Raman 
spectrograms of various known mixtures of ethyl bromide 
and carbon tetrachloride. 

The relative intensities are only semi-quantitative since 
no correction was made for variation of the sensitivity of 
the Agfa film with wave-length. This correction was, of 
course, not necessary for the depolarization factors. 

The methylacetylene used for the preparation of the 
halogen derivatives was made by the method of Heisig and 
Davis’ using propylene bromide with potassium hydroxide 
in n-butyl alcohol. The methylacetylene was liquefied in a 
dry-ice trap. Excellent yields were obtained. An outline of 
the method of preparation of the other compounds follows: 


* Forrest F. Cleveland and M. J. Murray, J. Am. Chem. 
Soc. 62, 3185 (1940). 

_* Forrest F. Cleveland, M. J. Murray, J. R. Coley, and 
V. I. Komarewsky, J. Chem. Phys. 10, 18 (1942). 

*G. B. Heisig and H. M. Davis, J. Am. Chem. Soc. 57, 
339 (1935), 


OF ACETYLENES 


1-Iodo-1-propyne, 


CH;C=CI, C;H3;I: Two-fifths mole of magnesium was 
converted into Grignard reagent, using ethyl bromide in 
ethyl ether. The flask containing the Grignard reagent was 
provided with stirrer, a dry-ice cooled reflux condenser, and 
a tube through which a slight excess of methylacetylene was 
introduced. After evolution of ethane ceased, a saturated 
solution of iodine in dry ether® was added, with stirring, 
until a slight color remained. Water was added, the ether 
layer separated, and the methyliodoacetylene purified by 
distillation through a Podbielniak-type column. B. p. 109- 
110°C. The compound is colorless when freshly distilled, 
but turns slightly yellow upon standing. It is much more 
stable to air than is the corresponding bromo derivative. 


1-Bromo-1-propyne, 


CH;C=CBr, C;H;Br: This compound was synthesized 
by a modification of the method which Straus, Kollek, and 
Heyn® used for the preparation of several other bromo- 
acetylenes. Potassium hypobromite was prepared by addi- 
tion of 35 cm? of bromine to 150 g of potassium hydroxide 
in one liter of water at 0°C. Through this solution (in two 
flasks, in series) was bubbled the methylacetylene obtained 
from 150 g of propylene dibromide. Considerable methyl- 
bromoacetylene formed in the flask was swept over by un- 
reacted methylacetylene. Part of the product was caught 
in a tube cooled by ice water; the remainder and somewhat 
more than half of the methylacetylene were received in a 
dry-ice trap. The recovered methylacetylene was again 
passed through the hypobromite solution and the process 
repeated until almost all had reacted. The lower layer from 
the hypobromite and the condensate in the traps were com- 
bined and distilled in an atmosphere of nitrogen. B. p. 
64.5°C at 750 mm. This agrees with the boiling point of a 
sample prepared by Loevenich, Losen, and Dierichs!® by 
dehydrobromination of 1,1,2-tribromopropane. 

The freshly distilled product is not spontaneously com- 
bustible, but the residue, or the pure compound, ignites 
quickly when it comes in contact with air and gives dense 
clouds of black smoke having an irritating odor. All trans- 
fers of the material from flasks to Raman tube, etc., were 
made in an atmosphere of nitrogen. 


1-Chloro-1-propyne, 


CH;C=CCl, C3;H;Cl: The authors were unable to dis- 
cover any reported synthesis of methylchloroacetylene. An 
attempt to prepare the compound by the action of methyl- 
acetylenesodium (freshly prepared sodamide and methyl- 
acetylene followed by removal of ammonia) on p-toluene- 
sulfonyl chloride in dibutyl ether was unsuccessful. Another 
synthesis, using methylacetylene-magnesium bromide in 
dibutyl ether, resulted in a very small yield. The Raman 


8 V. Grignard and H. Perrichon, Ann. chim. 5, 5 (1926). 

°F. Straus, L. Kollek, and W. Heyn, Ber. 63B, 1868 
(1930). 

10 J. Loevenich, J. Losen, and A. Dierichs, Ber. 60B, 950 
(1927). 
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TABLE I. Raman spectra of 1-iodo-1-propyne and 1-bromo-1-propyne.* 








CH;C =CI 


Av Assignment 


CH3C =CBr 


Av Assignment 





163 E 
+343 ’ E 
Ai 
a “ 163-+343 =506, 
Ai+A2+E 
691 Ay 


f E 
| 
eo | 1438 — 343 = 1095, 

Ai+A2+E 
1375 3 bs 


1438 E 
oes [ 1375-4691 = 2066, A, 
ead 2(1030) = 2060, A: +E 

{ 2892—691 =2201, A; 
2197 ? 

2(1375) =2750, A, 

2730 
aaes ' 2(1438) =2876, Ai +E 
2913 E 
2958 





171 E 
+343 E 
464 0.3 A; 


692 Ay 
1026 
1154 


1375 
1443 


2060 
2213 
2230 


E 
1443 — 343 = 1100, 


Ai+A2+E 


1 


E 
1375 +692 = 2067, A; 
2(1035) =2070, Ai+E 
2900 — 692 = 2208, A, 
A; 


2853 MS } {2(1443) = 2886, 41+£ 
2921 , \ Ai 
2968 8 E 








* Av is the Raman frequency in cm~|, J is the relative intensity based 
upon 100 for the strongest Raman line on the spectrogram, p is the 
depolarization factor, vw (very weak) designates lines that were too 
weak to measure with the microdensitometer, P indicates lines for 
which the weaker component of the polarization doublet was too weak 


spectrum of the small amount of product which was ob- 
tained (b. p. 37-41°C) showed that ethyl bromide was the 
main impurity. In working up the product, considerable 
quantities of methylacetylene and of unreacted p-toluene- 
sulfonyl chloride were recovered. The yield of the product, 
however, was not increased by heating on a water bath for 
several hours the vigorously stirred reaction mixture con- 
taining the acetylenic Grignard and the sulfonyl chloride. 


3,3-Dimethyl-1-butyne, 


(CH3)sCC=CH, Cs.Hno: Fifty grams of pinacolone (East- 
man Kodak Company product) were converted to the 
dichloride by action of 110 g of PCs." After the addition of 
water the dichloride was isolated by steam distillation. To 
obtain ¢-butylacetylene, the method of Heisig and Davis’ 
was again employed, with entirely satisfactory results. The 
solid pinacolone dichloride was introduced, a small amount 
at a time, into 200 cm’ of n-butyl alcohol and 90 g of po- 
tassium hydroxide. The ¢t-butylacetylene was caught in a 
dry-ice trap. Distillation through a Podbielniak-type 
column yielded a product whose boiling point was 37.5- 
38.0°C. A small amount of olefinic impurity was present in 
the final sample, as indicated by a weak Raman line at 
1630 cm, 


RESULTS 


The experimental results for 1-iodo-1-propyne 
and 1-bromo-1-propyne, which were obtained in 
the pure state, are listed in Table I. It was not 


11 A, Faworsky, J. Prakt. Chem. 37, 393 (1888). 


to measure with the microdensitometer, D means depolarized, + indi- 
cates lines that were obtained both as Stokes and as anti-Stokes lines, 
and the parentheses inclose data in regard to which there is some 
uncertainty. 


possible to obtain pure 1-chloro-1-propyne, but 
the data for the impure sample, which contained 
ethyl bromide and a trace of some olefinic im- 
purity, are listed in Table II. The lines that are 
certainly due to 1-chloro-1-propyne are indicated 
by asterisks. Also included in Table II are the 
data obtained for the 3,3-dimethyl-1-butyne. 


DISCUSSION 


The symmetry of each of the methyl halogen 
acetylenes concerned in the present paper is C3); 
the t-butylacetylene also has this symmetry if 
one of the hydrogen atoms of each methyl group 
lies in one of the vertical symmetry planes. Using 
group theory methods, one can work out not only 
the selection rules for the fundamentals, com- 
binations, and overtones but also the number of 
fundamentals of each vibration type.” The re- 
sults are summarized in Table III, from which it 
is seen that five totally symmetrical A, and five 
doubly degenerate E vibrations are permitted by 
the selection rules for the methyl halogen acety!- 
enes. The A, vibrations should be polarized, the 
E vibrations depolarized, in the Raman spectra. 


For an elementary account of these methods see 
Arnold G. Meister, Forrest F. Cleveland, and M. J- 
Murray, Am. J. Phys. 11, 239 (1943). 
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With the aid of the polarization data it is not 
difficult to pick out the A; and E fundamentals 
for the 1-iodo-1-propyne and 1-bromo-1-propyne; 
the assignments are indicated in Table I. One 
selection which may be somewhat questionable 
is the choice of the 1375-cm™ line as an A, 
fundamental, since it is not very highly polarized. 
However, there are five E type fundamentals and 
only four A, fundamentals, if this line is excluded, 
and it thus seems fairly certain that it must be 
of the A, type. Crawford was confronted with 
the same difficulty in his analysis of methy]l- 
acetylene (another symmetrical-top acetylene) ; 
there the corresponding line was listed as 
‘1380(2)D?”, but was nevertheless chosen as an 
A, fundamental. 

In the 2200-cm~ region, three lines were ob- 
served for each of the compounds. These could 
arise from the resonance interactions of the 
triple-bond fundamental with the combinations 
and overtones indicated in Table I. In addition 
to the approximate coincidence of the combina- 
tion or overtone frequency with the fundamental 
frequency, it is necessary that a symmetry condi- 
tion be met in order for resonance doubling to 
appear. This symmetry condition is that the 
overtone or combination must have the same 
symmetry as the fundamental, or at least have 
one component which is the same as that of the 
fundamental.!? Thus the combinations and over- 
tones listed in Table I for the 2200-cm- region 
can produce resonance interaction with the A, 
triple-bond fundamental because they are either 
of type A; or contain A; as a component. The 


TABLE IT. Raman spectra of 1-chloro-1-propyne and 
3,3-dimethyl-1-butyne. 








CH3C =CCI+ 


(CH3)sCC =CH C2H;Br 


I p 


(CHs)sCC =CH 





8 


- 
OR AR W AUN 


*1380 











* These lines belong to the 1-chloro-1-propyne molecule. 
** This line is due to a trace of some olefinic impurity. 


SB. L. Crawford, Jr., J. Chem. Phys. 8, 526 (1940). 
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TABLE III. Selection rules for any molecule of sym- 
metry C3, and the number of vibrations of each type for 
CH;C =CX and (CH;3);CC =CH. 








Activity* No. of fundamentals 


Raman Infra-red | CHsC=CX (CHs)sCC=CH 





Fundamentals i a : 5 10 
ia i 0 4 
a F 5 14 


Combinations 
tAiXAi 
AiXAe2 
AiXE 


Overtones 
ttAi" 
A (n even) 
A (n odd) d 
E2 Ait+E 


Bs Art+Act+E 

















* a =active, ia inactive. 

+ Ai XAj1 represents a combination frequency (vitvj) of an A1 
frequency vi with another A: frequency »j. 

tt Ai" represents the (7 —1)th overtone of an A: frequency. 


combinations (1438+691)=2129 and (2958— 
691) =2267 for 1-iodo-1-propyne and (1443+ 
692) = 2135 and (2968 — 692) = 2276 for 1-bromo- 
1-propyne are of type E, and hence by the 
symmetry restriction cannot cause resonance 
doubling of the triple-bond fundamental. Several 
other cases of resonance doubling appear in 
Table I and it will be noted that in each case the 
symmetry condition is fulfilled. 

Since the CH deformation vibration near 1375 
cm, the CH stretching frequency near 2920 
cm~, and the carbon stretching frequency near 
700 cm~! are present in a great number of di- 
substituted acetylenes, the combinations (1375 
+700) and (2920—700) represent two other 
possibilities for the resonance doubling in addi- 
tion to the overtone of the carbon stretching 
frequency near 1100 cm™ as proposed by Badger 
and by Crawford.” In addition, as suggested by 
Badger,! the second overtone of the 700 cm fre- 
quency, which would be of type Ai, may some- 
times be involved. Thus there exist four overtone 
or combination frequencies which are satisfactory 
for explaining the multiplicity of lines usually 
observed in the 2200 cm“ region for the disubsti- 
tuted acetylenes. Whether two or more lines are 
observed depends upon the frequency of the 
triple bond fundamental in the compound under 
consideration. 

These possibilities would not apply to mono- 
substituted acetylenes since for them the 700- 





rs. Ox 


TABLE IV. Fundamental frequencies.* 








CH;C =CI CH3C =CBr CH3C =CI CH3C =CBr 
Av D Av Av D Av D 


171 2 1375 1375 
343 1438 1443 
464 2192 2225 
692 2892 2900 
1035 - 2958 2968 





163 2 
343 2 
410 1 
691 1 
1030 2 








* Av =Raman frequency in cm~!. D =degeneracy. 


and 1100-cm™ lines are replaced by lines near 

925 and 3300 cm™, respectively. Hence the few 

cases in which doubling has been observed for a 

monosubstituted acetylene must involve some 

other vibration. Dimethylethynyl carbinol which 
O 


| 
has the group C—C— adjacent to the triple bond 


C 


was previously found® to have a double line at 

2089 and 2112 cm™ and the ¢-butylacetylene 

concerned in the present study also had a double 

line at 2105 and 2138 cm™; the latter has the 
C 


| 
group C—C— adjacent to the triple bond. The 
i 


symmetrical, ‘‘breathing’’ frequency of these 
groups falls at 709 for dimethylethynyl carbinol 
and at 690 cm for ¢-butylacetylene. In both 
cases this frequency is highly polarized and thus, 
at least for t-butylacetylene, corresponds to an 
A, type fundamental whose second overtone 
could produce resonance doubling of the triple- 
bond frequency, which for most monosubstituted 
acetylenes appears at 2118 cm“. 

For t-butylacetylene, the selection rules (see 
Table III) permit 24 frequencies to appear in the 
Raman spectrum. Twenty-three were observed, 
but not all of these can correspond to funda- 
mentals. The polarized lines at 690, 885, 2105, 
and 2923 as well as the lines at 1390 and 3307 
cm! can rather safely be taken as A; fundamen- 
tals, and the depolarized lines at 183, 931, 1205, 
1453, and 2968 as E fundamentals. The lines at 
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2138, 2865, 2887, and 2905 cm- probably owe 
their appearance to resonance interactions. A 
complete interpretation of this spectrum is not 
possible upon the basis of the present data. 


HEAT CAPACITY CALCULATIONS 


Since all of the fundamental frequencies were 
obtained for 1-iodo-1-propyne and 1-bromo-1- 
propyne, it is possible to calculate the heat ca- 
pacities of these compounds for the ideal gaseous 
state at one atmosphere pressure. The tables 
given by Wilson" were used in the calculations. 
The results are only approximate since the fre- 
quencies were obtained for the liquid rather than 
the gaseous state. Because of the several cases 
of resonance interaction, it was necessary to 
make slight adjustments of the frequencies in 
these cases; the adjustments were made by con- 
sideration of the relative intensities of the com- 
ponents of the resonance doublet or multiplet. 
The fundamentals used for the calculations are 


TABLE V. Heat capacities at constant pressure. 








CH;3C =CBr 


16.77 
18.08 
20.56 
24.47 
27.42 
29.65 


CH:C =CI 


16.95* 
18.23 
20.67 
24.56 
27.49 
29.71 


Temperature 


250.0°K 
298.16 
400. 
600. 
800. 
1000. 











* Calories/(mole-degree). 


listed in Table IV; the numbers appearing under 
the letter D are the degeneracies. The values of 
the heat capacities calculated from these fre- 
quencies are given in Table V. 
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14 E. B. Wilson, Jr., Chem. Rev. 27, 17 (1940). 
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Previous discussions of the kinetic theory of rubber elas- 
ticity have dealt with individual long chain molecules, 
but the theory of the structure of bulk rubber has been 
almost entirely undeveloped. The present paper goes be- 
yond earlier ones both in the more detailed treatment of 
the individual chains and in the development of a clear-cut 
model of the bulk material. From the consideration of 
familiar properties of rubber, it is concluded that in the 
lightly vulcanized state it consists of a coherent network of 
flexible molecular chains (this involving a considerable frac- 
tion of the total material) together with other molecules 
not actively involved in the network, but acting like a 
fluid mass through which the network extends and in 
which it moves with Brownian motion. The idea of an 
effective internal pressure is advanced and discussed. A 
simplified model for bulk rubber is proposed, consisting of 
a network of idealized flexible chains extending through 
the material and a fluid filling it, the bounding surfaces 
being in equilibrium under all the forces acting on them— 
internal pressure, pull of the molecular network, and any 
external forces. For the quantitative treatment of this 
model the principal problem is that of computing the forces 
exerted by the molecular network, because of its thermal 
agitation, on the bounding surfaces. Two models of the 
flexible molecular chains have been used—one with a linear 
stress-extension relation (Gaussian chain), the second with 
independent links of fixed length, having, like real molecu- 
lar chains, a definite maximum extension. Methods for the 
statistical treatment of chains of independent links at all 
extensions are developed and applied to the computation 
of stress-strain relations for the second of the above models. 
It is shown that an irregular network of Gaussian chains is 
equivalent to a simple set of independent chains; the corre- 
sponding reduction in the case of non-Gaussian chains is 
only approximate. Making this reduction in all cases, the 
model is applied to the quantitative computation of stress- 
strain curves for unilateral stretch of rubber and for stretch 
in two directions. A prediction is made of the change in 
elastic properties of rubber in the swollen state. The thermo- 
elastic properties of rubber are worked out, and the change 
in slope of the nearly linear stress-temperature curves from 


1. INTRODUCTION 


UBBER is the best known representative of 
a class of materials set apart from other 
*On leave of absence at the Radiation Laboratory, 


Massachusetts Institute of Technology, Cambridge, Massa- 
chusetts. 


negative to positive values with increasing strain is ex- 
plained. The location of the thermoelastic inversion point, 
at which this slope is zero, is shown to depend only upon 
the cubic expansion coefficient of unstretched rubber. The 
linear thermal expansion coefficient of stretched rubber in 
the direction of the stress is shown to be positive (and of 
the order of the same coefficient for unstretched rubber) 
below the thermoelastic inversion point. Above the thermo- 
elastic inversion point this coefficient becomes negative 
and is—for appreciable extensions—of the order of the 
cubic expansion coefficient for a gas, independent of the 
composition of rubber. The linear coefficient of thermal 
expansion perpendicular to the direction of stress is always 
positive. Stress-strain curves for bulk rubber at constant 
temperature have been compared with the theory for exten- 
sions up to 400 percent. The agreement is particularly good 
as concerns that part of the total stress due to changing 
entropy. In some rubber compounds minor deviations, 
variable from sample to sample, are found for small exten- 
sions. These are attributed to van der Waals forces not 
taken into account by the theory. The characteristic knee 
in rubber stress-strain curves is shown to be due to change 
in internal pressure in the material. This varies inversely 
with the extended length, and in unstretched rubber is of 
the order of 5 kg/cm?. The upward curvature of the stress- 
strain curve for larger extensions, which appears even be- 
fore the onset of crystallization, is explained as due to the 
approach of the molecular network to its maximum exten- 
sion, which is of the order of 10 times its extension in the 
unstretched bulk rubber. Crystallization may in general en- 
hance the S-shape of the stress-strain curve for natural 
cured rubber, Neoprene and butyl rubber. Synthetic rubber 
of the Buna type, while having an S-shaped stress-strain 
curve, does not exhibit any crystallization at all. It is esti- 
mated that, in the materials considered, roughly a fourth 
of the rubber chains are actively involved in the network. 
Internal pressure, maximum extensions, and the fraction 
of active material all change with the state of vulcanization. 
The general theory developed in this paper provides a basis 
for the treatment of many other physical properties of 
stretched rubber. 


solids by their characteristic elastic and thermo- 
elastic behavior. The most striking property of 
rubber is its long-range reversible elasticity. This 
property is so characteristic that any material 
possessing it may be called a ‘‘rubber’’ regardless 
of its cheniical constitution. With soft rubber, 
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elongations of some 1000 percent (greater than 
those possible with ordinary solids by a factor of 
10°) are obtainable with relatively small devia- 
tions from complete reversibility; the elastic con- 
stants are, correspondingly, of a quite different 
order of magnitude. The stress-strain curves 
have, also, an unusual S-form, the differential 
Young’s modulus at first decreasing with in- 
creasing strain until it is of the order of one-third 
its value for small strain, and thereafter increas- 
ing more and more rapidly (cf. Fig. 1). 

The thermoelastic properties of rubber-like 
substances are equally striking, and perhaps 
more suggestive as to the mechanism involved. 
Rubber extended under constant load contracts 
when heated; the stress of rubber kept at con- 
stant length increases with rising temperature, 
under not too extreme conditions very nearly 
proportionally to the absolute temperature. The 
striking generation of heat when rubber is adia- 
batically stretched is thermodynamically related 
to these other elements in the behavior of the 
material and needs no special explanation. 

In a general way the existence of long range 
reversible elasticity in rubber and the tempera- 
ture dependence of stress for constant strain have 
been understood for some years;? an attempt? 
has even been made to get an estimate of the 
order of magnitude of the elastic constants. These 
previous discussions have, however, been of a 
very general character, or have involved the use 
of models which were vaguely defined, physically 
unsatisfactory, or incompletely worked out.‘ It 
is Our purpose in the present paper to present a 
definite and reasonable model for rubber-like 
materials in bulk, and to employ it in the quan- 
titative explanation of various properties of 
rubber, especially the hitherto little understood 
form of the stress-strain curves at constant 
temperature. 

Many qualitative conclusions about the mech- 
anism of rubber-like elasticity can be drawn from 


(1938) Guth and H. Mark, Monats. f. Chemie 65, 93 

2E. Guth, Kautschuk 13, 201 (1937). 

*W. Kuhn, Kolloid Zeits. 76, 258 (1936). 

*For a more detailed discussion, with reference to the 
literature, see E. Guth and H. M. James, Ind. Eng. Chem. 
33, 1624 (1941); and H. M. James and E. Guth, ibid. 34, 
1365 (1942). In these two papers a preliminary discussion 
is given of the new theory developed in detail in the present 
paper. 
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Fic. 1. Typical stress-strain curve for rubber. 


the most elementary features in the behavior of 
these materials. 

The behavior of rubber as a solid, the possi- 
bility of elastic deformations, shows that bulk 
rubber possesses a definite structure which de- 
termines its form in the absence of external 
forces. The elements of structure must, however, 
be of a quite different sort from those in normal 
solids, in which bonds of electrostatic, chemical, 
or van der Waals type connect atoms or mole- 
cules into a definite and rigid three-dimensional 
structure. The long-range extensibility suggests 
that the structure of rubber must be of a much 
more flexible, more open, type with the relation- 
ships between atoms less completely defined. The 
small forces needed to produce considerable dis- 
tortions suggest that the elastic restoring forces 
cannot be of any of the types listed above. One 
must suppose that ordinary distortions of bulk 
rubber do not appreciably change the distances 
between atoms linked by chemical or other strong 
bonds. That is, in distorting rubber one does not 
act against the strong bonds which suppress the 
plasticity and maintain the coherence of the ma- 
terial. Only for the very large extensions, at 
which rubber begins to crystallize or approaches 
the breaking point, do these bonds come into 
play. For instance, ordinary chemical bonds ap- 
pear to contribute to the tensile strength of 
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rubber, which (taken with regard to the cross 
section at the breaking point) is about that of 
copper. 

Rubber resembles a liguid in possessing high 
volume rigidity together with small rigidity of 
shape. The compressibility of rubber is that of an 
ordinary liquid. In consideration of the stretching 
of rubber this is quite negligible, and one may 
treat rubber as incompressible for extensions up 
to 300 or 400 percent; Poisson’s modulus is thus 
+ for small deformations. Furthermore, the co- 
efficient of thermal expansion for rubber is of the 
same order of magnitude as for liquids. These 
points suggest that the interactions between ad- 
jacent molecules in bulk rubber are, for the most 
part, similar to those which maintain the fixed 
volume of ordinary liquids. In addition to these 
forces, which do not maintain a fixed coupling of 
molecule to molecule or atom to atom, there 
must be strong bonds between the molecules, 
constituting the elements of structure previously 
mentioned. These bonds may, however, affect 
only relatively small portions of the very large 
rubber molecules, and thus have relatively small 
effect on the above properties of the material. 

Rubber resembles a gas very strikingly in its 
thermoelastic behavior. The proportionality of 
stress to absolute temperature under constant 
strain suggests that the tension in stretched rub- 
ber, like the pressure in a gas, is associated with 
a change in entropy of the material when it is 
deformed, rather than with a change in internal 
energy. Quantitatively, the behavior of stretched 
rubber, or of a gas under pressure, is determined 
by its free energy, which in the case of rubber 
under simple tension may be written 


F(L, T)=U(L, T)—TS(L, T). 


(1.1) 


Here Z represents the length of the stretched 
material and the other symbols have their cus- 
tomary significance. The tension is 


Z(L, T) =(0F/dL)7=(dU/dL) r—T(dS/dL) 7 
=(0U/dL)r+T(0Z/dL)r. (1.2) 


In the case of a gas one has only to replace L by 
the volume V, and Z by the negative pressure 
(—p). The simple and uniform thermoelastic be- 
havior of ideal gases, independent of their chem- 
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ical composition, is due to two characteristics of 
these gases: (1) the internal energy is in all cases 
nearly independent of the volume; (2) the en- 
tropy consists of two parts, one associated with 
the thermal capacity of the gas, dependent on 
the chemical composition, but making no con- 
tribution to (0S/d V)7, and the second associated 
with the number of positional configurations 
available to the system, independent of the com- 
position, and making a contribution to (0S/dV)r 
which is independent of T. The uniform and 
simple thermoelastic behavior of ‘‘ideal’’ rubber- 
like materials, independent of their composition, 
may be understood in a similar way, consistently 
with our previous discussion of their properties. 
The internal energy of the liquid-like mass of 
molecules will be practically independent of its 
form, such strong chemical bonds as may exist 
being but little stretched by the relatively small 
forces which produce moderate deformations of 
the material; (9U/dL) 7 will be small. As we shall 
indicate later in more detail, the entropy will 
consist of two parts, one associated with the heat 
capacity and not contributing to (0S/dL)7, the 
other associated with the number of configura- 
tions available to the molecules within the ma- 
terial. The number of these configurations will 
depend on the external form of the material but 
not on the temperature, on the general structure 
of the masses of molecules but not on their pre- 
cise chemical composition. (0S/0L) r will, accord- 
ingly, be independent of temperature, and we 
will have, approximately, 


Z(L, T)=Z,(L)T. (1.3) 


The setting up of models for bulk rubber con- 
sistent with all these ideas, and the calculation 
of their entropies and the associated forces, is the 
main problem considered in this paper. 


2. MODELS OF THE STRUCTURE OF RUBBER 


In earlier discussions of rubber elasticity as an 
entropy effect, attention has been concentrated 
not on bulk rubber, but on a single rubber 
molecule. 

The rubber molecule is a very long unbranched 
hydrocarbon chain composed of isoprene groups 
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linked end to end: 
| 


JAMES AND E. GUTH 


CH,+-CH:—C=CH—CH;--CH,—C=CH—CH:~-CH: =. 


| 
| CH; 


The first essential characteristic of this molecule 
is its great length. The second is its flexibility. 
Of the carbon-carbon bonds in the chain, three- 
fourths are single, one-fourth double. Now the 
molecule can have its two parts rotated with 
respect to each other at any one of the single 
bonds without change in the energy, to the ap- 
proximation in which interactions are neglected 
between atoms in the chain farther apart than 
third neighbors. Indeed it possesses the possi- 
bility of free rotation, in the usual quantum- 
mechanical sense,’ to the better approximation 
which does take account of interaction of the 
atoms directly bound to adjacent carbon atoms 
in the chain. It is sufficient for the following 
general discussion that there exist configurations 
of the molecule differing by rotations of one part 
with respect to the other about any given single 
bond and having the same minimum energy; the 
configurations may be separated by energy bar- 
riers, provided only that they are not so high as 
to prevent transitions from one configuration to 
another in a time small compared to that in 
which a distortion is applied to the bulk rubber. 
By twisting about all the single bonds, the mole- 
cule can take on® an enormous number of con- 
figurations, all of the same energy, to the ap- 
proximation just mentioned. Given the direction 
of one carbon-carbon bond in the chain, that of 
the next bond is restricted to the neighborhood 
of two or three definite directions. That of the 
second bond along the chain will be less definitely 

5 Cf. for instance, R. H. Fowler and E. A. Guggenheim, 
Statistical Thermodynamics (Cambridge, 1939), pp: 107- 
109. 

® Measurements of the heat capacity of ethane indicate 
that the energy barrier hindering the free rotation of the 
two CH; groups.relative to each other is of the order of 
3 kcal./mole. However, for isoprene the presence of the 
double C=C bonds may very well lower the barrier by a 
factor of 3 or more. In fact for dimethylacetylene 
(CH;—C =C—CH,)), the presence of the triple C=C bond 
lowers the barrier so much that the two CH; groups 
rotate practically freely relative to each other. These 
points seem to have been overlooked by S. Bresler and 
J. Frenkel, Acta Physicochem. U.S.S.R. 11, 485 (1939) in 
their attempt to calculate the influence of the potential 


barrier on the contracting force in single long chain hydro- 
carbon chains. 


CH; 





fixed, there being more different orientations of 
maximum probability, while bonds only a little 
farther down the chain will take on their orienta- 
tions practically independently of the prescribed 
orientation of the first. So long as it is not neces- 
sary to consider in detail the configurations of 
short segments of the molecule, one may treat it 
like a perfectly flexible chain, which can take on 
any configuration without change of energy.’ 
One may treat such a long chain molecule, 
consisting of very many units, by the methods 
of statistical mechanics. To a chain in thermal 
equilibrium with a temperature bath and subject 
only to the constraint that the ends have a 
separation L, one may assign the entropy 


S=k log C(L)+A(T), (2.1) 


where & is Boltzmann’s constant and C(L) is the 
relative number of configurations of the chain 
consistent with the separation L of its ends. The 
term h(T) arises from the thermal capacity of 
the chain, and, because of the flexibility of the 
chain, is independent of L. For a flexible chain 
the internal energy is completely independent of 
the configuration of the chain, and thus of L; 
hence, we have 


u=u(T). (2.2) 


The free energy becomes 


F=u(T)—kT log C(L) —Th(T) 


7For the mathematical development of the theory one 
may introduce the flexibility of the chain in a variety ol 
idealized ways, as by treating the chain as involving only 
single bonds about which perfectly free rotation is pos- 
sible, or as consisting of rigid links each of which can take 
on its orientations independently of its neighbors. These 
are the models of the rubber molecules employed in the 
present work. The general character of the results will be 
independent of the way in which the flexibility is intro- 
duced, but there will be differences in details, and to some 
degree in estimated magnitudes, from model to model. 
Concerning such details, the theory will be reliable only in 
roportion as the model is well chosen for the particular 
ong chain molecules to be treated. In the present paper 
we have chosen the second and simpler model indicated 
above, as it has all the essential characteristics of the real 
chains. The other, and quantitatively more reliable, model 
will be treated in a second, more mathematical, paper. 
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and the tension in the chain 


Z(L, T) = —k-(1/C(L) ]-(dC(L) /dL]-T. (2.4) 


The functions u(T) and h(T) drop out of the 
computation of Z, which is proportional to the 
absolute temperature. This will clearly be the 
case with any system of flexible chains which one 
might use as a model for bulk rubber. 

This problem could also be discussed from the 
point of view of kinetic theory. A chain molecule 
immersed in a temperature bath will undergo 
constant changes in configuration through its 
thermal agitation—an internal Brownian motion. 
For a perfectly flexible and unconstrained mole- 
cule, all configurations will have the same energy 
and the same probability of occurrence. However, 


since knotted and coiled configurations are far’ 


more numerous than relatively straight and ex- 
tended ones, the probability is overwhelming that 
at any given time the molecule will be found in 
a highly twisted state, with its ends separated by 
a distance small compared to the total length of 
the chain. If the ends of the molecule are con- 
strained to lie at fixed points, the thermal agita- 
tion will result in a constant jerking of the mole- 
cule against the constraints. Unless the ends of 
the molecule are coincident, jerks in the different 
directions will not be equally probable; they will 
result in a net force tending to bring the ends of 
the molecule closer together—an effective ten- 
sion. If the temperature is raised, the thermal 
agitation will be more violent and the effective 
tension will increase. If the distance between the 
ends of the chain is slowly increased at constant 
temperature, work must be done against the 
forces in the jerking chain. The energy of thermal 
agitation of the chain will, however, be un- 
changed, and the energy given to the chain will 
appear as heat. In all these respects the single 
molecular chain would show the behavior ob- 
served with bulk rubber, and the same would be 
true of any system of chain molecules so con- 
nected together that the intramolecular Brown- 
ian motion is not suppressed. This approach to 
the understanding of rubber elasticity as due to 
the tendency of long chain molecules to curl up 
under the influence of thermal agitation is quali- 
tatively useful. For quantitative purposes, how- 
ever, methods involving the enumeration of the 


configurations of the chain under the imposed 
constraints are usually the only practicable ones. 

In order to investigate the dependence of Z 
upon L, one must determine the form of the 
function C(L), which will vary as to details from 
model to model. Guth and Mark! treat the case 
of chains of N independent links, each of length 
l, in both one and three dimensions. With L de- 
noting the component of the extension in a fixed 
direction, they find in either case that for small 
L a Gaussian distribution applies: 


C(L) =A exp (—L?/2(L?)w). (2.5) 


For the one-dimensional case (L*)y = N/; for the 
three-dimensional case (L?)y=(1/3)NP. Kuhn® 
considers a three-dimensional chain in which each 
independent link can contribute +//v3 to the 
extension in any or all of the coordinate direc- 
tions, and obtains the same result. Both make the 
extension to the case of a chain of links each 
making a fixed valence angle 6 with its neighbors, 
but otherwise independent of them, by assuming 
the same form for C(L) and noting that in this 
case (L?)y=(1/3)NF tan? 6/2. The limited ac- 
curacy of such results will be discussed in Part 6 
of this paper. 
With this form for C(L), one obtains 


Z=[kT/(L)w]L. (2.6) 


In order to extend this discussion to bulk 
rubber, it is necessary to make some assumption 
as to how the material is built up from single 
chains. The implicit assumption of Guth and 
Mark is that a unit cube of rubber can be treated 
as a bundle of M similar, parallel, and inde- 
pendent chains. These chains would exert equal 
restoring forces; one has at once as the restoring ° 
force in stretched rubber, per unit cross section 
of the unstretched material, 


Z=MURT/L*)w]-L. (2.7) 


This model is much oversimplified. It is ani- 
sotropic, being stretchable only in the direction 
of the bundle, and has length and volume zero 
when acted upon by no external forces. The pro- 
portionality of tension to length given by this 
model is, correspondingly, far simpler than the 
observed relation for rubber-like materials. 
Kuhn gives a rather vague description of his 
model and his formalism is actually inconsistent 
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with his textual picture of the situation. He does 
picture a transition from a hypothetical melted 
rubber, in which the molecules take on all the 
different configurations possible for isolated mole- 
cules, to the actual solid, in which each molecule 
is immobilized in the configuration it possessed 
at the instant of solidification. The picture of a 
complete immobilization of the molecular chain 
is certainly incorrect; such a suppression of the 
internal Brownian motion would destroy the 
rubber-like properties of the material. However, 
Kuhn’s mathematical treatment actually corre- 
sponds to the assumption that, after the freezing, 
each molecule is fixed at the end points and at 
points 4, 3, and 3 of the way along its length, with 
the Brownian motion otherwise unhampered. He 
assumes that, when the bulk material is distorted, 
the configurations of the molecules (in his mathe- 
matical treatment the fixed points) have the mo- 
tion which would be assigned to them by a con- 
tinuum theory of elasticity. He provides no 
adequate picture of the mechanism by which this 
might be brought about; the basis of his discus- 
sion is essentially the assumed variation of the 
distribution function with distortion of the ma- 
terial. This consideration of what happens to the 
material perpendicular to the direction of stretch 
makes the model, in contrast to that of Guth and 
Mark, a three-dimensional one; a complete treat- 
ment of his formalism would, in fact, lead to our 
preliminary form of the stress-strain curve as 
discussed in Section 3. Kuhn, however, deter- 
mines only the slope of the stress-strain curve for 
very small stresses and not its complete form. 
Wall® gives a formal discussion of bulk rubber 
_similar to Kuhn’s, based on an assumed distribu- 
tion function without an associated physical pic- 
ture, and obtains in his second paper the stress- 
strain curve previously described by the authors.‘ 
Kuhn and Wall treat, in an identical way, the 
problem of the ‘‘molecular weight’? of rubber 
molecules in bulk rubber. A different view of the 
situation is presented in Section 7 of the present 
paper. Kuhn’s model and its relation to our 
model is further discussed in that section and 
in the Appendix. 


8F. T. Wall, J. Chem. Phys. 10, 132, 485 (1942). In the 
first paper, the incompressibility of stretched rubber was 
not taken into account. 
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The picture of the structure of rubber on 
which we base our model is as follows.° Raw 
rubber consists of a tangled mass of long chain 
molecules, of which some may be linked together 
by chemical bonds, while others can move more 
or less freely past each other, impeded only by 
the tendency of the Brownian motion to keep 
the chains tangled. The possibility of an un- 
limited relative displacement of the molecules is 
responsible for the plasticity of raw rubber, while 
the tangling of the chains keeps the internal vis- 
cosity high—higher the longer the molecular 
chains. Raw rubber deformed rapidly and im- 
mediately released will return to its original shape 
with little plastic deformation because there has 
been time for only a little motion of the molecules 
with respect to each other.!° 

Vulcanization appears to consist of the rein- 
forcing of cross-links already present in raw rub- 
ber and the introduction of new cross-links."! 
These bonds will be formed at random, some- 
times connecting different parts of a single mole- 
cule, but more often linking parts of two different 
molecules which happen to be adjacent at the 
time. As more.and more bonds are formed, the 
whole mass of molecules will come to be linked 
together in a coherent network, very irregular in 
detail and quite open in structure, there being 
only a few bonds formed for each of the long 
molecules. These bonds introduce fixed relations 
between the molecules which will be undisturbed 
by later distortions of the body, and result in 
the suppression of the plasticity of the material. 
On the other hand, they affect only a small por- 
tion of each molecule; the intramolecular Brown- 
ian motion is largely undisturbed, and over the 


9A convenient summary of the great variety of ideas 
concerning the structure of rubber is provided in the 
monograph, Chemistry and Technology of Rubber, edited 
by C. C. Davis and J. T. Blake (New York, 1937). 

10 The recent work by L. R. G. Treloar, Trans. Faraday 
Soc. 36, 538 (1940), shows that the plasticity of raw rubber 
is largely of a semi-permanent character. On release after 
stretching, the material returns almost completely to its 
original unstretched length, though this recovery may take 
a long time, require additional heat treatment and the 
application of solvents. 

11 The exact chemical nature of vulcanization is irrel- 
evant for the purposes of this paper. The cross-links may 
be van der Waals bonds or chemical bonds. Where the 
latter are concerned, it is very likely that the double C=C 
bonds are involved in the formation of cross-links. At any 
rate, saturated long chain hydrocarbons (which may show 
the properties of raw rubber to a marked extent) cannot 
be vulcanized. 
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greater part of their lengths the molecules will 
be free to move past each other, as in a liquid or 
in raw rubber. Rubber vulcanized to this degree 
will show all the properties of raw rubber except 
the plasticity. If the vulcanization is carried fur- 
ther, the intermolecular bonds will eventually 
become so numerous as to suppress the intra- 
molecular Brownian motion to a great extent and 
to change markedly the relationship of adjacent 
molecules. The extensibility of the material then 
approaches that of ordinary solids. We must ex- 
clude from consideration such highly vulcanized 
rubbers, and even soft rubber under high strains, 
in which there is known to be induced a crys- 
tallization which radically modifies the properties 
of the material. Our discussion will be restricted 
entirely to the case of ideal soft rubber-like ma- 
terials, characterized by negligible plasticity and 
at least approximate proportionality of stresses 
to the absolute temperature. 

The first essential element in the structure of 
rubber of which our model must take account is, 
then, the existence of a complicated network of 
flexible molecular chains extending throughout 
the material. Such a network will, as already 
shown, exert restoring forces proportional to the 
absolute temperature. An idealized network of 
this type would alone, however, be quite inade- 
quate as a model for rubber. For instance, such 
a model would extend through a very small vol- 
ume in the absence of external forces. It is not 
possible to neglect completely, as one does in the 
treatment of idealized flexible chains, the volume- 
filling property of the rubber molecules. 

The volume of the bulk rubber will be regu- 
lated by lateral forces between the chains quite 
unrelated to those effective tensions along the 
chains which tend to maintain the shape of the 
material. The chains, for the most part moving 
freely past each other, will jostle each other, and 
any bounding surfaces, as they would in a liquid. 
The bounding surfaces may be those of other 
solids, or free surfaces of the material at which 
the outward forces of the jostling molecules are 
balanced by the inward pull of the coherent net- 
work. The effect of the sidewise thrusts and pulls 
will be the same as in an ordinary liquid—the 


production of an effective hydrostatic pressure 
in the material, which will depend on its volume. 
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This will clearly be uniform throughout the ma- 
terial so long as this undergoes no irregular dis- 
tortions. It is likely that the mobility of the 
rubber molecules within the bulk rubber is suffi- 
cient to assure the usefulness of the idea of a 
hydrostatic pressure under somewhat more gen- 
eral conditions, but such an extension of our 
model is beyond the scope of the present paper. 

In our model of bulk rubber, we separate and 
idealize the two aspects of the molecular net- 
work. We replace the real molecular network by 
a similar network of idealized flexible chains, 
very irregular in detail, but homogeneous and 
isotropic on the average. The space-filling aspect 
of the network we represent by filling the model 
with a fluid. This fluid can, as a first approxima- 
tion, be taken as incompressible; further refine- 
ments may be introduced by attributing to it an 
appropriate compressibility and thermal expan- 
sion. To prevent the collapse of the idealized 
network into a small volume, we must then 
constrain it to stretch throughout the whole of 
the model by constraining appropriate points of 
it to lie on the bounding surfaces. This completes 
the representation of volume-filling properties of 
the network. Under equilibrium conditions, to 
which we restrict ourselves, every surface of the 
model must of course be in equilibrium under all 
the forces which act on it—the push of the hydro- 
static forces, the pull of the network, and any 
external forces. 

It should be pointed out that not all of the 
molecular chains will be involved actively in the 
network. A chain connected into the network at 
a single point will not contribute to its tendency 
to pull against constraints. The same may be 
said of closed loops in chains in the network, and 
chains threaded through the network but not 
linked to it. Finally, some of the chains may be 
completely dissociated from the main network, 
either removable by solvents, or, if connected 
into large and complicated snarls, not so easily 
separable from it. All these ‘‘inert’’ molecular 
chains will play just the role of the fluid filling 
our model and can be omitted from the idealized 
network. In the final section of this paper we 
shall arrive at an estimate of the fraction of 
the 


the molecular chains actively involved in 


network. 
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It may be objected that, in idealizing the mo- 
lecular chains and neglecting their interactions 
with other chains, we neglect a factor which 
affects the configurations possible to the network, 
and therefore the forces which it will exert. For 
the most part, however, the effect of these other 
chains will be that of a foreign fluid, even though 
they are not of the ‘‘inert’’ sort just described. 
They will not modify the probability distribution 
of the configurations eventually available to the 
chain, since they do not occupy fixed positions, 
and they will thus not modify the forces which it 
exerts. An obvious exception occurs in the case 
of portions of chains near a junction, where rela- 
tive positions are more or less fixed. However, if 
the material is lightly vulcanized a relatively 
small fraction of the total chain length will be 
thus affected, and the approximation should not 
be serious. 

We begin the treatment of this model by a dis- 
cussion of the stress-strain relation which it ex- 
hibits when the distribution function C(L) for 
each segment of molecular chain between two 
junctions is taken to have the Gaussian form 
(Section 3). The principal calculation required 
for this discussion, the reduction of the forces 
exerted by an irregular network of flexible chains 
to those of an equivalent regular network, is to 
be found in the Appendix. Other properties of 
the model with Gaussian chains are discussed in 
Section 4. Agreement with experiment is satis- 
factory for moderate strains, testifying to the 
value of the idea of internal hydrostatic pressure. 
However, the deviation of the theory from experi- 
ment at large strains is of such a character as to 
suggest that we have here to deal with flexible 
chains approaching their maximum extension, 
for which the Gaussian distribution does not hold. 
We shall then develop methods for discussing the 
stress-strain relations in flexible chains when the 
extensions are comparable to the maximum 
length, in Section 5. Section 6 deals specifically 
with chains consisting of completely independent 
links, and describes an improved model involving 
such chains. The problem of imperfectly flexible 
chains in general, and of chains with free rotation 
and fixed valence angle jn particular, is deferred 
to a later paper. Finally,’ in Section 7, a compari- 
son of theory with experiment is carried out, and 


JAMES AND E. GUTH 


further conclusions are drawn as to the conditions 
existing in soft rubber compounds. 


3. THE NEW MODEL FOR RUBBER: 
GAUSSIAN CHAINS 


We shall now discuss the stress-strain relation 
for our model of rubber, treating the material as 
essentially incompressible. This assumption is 
well substantiated by the experiments of Holt 
and McPherson” who have found that the vol- 
ume is essentially unchanged when soft gum 
compounds are given extensions of up to 300 
percent. 


We shall consider a unit cube of rubber 


stretched by forces of magnitude Z and having 
the z direction, and applied uniformly to cube - 
faces perpendicular to the z axis. The cube will 
be deformed into a parallelopiped of dimensions 
L,, Ly, L.z, with L,=Ly,. According to the as- 
sumption of incompressibility we have then 


L,L,L,=L2L,=1. (3.1) 


In the treatment of the model thus stretched 
the first task is to compute the pull of the irregu- 
lar molecular network. This may be thought of 
as a system of junction points connected by seg- 
ments of flexible chains which are independent 
of each other. The theoretical treatment of this 
network may be split into two parts: 

(a) Determination of the configuration function 
for a single chain segment, C(ra, fs). This func- 
tion of the positions of the endpoints of the chain 
will depend only on the separation rag of the ends. 
(b) Determination of the -configuration func- 
tion for the whole net when a number of the junc- 
tion points, specified by the vectors 1, fe - 
are kept fixed, as at the surface of the model. 
This configuration function will be a function 
P(r,, 1 -+-) of the positions of the fixed junction 
points, obtainable by integration of a product of 
the functions C(rz, fs) over the positions of all 
junctions which are not fixed. 

Problem (a) is solved in this paper for several 
types of flexible chain, and for all extensions. 
Problem (b) has been solved rigorously only 
when each chain segment can be treated as of 
Gaussian type. This assumption requires that the 


2W. L. Holt and A. T. McPherson, Rubber Chem. 
Tech. 10, 412 (1937). 
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Fic. 2. Comparison of experiment with theory using 
Gaussian chains. Solid line, experiment. Circles, theoretical 
values. Dashed line, asymptote approached by theoretical 
values at large extensions. 


segments of chains between junctions all contain 
very many links and (as will be apparent later) 
that these segments normally have extensions 
small compared with their maximum length. For 
the purposes of the present section it is sufficient 
for us to make use of results established in the 
Appendix: The forces exerted by a network of 
Gaussian segments, irregular in detail but homo- 
genous and isotropic in the average, with fixed 
points constrained to lie on the six planes x=0, 
x=L,, y=0, y=Ly, 2=0, z=L,, are the same as 
those exerted by a simplified network in which 
one set of Gaussian chains runs between the 
planes x=0, and x=L,, a second completely in- 
dependent set runs between the planes y=0 and 
y=L,, and a third between the planes z=0, and 
z=L,. Each set of chains can be taken to include 
M chains, each of WN links, and each chain may 
be thought of as parallel to an edge of the paral- 
lelopiped, in the sense that its fixed ends lie on 
a straight line parallel to that edge. 

The Appendix establishes this result by a treat- 
ment which gives a picture of the behavior of the 
network in both unstretched and distorted con- 
ditions and indicates a method of computing the 
forces from the detailed picture of the network. 
There is also possible a more direct proof, which 
consists essentially of the reduction of a Gaussian 
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distribution in many dimensions to a Gaussian 
distribution in a smaller number of dimensions 
by integrating over certain of the variables. 

Replacing the irregular network of our model 
by the equivalent system of three sets of M 
chains of N links acting independently in the 
three coordinate directions, we can determine the 
internal pressure in the material by a considera- 
tion of one of the free surfaces, for which the 
outward push of this pressure is balanced by the 
pull of M molecular chains, each extended to 
length L, and thus exerting an average force of 
CRT /(L?)\y |Lz. That is, 


PL,L.= M[RT/(L)w JL, 


P=M-[RT/(L)wL-". (3.2) 


Thus the internal pressure in bulk rubber de- 
creases when the material is subjected to stretch. 

To obtain the stress-strain relation for the 
cube, one needs only to consider similarly the 
equilibrium of one of the end faces, equating the 
outward push of the pressure plus the external 
force to the pull of the molecular chains. One 
finds at once 


Z=M[RTMLw[L:-1/L2]. (3.3) 


The first term on the right is the force given by 
the bundle model of rubber already mentioned, 
while the second term is the correction for the 
effect of the decreasing internal pressure as the 
material is stretched. 

Figure 2 shows a curve of the form given in 
Eq. (3.3), in comparison with an experimental 
curve.” It will be seen that the nature of the 
characteristic knee in the experimental curve is 
well reproduced by this simple theory, in which 
the only adjustable parameter governs the ver- 
tical scale of the figure. The form of this knee is 
determined by the decreasing internal pressure 
in the material. Figure 2 thus indicates immedi- 
ately the importance of the internal pressure in 
the understanding of the form of the stress-strain 
curve of rubber. We defer more detailed compari- 
son between theory and experiment to Section 7, 
after the theory has been further developed. 

This theory gives a relation of Z to L, which 
approaches linearity as L, becomes large and the 


1 R, L. Anthony, R. H. Caston, and E. Guth, J. Phys. 
Chem. 46, 826 (1942). 
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internal pressure becomes negligible in effect. 
The asymptote for large L, passes through the 
origin of the (Z,Z.) plane, a feature which is 
possessed by the tangents to many rubber stress- 
strain curves in the nearly linear portion for 
intermediate stresses. 

The differential Young’s modulus is, of course, 
not a constant, since Hooke’s law does not hold. 
Referred to the original cross section, we have 


i] 
Z 3 


For large values of LZ, this approaches 4 of its 
value for L,=1—another feature for which one 
can easily check in experimental stress-strain 
curves. Referred to the actual cross section by 
use of (4.1), the stress is given by 


kT 1 
rl ae | 
(L*) ny L, 


and the differential Young’s modulus is 
kT 
M- 
(L? hw 


dZ kT 
E(L, =—=M- E (3.4) 
dL 


(L?) wy 


f“—-/- (3.5) 


L,? 


[22+ (3.6) 


Unlike E(L,) this increases indefinitely with L,.. 
Introducing the perhaps more customary strain 
variable, the extension e=L,—1, one finds 


Z=3M(RT/(L?)w le(L1—e+$e—---]. (3.7) 


This form exhibits well the large deviations from 
Hooke’s law shown by rubber. 

The internal pressure in the unstressed rubber 
is simply related to the differential Young’s 
modulus for zero stress: 


Py= MRT /(L)y=4E(1). (3.8) 


For the rubber sample of Fig. 2, Po is thus 
about 5 kg/cm?, a value which can be taken as 
representative for soft gum compounds. When 
this pressure is reduced by stretching one might 
expect an increase in the volume of the material. 
Since the compressibility of rubber is about 10~* 
cm?/kg, the increase in volume due to complete 
removal of the internal pressure would be some 
0.05 percent, an amount quite compatible with 
the observations of Holt and McPherson. Gen- 
eralization of the theory to take account of the 
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compressibility of the material, while simply ac- 
complished, is, accordingly, of little interest. 


4. FURTHER APPLICATIONS OF NEW MODEL 
WITH GAUSSIAN CHAINS 


We now illustrate the application of our model 
to the treatment of a number of other properties 
of rubber. 

The correction for the volume thermal expan- 
sion of the material is of considerable interest. If 
one deals with a cube of rubber having unit vol- 
ume at temperature T> and a coefficient of ther- 
mal expansion a, the volume becomes at tem- 
perature 7 


V=i+a(T—T). (4.1) 


Generalizing (3.3) to apply to temperatures other 
than the original one, one finds 


kT 1+a(Z'— T) 
Z=M. —|L.— | (4.2) 


(Lm 


L,* 


the correction for the internal pressure being 
changed by a factor V. According to this, the 
stress at constant strain will not be exactly pro- 
portional to TJ, but will involve also a small 
negative term proportional to J’. The tendency, 
described by this quadratic term, of thermal 
expansion to decrease the stress in material 
heated under constant strain is of course the 
dominant tendency in the case of a stretched 
wire. It should be noted, however, that in or- 
dinary materials the temperature dependence 
of this term is different from that in rubber, 
the stress going down linearly with increasing 
temperature. 

The effect of this correction term appears most 
clearly when one plots stress against tempera- 
ture at various constant strains. Each stress- 
temperature curve will appear nearly straight if 
the range of temperature considered is not ex- 
cessive; the most striking effect of the new terms 
in (4.2) is the change in slope of these curves with 
changing strain. If the strain is small (0Z/dT) 
is negative, and the behavior is that of an ordi- 
nary body, while if the strain is large (0Z/07) 
has the large positive value usually considered 
characteristic of rubber. From (4.2) one finds 
(0Z/0T) rr, 


=[(Mk/(L*)wLL.—(4+aTo)/L2]. (4.3) 
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Thus the critical extension L,, at which the stress- 
temperature curve changes its slope is given by 
the isometric thermoelastic inversion point 


zc’ =1+aT>. (4.4) 


For a soft unaccelerated gum compound with 
8 percent sulphur content, a=6.2X10~*. For 
T =300°K one then obtains L,.=1.065—the in- 
version should take place for a strain of about 7 
percent. 

Since 


(0L./8T)z= —(0L./0Z)7(0Z/dT)1,, (4.5) 


(0L./dT)z and (0Z/dT);,, must change sign to- 
gether. Thus the thermoelastic inversion point is 
the same for isotonic curves as for isometrics. 
We now compute the linear thermal expansion 
coefficients of stretched rubber, both parallel and 
perpendicular to the direction of stress. The ex- 
pansion coefficient parallel to the stress is 


a, (L,) = (1/L,)(0L./0T)z. (4.6) 


Using (4.5) together with (4.2) we obtain for 
T= To 
1 L-1 
aun(L-) se 


+ . 
TL3+2 L342 





The expansion coefficient perpendicular to the 
stress is 

a,(L,)=(1/L,)(dL,/dT)z. (4.8) 
Since 


VuLL.=-1+e(T—-T2, (4.9) 


one has at T=T», 


117-1 aLl+i1 
a,(L,) ema 


~ . (4.10) 
2T LA+2 2L,3+2 





In the unstretched state, L,=1, we have 
a,(1) =a,(1) =a/3, (4.11) 


as for any material. But as soon as the rubber is 
appreciably stretched the terms involving a be- 
come negligible: 


a, (L.)2—(1/T)(L2—-1)/(LP +2), 
a, (L.)A(1/2T)(L3—1)/(L2+2). 


When these equations hold, the changes of di- 
mensions of the material are quite unrelated to 
its volume expansion; they are shifts in relative 


(4.12) 
(4.13) 
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dimensions needed to restore equilibrium as the 
effective stresses in the flexible chains vary with 
temperature. The coefficient of ‘‘expansion’”’ 
parallel to the stress is negative, while the other 
is positive, a striking illustration of the anisotropy 
of stretched rubber. These linear expansion co- 
efficients are large, of the order of magnitude 
1/T of the volume expansion coefficients for 
gases. Equation (4.7) will make it clear why rub- 
ber stretched by a weight will contract visibly 
when it is warmed. 

The change in temperature AT on adiabatic 
stretching is given according to thermodynamics 
by the following formula: 


aT=(T/C:) f (0Z/dT)tdL., (4.14) 


where C,, designates the specific heat at constant 
length. 

Inserting for (0Z/dT)z, its value from (4.3) 
and integrating we obtain 


AT =(KT/2C,)(L24+L.—2(i+aT») ] 


According to this formula AT starts with the 
value zero at L,=1, reaches a minimum for 
L.=L. as given by (4.4), then becomes zero for 
the adiabatic thermoelastic inversion point, de- 
termined by 


Lift+La- 2(1+aT») =0, 


(4.15) 


(4.16) 


and takes on positive values for L,>L.. For 
T=300°K and a=6.2X10~ one obtains L.o= 
1.13—-AT becomes zero for a strain of 13 percent, 
i.e., twice the strain at the isometric thermoelastic 
inversion point. 

This model permits consideration of other dis- 
tortions than simple stretch of rubber. We shall 
discuss briefly three of these: stretching in two 
directions, shear, and swelling followed by 
stretching. 

Let a piece of rubber, a unit cube in the un- 
stretched state, be stretched to length L, in the 
x direction, and thereafter stretched to length L, 
in the z direction. Neglecting compressibility 


effects, one has then 
L,=1/LeLs. (4.17) 


The internal pressure is found, as before, by con- 








466 H. M. JAMES 
sideration of the equilibrium of the free faces of 
the material, to be 


P=(MRT/{L*)w ](1/L2L). 


Thus by consideration of the faces perpendicular 
to the z axis, 


Z=(MRT/{L*)y |[L.—1/L2L 2 |. 


To make a proper check on this expression care 
must be taken that LZ, and L, are not made ex- 
cessively large, for the deviations from Gaussian 
behavior of the molecular chains in real rubber, 
which will be discussed in later portions of this 
paper, will be particularly important in this case. 
It is a familiar fact that rubber stretched strongly 
in one direction is difficult to stretch in a perpen- 
dicular direction. As ordinarily observed this be- 
havior is, however, probably largely conditioned 
by crystallization induced in the stretched rubber. 
Examination of Section 3 and of the Appendix 
will show that the reduction from the irregular 
network to the model with parallel molecular 
chains will apply also for the discussion of shears. 
If a shear slides the xy planes of the model over 
each other in the x direction, both x chains and 
y chains will have their extensions unchanged 
and will exert no restoring forces. The internal 
pressure will also give no restoring forces. The 2 
chains, however, will be stretched and will exert 
a component of force opposing the shear. Con- 
sidering, as before, a unit cube, now sheared 
through an angle @, the mean force exerted by a z 
chain, stretched in length by a factor 1/cos @, 
will be [k7/(L?)w |[1/cos 6], and the total com- 
ponent of force opposing the shear will be!* 


(4.18) 


(4.19) 


F,= MURT/{L*)w | tan @. (4.20) 
The shear modulus is thus 
G=dF,/d0= MLRT/{L*),, | sec? 6, (4.21) 
or, for infinitesimal deformations, 
G= M[RT/{L*)w J. (4.22) 


“Jt will be noted that the shearing force is strictly 
proportional to the shear displacement, this being equal 
to the thickness of the sample times tan 6. M. Mooney, 
J. App. Phys. 11, 582 (1940), has discussed from the point 
of view of phenomenologic elasticity theory the form which 
may be taken by stress-strain curves of elastic incom- 
pressible bodies having this property, and our relation is 
of course to be found as a special case of his general result, 
as pointed out by Dr. Mooney to us some time ago. The 
dependency upon temperature, however, cannot be ob- 
tained from elasticity theory. 
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This, of course, satisfies the relation with 
Young’s modulus, 


G=E/2(1+u), 


as found in the theory of elasticity for infinitesi- 
mal deformations, when Poisson’s modulus yu has 
the value 0.5 valid for incompressible materials. 

Rubber has the remarkable property of being 
able to absorb large quantities of liquid hydro- 
carbons without loss of shape or elastic proper- 
ties. The behavior of the resulting swollen rubber 
suggests that the hydrocarbon molecules do not 
greatly affect the structure of the material, but 
merely serve to increase its bulk and tend to 
prevent crystallization. If this assumption is cor- 
rect our model should be able to explain the 
elastic properties of swollen rubber. 

Suppose a unit cube of rubber to be swollen by 
the addition of an inert hydrocarbon until its 
linear dimensions are increased by the swelling 
factor o. If it is then stretched or compressed 
until its zg dimension is L,, one has 


L,L/=0'. 


(4.23) 


(4.24) 
Proceeding as before, one finds 
P= M[RT/(L?)w \[1/Lz], (4.25) 


which is identical with (3.2). It is to be noted, 
however, that in the absence of external forces 


Po= MURT/(L?)w ](1/o). (4.26) 


Comparing with (3.8), one sees that, as rubber is 
swelled by the introduction of inert liquid mole- 
cules, the internal pressure falls, varying in- 
versely with the linear dimensions of the material. 
The stress-strain relation for the swollen cube is 
found to be 


Z=MURT/L)w][L:—0°/L2]. (4.27) 


This equation is not entirely comparable to (3.3), 
which relates to a portion of material which, in 
the unstretched state, is a unit cube. Modifying 
(4.27) to relate to a unit cube of the swollen 
material, rather than a swollen unit cube, we have 


Z=M(RT/(L*)w](1/o)[L2—-1/L2]. (4.28) 


That is, swelling of rubber should reduce its 
rigidity by a factor of 1/o, but leave the form of 
the stress-strain curve unchanged. 
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Results (4.26) and (4.28) are simply under- 
stood as follows: Swelling of the material di- 
minishes the number of chains crossing unit area 
by the factor 1/o?, and also the number of links 
per unit stretched length of chain by the factor 
of 1/a. The latter change increases the force per 
chain by the factor o. Thus, finally, the forces 
per unit area exerted by the chains, as well as 
the pressure forces which resist them, will be 
changed by (1/0*)o=1/c. 

The decrease in the rigidity of swollen rubber 
described by (4.28) is well known to occur. Un- 
fortunately there does not seem to be available 
any experimental material permitting satisfac- 
tory quantitative comparison with this theory. 
Tiltman and Porritt!® have measured the effect 
of benzene on the elastic properties of one par- 
ticular vulcanizate, but the high percentage of 
zinc oxide in their material makes their data un- 
satisfactory for our purposes. It may be noted 
that, at very high strains, they observe decreases 
of stress on the addition of relatively small 
amounts of benzene which are quite out of pro- 
portion to those indicated by this theory. Here 
the effect is doubtless that of a reduction of the 
strain-induced crystallization. 

Now it is obvious that the results of Sections 3 
and 4 can be valid only when the extensions of 
the chains in the molecular network are not too 
large compared with their total lengths, since 
only in this case can the assumed Gaussian con- 
figuration distribution be an adequate approxi- 
mation. The Gaussian distribution assigns a non- 
zero probability to any extension of the chain, 
however great, whereas the actual function for 
any finite chain must become zero as soon as L 
exceeds the total length of the chain. As L ap- 
proaches this limit the true distribution function 
must fall toward zero faster than the Gaussian 
function. Then —log C(L) must increase more 
rapidly than quadratically, and Z(L) for a single 
chain must increase more rapidly than linearly, 
becoming infinite as ZL approaches Lax. This 
deviation from the linear stress-strain relation 
will make itself felt when bulk rubber is subjected 
to large strains; the effect will be the introduction 
of an upward curvature into the stress-strain 


4 A. H. Tiltman and B. D. Porritt, India-Rubber J. 78, 
345 (1929). 


curve, the stress increasing indefinitely as the 
stretched length approaches the finite maximum 
extension of the network. It is thus evident that 
with a more careful treatment of the stress-strain 
relation in molecular chains our model will show 
an upward curvature in the stress-strain curve 
similar to that shown in the experimental curve 
of Fig. 2, but lacking in the theoretical results 
of the present section. (The observed upward 
curvature is not to be explained entirely as due 
to crystallization effects, as it is evident even for 
strains insufficient to induce crystallization, or 
when this is suppressed by swelling.) 

We turn, then, to a consideration of the stress- 
strain relation in a flexible thermally agitated 
chain when its extension is comparable to its 
total length. ; 


5. STATISTICAL TREATMENT OF HIGHLY 
EXTENDED FLEXIBLE CHAINS: 
GENERAL METHODS 


In applying the equation 
Z=—kTd{log C(L, N) |/dL (5.1) 


to the computation of the stress-strain relation 
for a flexible chain under thermal agitation, one 
is, essentially, considering the chain to have a 
fixed extension and computing the average force 
exerted by it. That is, one determines 


Z=2Z(L). (5.2) 


An alternative approach which offers marked ad- 
vantages for our problem is that of considering 
the chain with constant forces exerted on its ends 
and computing the average extension: 


L=L(Z). (5.3) 


For chains of few links, the instantaneous value 
of L may differ markedly from L, and the same 
is true for Z and Z. Under these circumstances, 
Eqs. (5.2) and (5.3) are by no means converse 
relations; only in the limit of very long chains 
will the methods lead to essentially identical 
results. The advantage of the approach through 
(5.3) will be made evident by the example at the 
end of the present section. 

We consider first a one-dimensional chain. If 
a constant force Z be applied to the ends of the 
chain there will be a potential energy —ZL asso- 
ciated with the extension ZL. Let the number of 
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configurations available to the chain when its 
extension is L to within dL be 


dC=C(L)daL. (5.4) 


- The probability that the chain acted on by the 
force Z will have the extension Z to within dL 
will then be proportional to 


dP =et2ZLlkTGC = C(L)etZ“/*TdL. (5.5) 
Let 
f=Z/kT. (5.6) 


Then 


Lf) = f LC(L)e/“d. / f C(L)eftdL, (5.7) 


L(f) =d{ log f C(L)e!at jf: df. (5.8) 


If Z can take on only discrete values we have, 
instead, 


L=d(log ¥1 C(L)e’“) /df. (5.8’) 
Let G(x) be the generating function for C(L). 
That is, for the continuous case 


G(x) = f Cera, (5.9) 


and for the discrete case 


G(x) = Dir C(L)x*. (5.10) 


Then in either case!® 


L =d[log Ge’) ]/df. (5.11) 


Often the most convenient method of determin- 
ing C(L) is by use of the generating function. 
Equation (5.11) offers a method of passing from 
G(x) directly to a relation of the form (5.3), thus 
obviating the necessity of determining C(L) 
and applying (5.1) to obtain a relation of the 
form (5.2). 

Another point of view is often useful. We have 


lefty, a f eLC(L)aL / f C(L)dL. (5.12) 


16 The discrete case can be formally unified with the 
continuous case by taking the integrals in the Stieltjes 
sense. 
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Thus a 
L=d (log (e")w)/df. 


As a third approach, we expand the exponen- 
tial in the integrand, Eq. (5.12). Then 


(5.13) 


(2), =, ft/r! f LC(L)dL 4 f C(L)dL. (5.14) 
Let 
(L\y= { L'C(L)aL f f C(L)dL, (5.15) 


the average of L’ in the force-free case. Using 
these relations, (5.13) becomes 


L=d[log Ys (L*)wf*/7!]/df. 


It is instructive to relate the statistical treat- 
ments of flexible chains at constant extension, 
Z=Z(L,T), and at constant stress, L=L(Z, T), 
with the corresponding thermodynamic ap- 
proaches to the problem. 

When temperature and extension are chosen 
as independent variables in the treatment of a 
thermodynamic system it is natural to introduce 
the Helmholtz free energy function F(L, T) and 
to determine stress from this: 


Z=(aF/aL)r. 


(5.16) 


(5.17) 


As we have already noted, the free energy is 
related to the number of configurations accessible 
to the system, C(L). (C(L) is the partition 
function.) 


F(L, T) = —kT log C(L). (5.18) 


Thus Eq. (5.1) is simply the statistical form 
of (5.17). 

Now if temperature and stress are chosen as 
independent variables in treating a thermo- 
dynamic system, one introduces the Gibbs func- 
tion, which, for our system, becomes 


@(Z, T)=U-ZL-TS. (5.19) 


From this one determines the extension of the 
system, 


L= —(d@/0Z) r. 


Statistically the Gibbs function is related to the 
sum over all accessible states of the statistical 
assembly, weighted with the factor e?”/*?, In 
analogy to (5.18) we have 


@(Z, T)=—kT log D(Z), 


(5.20) 


(5.21) 
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where D(Z) is a generalized partition function 


D(Z) = f C(L)e24/*TdL, (5:22) 


which is essentially the Laplace transform of the 
partition function C(L). Since from (5.9) we have 


D(Z) =G(e?!*?) =G(e/), (5.23) 


it is evident that (5.11) is the statistical version 
of the thermodynamic relation (5.20). 

That Z=Z(L) and L=L(Z) are converse rela- 
tions only in the limit of very long chains be- 
comes, in thermodynamic terms, the statement 
that only in this limit does there exist an equation 
of state 


Z=Z(L, T) or L=L(Z, T) 


describing the system. 

Depending on the circumstances, each of 
(5.11), (5.13), and (5.16) may serve as useful 
alternatives to (5.1). 

The extension of these results to the case of a 
chain in three dimensions involves no difficulty. 
Let one end of the chain be fixed at the origin. 
Then the number of configurations of the chain 
for which its free end will lie in the infinitesimal 
volume dV will be 


dC=C(x, y, 2)dV. (5.24) 


If there is no preferred orientation of the system 
C will be a function of x?+y?+2? only. 

Let a constant force of components Z,;, Z,, Zz 
be applied to the free end of the chain. Then the 
probability of finding the end of the chain in 
volume dV will be proportional to 


dP=exp {(xZ,+yZ,+2Z.)/kT}C(x, y, 2)dV. 


(5.25) 
As before, letting 


f.=Z./kT, 
fy=4Zy/kT, 


(5.26) 


et.. 
one has 


0 
sa as f elettvt9C(x, y,2)dV, (5.27) 
Of 


dj 
die, f eletulytt C(x, y,2)dV. (5.28) 
Ofy 
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Ordinarily there will be no preferred orienta- 
tion of the chain in space. Then, if a force is ap- 


plied in the x direction, one will have 7=Z=0. 
Defining 


C,(x) =f fice, y, 2)dydz, (5.29) 


one finds 


(5.30) 


d 
~=—lo ferc. x)dx. 
i g ( 


Thus, in this case, one can reduce the three- 
dimensional problem to a one-dimensional prob- 
lem by projecting the configuration function onto 
the direction of the applied force. In particular, 
if the chain consists of independent links, the 
problem can be reduced to a one-dimensional 
case at the very start, by projecting the proba- 
bility distribution for the individual links onto 
this direction, and employing this projection in 
the treatment of the chain as a whole. If the 
links of the chain are not independent—i.e., the 
probability distribution for each link depends on 
the orientation of the adjacent links—the pro- 
jection method can still be applied. One has then 
to deal with probability distributions for the 
projections of link extensions onto the force di- 
rection, these depending on the values of the 
projections for adjacent links. Examples of such 
reductions will be found in succeeding sections. 
As a simple illustration of the application of 
these methods, we consider a one-dimensional 
chain of N links, each of which will, in the ab- 
sence of external forces, contribute to the exten- 
sion of the chain an amount +1 or —1, with 
equal probability. The relative number of ways 
in which an extension L of the chain can be 
obtained is clearly the coefficient of x” in 
(1/2) (x+1/x)*. We have, thus, 
G(x) =(1/2%)(x+1/x)”, (5.31) 
and 
N! 
2% [(N+Z)/2]![(N-—L)/2]! 


To apply Eq. (5.1) one inust first use Stirling’s 
approximation, V!2(27N)!N%e-*, for the fac- 
torials. In the case of the factorial [(N—L)/2]! 





C(L, N) (5.32 
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this will not be legitimate for LZ arbitrarily close 
to N, but, by increasing N, this approximation 
can be made good for values of the fractional 
extension 


$= L/Lesx=L/N (5.33) 


arbitrarily close to 1. In terms of the fractional 
extension one has 


C(L, N) = 1+4)-/D a+) 
(nny ) 


X(1—-#)--9(1—2)-F (5.34) 


and 


kT 
Z=— log 
2 





I+t 2 ¢ 1 
( 
N 1-#? 


ar —). (5.35) 


To apply (5.11) we have simply to note that 
according to (5.31) 


G(e’) =[coshf |. (5.36) 


Then 


(5.37) 


” d 
L= ‘iy log coshf= WN tanhf, 


t=tanhf (5.38) 


exactly. Equation (5.38) is the inverse of (5.35) 
only in the limit of WV large, in which all terms on 
the right, except the first, are to be neglected. Its 
derivation is almost trivial, there being no terms 
of order 1/N to be separated out and dropped in 
passing to the limit. This illustrates the essential 
advantage in computing L(Z) instead of Z(L): it 
is the former quantity which varies the more 
simply with the number of links in the chain. In 
less simple problems the importance of this ad- 
vantage is greatly increased, as it becomes much 
more difficult to construct C(L, N) explicitly. 
Usually it will be necessary to derive C(L, N) by 
an inverse Laplace transformation from G(e’), 
whereas, by use of (5.11), the desired stress- 
strain relation can be derived directly from G(e’). 

The next section, and especially a later paper, 
will provide more elaborate illustrations of these 
methods. 


6. EXTENDED MODEL OF BULK RUBBER: 
CHAINS WITH INDEPENDENT LINKS 


The simplest type of flexible chain which one 
may use to represent a rubber molecule consists 
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of links of fixed length /, each of which will tend 
to take on all orientations in space with equal 
prabability, whatever may be the orientations of 
the neighboring links. In this model a link will 
correspond, not to a single carbon-carbon bond, 
but to a section of the chain so long that con- 
straints on the relative orientations of the ends 
are negligible. Just what this length should be is 
not well defined, and there will be a correspond- 
ing uncertainty in the magnitude of the forces 
computed using this model. In attributing a fixed 
length to each link we are also making a simplify- 
ing assumption which will somewhat affect the 
form of the stress-strain curve. Nevertheless, 
this model constitutes so marked an advance 
over that with Gaussian chains as to deserve 
consideration. In a later paper we shall compare 
this model with a more refined, but more in- 
tractable, model. 

A three-dimensional vector of magnitude / 
which takes on all orientations with equal proba- 
bility will have all x components between —/ and 
+/ with equal probability. Our three-dimensional 
chain is thus equivalent to a one-dimensional 
chain in which the ith link may make any contri- 
bution x between +/ and —/ to the extension, 
with equal probability when external forces are 
absent. 

To treat this system one may use Eq. (5.13). 
Note that 


(e7h) ny =(exp f(xitx2+:++))w 
= [Ii (exp fxibu, 
since the x’s vary independently. Now 


1 pt 
(exp fiu=— f exp fxidx,;=sinh Lf/Lf. 


Thus for a chain of JN links 
(ef”) = {sinh Lf/1f}%. 
Applying (5.13), we have at once 
t=coth /f—1/lf=L(/f), 

where £(x) is the Langevin function. 

It is not possible to pass rigorously from the 
stress-strain relation thus derived for a single 
chain to that for an arbitrary three-dimensional 


network, since the reduction of the network to an 
equivalent set of independent chains as carried 


(6.3) 


(6.4) 





ELASTIC PROPERTIES OF RUBBER 


out in Section 3 is not exact for these non- 
Gaussian chains. So long, however, as the condi- 
tions of stretch of the network are not too extreme 
we will expect that reduction to lead to a satis- 
factory approximation. Without further justifica- 
tion of the procedure we pass to the consideration 
of a model of rubber similar to that previously 
considered, except that in the unstretched state 
of the material each independent chain is taken 
to have an extension equal to x times its max- 
imum length. The maximum stretch for the 
model is thus given by a factor of 1/x, a quantity 
which will vary with the sample of rubber and 
will, accordingly, enter our treatment as a vari- 
able parameter. 

As before, we consider a unit cube of the ma- 
terial, taken, in a first approximation, as incom- 
pressible. In the stretched condition L,L,L.=1, 
as before. If a chain is given an extension L its 
fractional extension will be t=Lx. The stress- 
strain relation for the individual chain may then 
be expressed as 


Le=S&(1Z/kT). (6.5) 


This defines a relation between the fractional ex- 
tension of the model in a given direction and the 
tension of the chains extended in that direction 
which we express as 


IZ/kT =£-\(t). (6.6) 


One then finds, for the internal pressure in 
the model stretched only in the z direction, 
(Lz=L,=L-'), 


P=(MkT/IL,L)£“(t,). (6.7) 


The stress-strain relation for the model as a whole 
becomes 


Z=(MRT/I)[£L-"(t.) —(Ly/LL"(ty) ] (6.8) 


or 
Z=(MkT/I)[2-\(L-x) —L12-(L-'*)]. (6.9) 


This relation involves essentially two parameters; 
Mk/l, which determines the scale of the forces 
just as in the Gaussian case, and the new max- 
imum stretch parameter x. 

If the argument of the second term in the 
stress, that due to the internal pressure, is small, 


(6.10) 
(6.11) 


L(t) 238, 
Z2(MkT/I)[2-(Lx) —3x/L 2]. 
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When this approximation can be applied to both 
terms—all extensions small compared to the 
maximum—one has simply 


Z2(3MkkT/l)(L.—1/L 2], 


as for a Gaussian model. Putting x=1/N/ the 
relation (6.12) reduces exactly to (3.3), remem- 
bering that (L?)y=3NP. 

It should be noted that predictions of the 
theory for shear, stretch in two directions, swell- 
ing, and so on, as given in Section 4 for Gaussian 
chains will also be modified for higher extensions. 
Since these modifications are not formally inter- 
esting, and are quite insignificant so long as one 
does not deal with extreme distortions of the 
material, we shall not discuss them at this point. 

Comparison of these results with experiment 
is deferred to Section 7. 

The appearance of the Langevin function in 
these results is most easily understood by con- 
sideration of a method of treating stress-strain 
relations which is applicable only to chains of 
independent links. In cases in which the contri- 
bution of each link to the extension is inde- 
pendent of that of every other link, the reduction 
of the total problem to a one-dimensional equiva- 
lent may be accomplished by making a corre- 
sponding reduction for each individual link. We 
consider then a chain of N identical independent 
links, each of which will, in the force-free case, 
contribute to the x component of the extension 
of the chain according to the probability distri- 
bution c(x). Then 


(e/%2)w =L(e*)n J", 


the average on the right being computed for an 
individual link. Consequently, 


L,= Nd (log (e*\w)/df, 


(6.12) 


(6.13) 


(6.14) 


or 


t= [x exp («Z./kT)c(x)dx 


/ f exp (xZ2/kT)c(x)dx. (6.15) 


Here we do not deal with ¢, as before, since the 
maximum extension might be infinite, but rather 
with Z, the average contribution of each link to 
the extension. 
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Equation (6.15) states that the contribution of 
each independent link in a flexible chain can be 
computed as if it were the only link present and 
the external forces were applied directly to its 
ends. This is true whether or not all links in the 
chain are described by the same distribution 
function c(x). 

The situation of an independent link in an 
extended flexible chain is thus completely anal- 
ogous to that of a polar molecule in an electric 
or a magnetic field. The thermal agitation of the 
chain tends to disorient the links, to make the 
average contribution to the extension vanish, 
while the tension of the chain, in effect, is equiva- 
lent to a uniform field of force tending to align 
the links. The computations of the average con- 
tribution of a link to the extension and of the 
average contribution of a polar molecule to, say, 
the magnetic moment of a material are thus 
formally identical. If the link has a fixed length, 
the molecule a fixed moment, the computations 
are identical in detail; the Langevin function 
appears in each result. 

In the same way the function tanh f which 
appears in Eq. (5.38) for the extension of a chain 
of links for which only two configurations are 
possible is the function which appears in Ehren- 
fest’s classical theory of magnetism, which as- 
sumed the existence in crystals of polar atoms 
having two orientations of equal potential energy 
in the absence of a magnetic field. A chain model 
in which each link could take on a discrete set of 
orientations would similarly lead to the Brillouin 
function, of which Ehrenfest’s tanhx is the 
simplest special case.” 

We conclude this section with a more explicit 
statement of the role of the Gaussian distribution 
as an approximation in the case of chains of inde- 
pendent links. We assume that c(x) is an even 
function of x, as it must be if there is no preferred 
direction in space. Then, expanding the ex- 


17 P, Ehrenfest, Proc. Amst. Akad. Sci., Dec. 18, 1920 
or Leiden Communications, Suppl. 44b; cf. also P. Debye 
in Marx, Handbuch der Radiologie, Vol. VI, p. 712. For 
the Brillouin function, L. Brillouin, J. de phys. ev rad. 8, 
74 (1927); also P. Debye, loc. cit., p. 713. More recent 
accounts are given by J. H. van Vleck, The Theory of 
Electric and Magnetic Susceptibilities (Oxford, 1932), pp. 
257, 289; and R. H. Fowler and E. A. Guggenheim, Sta- 
tistical Thermodynamics (Cambridge, 1939), pp. 629-630. 


JAMES AND E. GUTH 


ponentials of Eq. (6.15), one has 


2 i 34-4 rr 
t= f |p + i" +++ lote)de 


1 
/f| : “s fea, (6.16) 


t= fen + SEH — H(0Ym) Eo 


or 
(6.17) 


For very small f and ¢, the relation between these 
quantities is determined by the value of (x*),, 
and is otherwise independent of the form of the 
distribution function. For larger extensions there 
will be deviations from the linear relation be- 
tween ¢ and f which will be determined by the 
relative values of (x*)w, (x*)w, -*+, for a single 
link. If these are all in the same ratio as for a 
Gaussian chain, the higher terms in (6.17) all 
vanish, and t=(x*)wf always. Otherwise, there is 
a deviation from the linear relation of f and ¢ 
which is quite unaffected by any increase in the 
number of links in the chain. That is, as the 
number of links in the chain increases, the linear 
relation of Z and L holds to greater and greater 
values of L, but not of L/N; deviations from the 
linear relation of Z and Z depend only on the 
fractional extension. Correspondingly, an in- 
crease in the number of links will make C(L) ap- 
proach the Gaussian form for any given range of 
L, but, in general, only for Z small as compared 
with the maximum extension of the chain. Simi- 
larly, by making the chain long enough, one can 
make the (L’)y of (5.15) approach the values for 
a Gaussian chain, for any given range of 7. There 
will, however, always be deviations of these quan- 
tities from the Gaussian values at larger 7, and 
these deviations will be important if the frac- 
tional extension of the chain becomes large. In 
short, the Gaussian approximation for the dis- 
tribution function is dependable only if the frac- 
tional extension of the chain is sufficiently small 
—smaller, as our results show, than extensions 
actually encountered in bulk rubber. 

The nature of the Gaussian form as an ap- 
proximation to the more exact form of the con- 
figuration function may be illustrated more ex- 
plicitly by reference to the one-dimensional prob- 
lem treated in Section 5. Neglecting a slowly 
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varying factor (1—#)~! in (5.34), we have 


C(L, N)= (1+4)-@™a+0(4 —~)- CVI) A-9, 
' (24rN)? ( ) 
(6.18) 


This can be written without further approxima- 
tion as 


C(L, N)= 


N 1+: 
(2eN)! : 


exp | —— | log ——dt 
2 1—t 


(6.19) 
Since 
log [(1+4)/(1-#)] 
=AUt+P/3+0/5+---), 


it is evident that, for small enough ¢, where all 
terms but the first are negligible, we have essen- 
tially the Gaussian distribution. Deviations from 
the Gaussian law set in when 


(6.20) 


(w/2) f aeeae = Nt*/4 


becomes appreciably different from zero—for 
greater absolute extensions Nt and smaller rela- 
tive extensions ¢ as N is increased—and are very 
great when ¢ is comparable to 1. 


7. COMPARISON OF THEORY WITH EXPERIMENT 


In explanation of the method to be used in 
comparing this theory with experiment we may 
survey briefly the approximations which have 
been made. A rather elaborate attempt has been 
made to treat the changes in entropy of the mo- 
lecular network when the bulk material is de- 
formed, as this accounts for much the greater 
part of the rigidity of the material. Here we were 
forced to neglect steric hindrances between chains 
—an approximation which we expect to be sig- 
nificant only for compact networks (highly vul- 
canized materials) or extreme deformations—and 
to carry out the reduction from an irregular to a 
regular one in an approximate way which will be 
increasingly in error as the deformation is in- 
creased. By including thermal expansion in our 
model phenomenologically we have introduced 
terms in the entropy which occur also in ordinary 
liquids, and also some terms in the internal 
energy. On the whole, however, the effects of 
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changes in the internal energy have been ignored. 
We have neglected van der Waals forces between 
the atoms in a single molecule, which may affect 
the behavior of the material for small extensions 
by changing both the internal energy and the 
entropy, and also van der Waals forces between 
essentially independent chains of atoms. At the 
smaller extensions we may expect these latter 
forces to be largely taken into account by inclu- 
sion of the thermal expansion, but at large exten- 
sions they give rise to crystallization, a large 
effect which is entirely outside the domain of the 
present theory. On the whole it is apparent that 
the theory is essentially one of entropy forces, 
valid at most for extensions which do not induce 
crystallization. It will be reasonable and desirable 
to resolve both theoretical and experimental re- 
sults so as to separate forces due to entropy 
changes from those due to changes in internal 
energy, and to deal with these separately, antici- 
pating the better agreement for the entropy 
forces. 

The division of the total force in unilaterally 
stretched materials into parts associated with the 
entropy and internal energy is expressed in Eq. 


(1.2). Since 
dZ 
-1(—) =Z-—Zs, 
OT/ 1 


0U 
tan(22) 22 
OL Ty 


the two terms in Z can be determined separately 
as soon as Z(L,7) is given, whether experi- 
mentally or theoretically. Proceeding from Eq. 
(4.2) and from the analogue of Eq. (6.11) which 
takes account of thermal expansion, for the cases 
of Gaussian and independent link chains, respec- 
tively, one finds for the entropy forces 


Zs=KT{L.—[1+a(2T- To) |L.7?} ’ 


Zs=K'T{\£“(L;k) 
—3«(1+a(2T—T») ]L;°}, 


(7.1) 


(7.2) 


(7.3) 


respectively. The forces associated with the in- 
ternal energy are 


Zu=KaT°L,’, 
Zu=3K'«aT°L,, 


(7.4) 
(7.5) 


in the two approximations. In both cases the 
thermoelastic inversion point is given by 


L,.=t+aTo. (7.6) 
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Fic. 3. Comparison of theory and experiment 
for natural rubber A. 


Experimental stress-strain curves obtained 
under conditions carefully enough controlled to 
permit satisfactory comparison with our theo- 
retical results are not abundant in the literature. 
In our preliminary report‘ comparisons were 
made between (7.2), with a=0, and total stress- 
strain curves from the literature, with very en- 
couraging results in the range of validity of that 
simplest theory. However, all the older experi- 
mental results are very incomplete, and, in par- 
ticular, they do not permit a resolution of the 
observed stresses into parts Zs and Zy. To make 
possible a more satisfactory comparison with 
experiment some very careful and complete ob- 
servations have been made and have been de- 
scribed elsewhere, by our experimental associates 
at Notre Dame, Dr. R. L. Anthony, Dr. R. H. 
Caston, Dr. L. E. Peterson, and Mr. S. L. Dart. 
We shall describe here the results with two types 
of vulcanized rubber and one synthetic rubberlike 
material of widely different composition. 


A. Unaccelerated Natural Pure Gum" 


This material was an unaccelerated pure gum 
stock containing 8 parts of sulfur and 100 parts 
of pale crepe rubber, cured for 3 hours at 147°C. 


AND E. GUTH 


It shows little permanent set and no crystalliza- 
tion for extensions up to 350 percent at room 
temperature. At higher temperatures greater ex- 
tensions are possible without crystallization. 

For this compound a=6.2X10-*. The experi- 
mental value of L., is approximately 1.10, as 
compared with 1.065 as given by (7.6). 

Zs was found to be proportional to T, to within 
the experimental error. 

Figure 3 shows the experimental values for 
both Zs and Zy in comparison with the values 
predicted by Eqs. (7.3), (7.5) for a model with 
chains of independent links, with K’=24, 
k=0.12. 


B. Accelerated Natural Pure Gum" 


This stock was an accelerated pure gum com- 
pound containing 5 parts sulfur, 3 parts diphenyl- 
guanidine, 1 part zinc oxide, 100 parts pale crepe. 
It was cured for 30 minutes at 147°C. This com- 
pound showed little permanent set. Crystalliza- 
tion began at 250 percent extension at room tem- 
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Fic. 4. Comparison of theory and experiment 
for natural rubber B. 


18 R, H. Caston, Notre Dame Thesis. 

R. H. Caston, R. L. Anthony, and E. Guth, Detroit 
meeting of the American Physical Society, February, 
1942; cf. Phys. Rev. 60, 542 (1942). 
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perature, and at somewhat higher extensions 
when the temperature was raised. 

For this compound, a=7X10~‘. This leads to 
a calculated value L..=1.07, against the value 
1.14 given in reference 18. More recent work at 
Notre Dame, still unpublished, gives an improved 
experimental value L..=1.09, in better agree- 
ment with the theory. 

Zs was found to be proportional to 7, to within 
the experimental error. 

Figure 4 shows experimental values for both 
Zs and Zy as functions of ZL, measured at 20°C. 
The indicated theoretical values were computed 
by Eq. (7.3), with K’=30, x=0.143. In this 
case the fit was made in the middle of the range 
of L, since the discrepancy between the observed 
and computed thermoelastic inversion points in- 
dicated the presence of either errors or additional 
entropy terms in the forces measured at low 
extensions. 

Figure 5 compares experimental!® and theo- 
retical values for the temperature rise in this 
material on stretching. These experiments give 
the isometric inversion point (AT a minimum) 
at L=1.07 and the adiabatic inversion point 
(AT=0) at L=1.15, in excellent agreement with 
the theoretical values of 1.065 and 1.13, respec- 
tively. Theoretical and experimental values of 
AT agree well over the whole range considered. 
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Fic. 5. Temperature changes of material B on stretching. 
+: Value measured after extension. O: Value measured 
after retraction. X: Theoretical value. 





19 Experimental work by S. L. Dart (unpublished). 
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Fic. 6. Comparison of theory and experiment 
for Hycar compound C. 


C. Accelerated Pure Gum Hycar” 


The composition of this stock was as follows: 
1 part sulfur, 1 part crude lauric acid, 5 parts 
zinc oxide, 100 parts Hycar OR. This compound, 
like all synthetics, exhibits considerable perma- 
nent set if stretched to higher extensions. Curing 
is not enough to restrict plasticity; it is necessary 
in addition to use reinforcing agents such as 
carbon black. We have wished to avoid consider- 
ation of compounds in which the reinforcing ma- 
terial adds a new element to the behavior of the 
material, and have accordingly been forced here 
to deal with extensions not exceeding 130 percent. 
Within this range the Gaussian approximation 
equations (7.2) and (7.4), are valid. No crystalli- 
zation whatever was observed on stretching this 
material. 

From the observed a=5X10~ one computes 
L.c=1.05, as compared with the observed 
L.z-=1.10. 

Zs was found to be proportional to 7, to within 
the experimental error. 

Figure 6 compares experimental and theo- 





20 LL, E. Peterson, Thesis, Notre Dame, 1942; cf. L. E. 
Peterson, R. L. Anthony, and E. Guth, Ind. Eng. Chem. 
34, 1349 (1942). 
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retical results for this material. Again the fit was 
made with reference to moderate rather than 
very small extensions of the material. No esti- 
mate of « is possible. 


DISCUSSION 


The proportionality of Zs to T at constant L, 
shown by all these samples, is in complete accord 
with the theory. The terms in 7? appearing in 
(7.3) are too small to produce an appreciable 
disturbance of this proportionality over the tem- 
perature range investigated. 

The observed entropy stress-strain curves are 
in very satisfactory agreement with the theory 
in all cases, over the whole range of extensions 
here considered. At the very highest extensions 
of compounds A and B the observed stress is 
beginning to increase more rapidly than the 
theoretical. This may be considered as an effect 
of incipient crystallization, or it may be due to 
the approximate replacement of an irregular net- 
work by bundles of chains all with the same 
fractional extension in the unstretched state of 
the material. A more satisfactory approximation 
would certainly require the use of chains with a 
variety of initial fractional extensions. This would 
lead to just the behavior here observed—a more 
rapid increase in stress at large extensions than 
one would expect from the behavior at moderate 
extensions, due to the increasing relative impor- 
tance of the shortest chains. 

Compound B shows in addition an appreciable 
discrepancy in Zs for L near 1, associated with a 
nearly equal and opposite discrepancy in Zv. 
Compound C shows a similar discrepancy which 
vanishes for L = 1.2—that is, rather more rapidly 
than in the case of compound B. The variability 
of this discrepancy and its rapid decrease with 
decreasing internal pressure in the material sug- 
gests that it is due to van der Waals forces ne- 
glected in the theory. The close agreement of the 
theoretical and experimental thermoelastic in- 
version points indicates that the actual discre- 
pancy is even smaller than it appears to be in 
Figs. 3, 4, and 6. 

Aside from the small effect just mentioned, 
synthetic rubber C gives excellent agreement 
with the theory. It is thus proved, apparently for 
the first time, that synthetics behave quantita- 
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tively in the same way as natural rubber under 
conditions for which thermoplasticity is negli- 
gible. It is to be emphasized that rubberlike 
elasticity requires the presence of a coherent 
network of flexible molecular chains, but does 
not depend on the special nature of the chains. 
The form of the stress-strain curves for small 
extensions (at least in the absence of excessive 
van der Waals forces) is then universal; varia- 
bility in the behavior enters when plasticity is 
not completely suppressed, or the network ap- 
proaches its maximum extension. 

Experimentally, the entropy forces account for 
80 percent or more of the total stress. As was to 
be expected, our theory gives a less satisfactory 
account of the smaller internal energy force. For 
very small extensions the predicted force is about 
half that observed. Whereas the theory predicts 
that this force will continually decrease as L in- 
creases, actually in the natural rubber com- 
pounds it passes through a minimum and then 
increases. In the case of the Hycar compound C 
this rise does not appear. The effect appears to 
be due to the increasing importance, as the chains 
are straightened and aligned, of those van der 
Waals forces which eventually lead to crystalliza- 
tion of the natural compounds. 

The values of « indicate immediately that in 
compounds of the type here considered the chains 
of the molecular network are, in the unstretched 
state of the bulk material, extended to some 10 
percent of their maximum extension. 

By combining this result with the observed 
Young’s modulus of the material one can gain 
some idea as to the density of the active coherent 
network. At small extensions we have 


E=9MkT/NP, (7.7) 


M being as before the effective number of chains 
in one coordinate direction and N the number of 
links per chain crossing a unit cube of unstretched 
material. Since each link contributes an amount / 
to the maximum extension of the chain we have 


N=1/l. (7.8) 


Now as a measure of the fraction of the material 
actively involved in the network we may take 
the quantity 


A=3MNm/Ap, (7.9) 
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A being Avogadro’s number, p the density of the 
material, and m the atomic weight per link in 
the chain. Combining these expressions, one 
obtains 


A=Em/3KApkT. (7.10) 


Each link in our model chain will correspond, 
not to a carbon-carbon bond, as these are not 
independent of each other, but to a segment 
of chain just long enough for the orientations of 
neighboring sections to be nearly independent. 
If we suppose that each isoprene group should 
be represented by the links we have m=35. 
Introducing experimental values into Eq. (7.10) 
we then find A~}. We thus estimate that some 25 
percent of the molecular chains are active in the 
formation of the network, the remainder forming 
a relatively inert medium through which the net- 
work extends. This result, which appears to us 
to be physically reasonable, of course applies only 
to compounds similar to those here considered, 
for it may be much affected by the degree of 
vulcanization of the material. Increasing vul- 
canization, forming new links within the network 
and tending to link new material into it will tend 
to increase x, E, M, and probably A. It would be 
very interesting to have more information on the 
way in which these quantities depend on the 
vulcanization of the material. 

The idea of a molecular weight does not appear 
in the network theory of rubber, which considers 
essentially a single enormous molecule. There is, 
however, an interesting quantity, the mean 
weight of the chain segments between junctions 
in the network, which plays essentially the same 
role in this theory as mean molecular weight does 
in the theory of Kuhn. To determine this quan- 
tity one must go beyond the proper bounds of the 
present work, setting up a more detailed theory 
of the structure of the network in unstressed 
rubber, presumably on the basis of consideration 
of the processes by which it is formed. 

Kuhn’s theory leads to an estimated molecular 
weight because it does involve a definite picture 
of the conditions in unstressed rubber, the as- 
sumption that the distribution of extensions of 
the molecules is then essentially that which would 
obtain at any instant if the molecules were freely 
suspended in a fluid. On this basis the mean con- 
tribution of a molecule to the elastic constant 
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becomes independent of its length, and one can 
determine from this constant the total number 
of molecules irrespective of length. With Kuhn’s 
effective distribution of extensions one has 


E(1) =7GokT, (7.11a) 


G» being the total number of molecules per unit 
volume. An at least equally reasonable theory 
which considers only the distribution of distances 
between the ends of the molecules would give 


E(1) =3GokT. (7.11) 


(Note that under these assumptions the valence 
angle does not enter into the expression for the 
molecular weight.) From Gop one can at once 
determine a mean molecular weight 


u=pA /Go. (7.12) 


Now we have shown in the Appendix that in 
the network theory of rubber (at least in the case 
of Gaussian chains) the mean positions of the 
junctions of the network move in just the way 
that Kuhn assumes the ends of the molecules to 
move when the material is distorted. Correspond- 
ingly the lengths of chain between these junctions 
will enter our theory in the same way that the 
molecular weight enters Kuhn’s. If we were to 
assume that the mean extensions of the chain 
segments were related to their lengths by the 
same distribution function which would apply to 
thermal agitation of the free segments, we would 
obtain the same result as in the above modifica- 
cation of Kuhn’s theory, with Go now the number 
of active chain segments between junctions. 
For the mean weight of chains between junctions 
we would then obtain 


u=pAA/Go, (7.13) 


the presence of the factor A marking an essential 
change from the point of view of Kuhn. The mean 
weight of a chain between junctions thus comes 
out on our theory to be of the order of one-fourth 
of Kuhn’s molecular weight. Combining this re- 
lation with (4.15) and (4.16), one has 


p=m/k’*, (7.14) 


while the number of links between junctions is 
simply 


u/m=1/*?. (7.15) 
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On this basis we would estimate that the number 
of links between junctions is of the order of 100 
in the materials considered. We do not wish to 
emphasize this result unduly, as it is not clear 
that it would arise from a more detailed theory 
of network structure, in which, we suspect, the 
relatively large extensions would appear more 
frequently than in the thermal distribution. 

One of us (H. M. J.) wishes to express his ap- 
preciation of the grant of a Guggenheim Me- 
morial Fellowship, during the tenure of which 
part of this work was done. We wish to acknowl- 
edge very helpful discussions with Professor R. 
H. Fowler during the early part of our work on 
this problem. We are also indebted to Professor 
G. E. Uhlenbeck for reading the manuscript of 
our paper and for constructive suggestions. 


APPENDIX. FORCES EXERTED BY AN IRREGU- 
LAR NETWORK OF FLEXIBLE CHAINS: 
REDUCTION OF AN IRREGULAR 
NETWORK TO A REGULAR 
NETWORK 


In this Appendix we shall discuss the average forces 
exerted by an irregular network of Gaussian flexible 
‘chains, by a method which involves rather detailed con- 
sideration of the statistical behavior of such a network. 
The method, as an approach to the final result, is not the 
most direct one possible, but it gives a particularly clear 
picture of the situation and provides useful intermediate 
results. It seems to offer a basis for further study of the 
behavior of these networks, and indeed of non-Gaussian 
networks, to which our basic theorem equally applies. 

Consider a network of flexible chains of arbitrary type, 
certain junction points of which are subject to external 
constraints which may, for instance, force the network to 
extend through a given volume in space. The constrained 
junctions (on which alone the network can exert a non-zero 
average force) will be indicated by roman letter subscripts, 
their positions being indicated by fa or xa, Ya, a, While 
free junctions will be distinguished by numerical sub- 
scripts; greek letters may indicate either type. 

As already noted, we can compute the average com- 
ponents of force exerted by the network on fixed junction 
a as 


Fra(x, y, 2) = ere log L(x, y, 2), (8.1) 


Fya(, 9,8) = -iTS log 2(x, 9, 2), 


where (x, y, z) is the number of configurations available 
to the network when the fixed junction a is at the point 
x,y, 2, all other constraints on the network being taken 
into account. 

One can also compute this average force as the sum of 
the average forces due to all chains reaching the ath 
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junction. The average force due to the chain running from 
the ath junction to the rth junction (fixed or free) can be 
computed by treating the rth junction as fixed and aver- 
aging the results over all possible positions of the rth 
junction, with weights proportional to the probabilities 
with which these positions occur. We shall show that this 
procedure is valid, in the absence of steric hindrances, by 
proving that it gives the same result as the preceding 
method. 

Let the number of links in the chain between the ath 
and rth junction be Na;, the number of configurations 
available to this chain be c(Na;, fa—fr). Then, neglecting 
steric hindrances, the number of configurations available 
to the network as a whole, when the position of every 
junction is specified, is 


C(f1, f2, °**, Tay To°**) 


=c(Nie, T1—f2)c(No3, f2—13)+--, (8.2) 


with one factor on the right for every chain joining a pair 
of junctions. The number of configurations available to 
the network for specified positions of the fixed junctions 
only is 


@(ra, to, -++)= fdV-Clt, ts, -++4a, t6°°*), (8.3) 


the indicated integration being over the positions of all free 
junctions. If we indicate by /dV, integration over the 
coordinates of all free junctions except the rth, then the 
probability distribution for the rth free junction is 


P;(f2; fa, To, °° +) 


= {av,C(n, fo, °° yf fa VC(rite, ---). (8.4) 


As indicated, this is a function of the position of all fixed 
junctions as well as of the rth free one. 

In this notation the x component of force on the ath 
fixed junction is by the standard method, 


re] 
Fizq(fa; Tp, ee log L(ra, Tp, a? (8.5) 
Xa 


We now compute this force as given by the proposed 
method. 

The x component of force exerted on the ath junction 
by the chain running to the rth junction, considered as 
fixed at fr, is 


fear™ ~-7T2- log c(Nar; Ta—Tr). (8.6) 
OXa 


By the proposed method the total x component of force 
is then 


Fea=2r,f dt,Pe(tr; fa, Tb, ***)fzars (8.7) 


the sum being over all chains ending at the ath junction. 
Using (8.4) and (8.6) this becomes 


0 
Fra= —ET f dVC(n, ma? i 2+ log ¢(Nar, fa—Fr) 


/ Save, fo,°+-). (8.8) 


Remembering the form of C as given in (8.2), we can 
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express this as 


o . 
Fon ATS faVC(ty 12, --ta:-) / faVC(ry 1 -+-), 
(8.9) 


which is identical with (8.5) as required. 

This method is particularly easy to apply when the 
chains are all of Gaussian type. It then becomes a trivial 
matter to average over the positions of the free junctions, 
all of which are subject to probability distributions of 
Gaussian form. For a Gaussian chain between the rth 
and vth junctions we have as the x component of force on 
the vth junction 


2kT 
fare = ay, — >). (8.10) 


If P;(r,) is of Gaussian form, 


f dtp%,P y(t) =Xor, (8.11) 
where Xp, is the x coordinate of the center of the distribu- 
tion. The average x component of force between the rth 
and rth junctions is then 


J far,de,P, (r+) Py(tv)feve (fr - Ty) 


ee , 
—EN,, See) = fory(Tor fad. (8.12) 


Thus the average force between the junctions is the same 
as if both were fixed at the centers of their respective 
Gaussian distributions, and is entirely independent of the 
extent of these distributions. Accordingly, for the special 
case of Gaussian chains we can get the mean force on any 
fixed point by determining (1) the elastic constants 
(2kT/I2N) of the chains attached to this point, and (2) the 
mean position of the other ends of these chains. 

It remains to prove the above statement that all junc- 
tions of an irregular network of Gaussian chains have 
probability distributions of Gaussian form. We have, when 
all chains are Gaussian, 


1f |ri—re!? 
C(ri, Te, +-)=exp { —f[ 


4 Bl |}: (8.13) 


or 


1 1 
C(r1, f2, ree ) = exp {= 22 Nr» 


X [xr — xv)? + (2 — v)? + (Gr -s,)*J}, (8.14) 


in which expression we can sum over all pairs of junctions 
(r, v) if we take N,,= © whenever there is no chain con- 
necting the rth and rth junctions. 

Let us now consider the most probable configuration of 
the junctions in the network—the positions of the free 
junctions which make C a maximum. The coordinates of 
these points we shall denote by %o;, Yor: Zor, and we shall 
let xoqg=%a, etc., in the case of the fixed junctions. Then 
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these quantities can be determined from the equations 


4 


1 
Ly, “or — Xo») =0 


Ea -Wor— Yor) =0 ry=1,2,3,++, (8.15) 





2 (or —Zoy) = 0 

TiYrp 4 
where the sum is over all junctions, free or fixed. These 
equations determine the coordinates Xo,, Yor, Zor Of the free 
junctions as linear expressions in the corresponding coor- 
dinates of the fixed junctions: 


Nor =CraXatCroxpt+ eileen 


Yor =CraVatCroYot***, 
Sor =CraatCrotpt*:. 


(8.16) 


The coefficients ¢ which appear in these equations are 
functions of the N’s only. 
If we introduce as new variables 


X¢=X_r—Xoz; 
(8.17) 


Y;=Jr—Yors 
2, =27—Zor, 
specifying the position of each free junction with respect 


to its position in the most probable configuration, the ex- 
pression for C becomes 


Ilys ! 
C(t1, f2, +++) =exp {-} ee om 
tT>¥ TY 


x (or _ Xo»)? + (Yor oad Yoo)*+ (Zor —s Zon)4}} 


- 1 > 1 
ae {-5 ee Ny 


x [(ke— Xo) (Fe Vo)*+ (Ze 25). (8.18) 


In reducing (8.14) to this form we have used (8.15) and 
the fact that for fixed junctions X= Y=Z=0. The first 
factor on the right is, by (8.16), a quadratic form in the 
coordinates of the fixed junctions only, while the second 
factor involves coordinates of both free and fixed junctions. 

If one integrates (8.18) over the coordinates of the free 
junctions only the second factor is affected. One obtains 
then 


L (ra, Tb, +-)=Kexp4 G22 5 
TV TY 


[lor —son)*+ (Yor—Yor)™+ or —te)*T}, (8.19) 


with K a function of the N’s only. Using this expression 
together with (8.1), one can proceed to a discussion of the 
forces exerted by the network. Instead, we pass on to the 
evaluation of P,(r,). This we can write as 

1 


1 
Py(ty) = 7a, Tb, ** fav, exp{ 22 N, 


X ((X,—X,)?+(¥,— Y,)?+(Z, -2,)}}, (8.20) 


P,(t,) = pitta ty:++)1(X,)1(Vp)1(Zy), (8.21) 
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with 
1(X,) =f dX1-+-dXpsdXp4a--- 


Iyyliy _yx): 
xexp { - RUZ yo (Xe—X) }. (8.22) 
Now one can write the sum in the exponent as 
1 
(23) (X1- Ls, Xa -)/-+ OKs, Xa, +++), 
vr 4iVir 


where Xj is singled out for special treatment, and L is a 
linear form, Q a quadratic form in the other variables. 
Integrating over Xj, one gets 


1(X,)=kf dXo-+-dX,adX pyr 


xexp { —J0(%, Xa =) (8.23) 


where k depends only on the N’s. Proceeding similarly to 
carry out the other integrations, one finally obtains 


11 
1Xp)=CK'Pexp{—57Xe}, (8.24) 
where N, is the number of links in a chain which would 
give a similar distribution; it depends on the number of 
links in all the various chains in the network. We have 
then, from (8.21) and (8.24), 


K’ 1 1 < 
Py (ty) = (ta, Tp, ***) exp {- 57 t I}, (8.25) 


a Gaussian distribution, as was to be proved. It will be 
noted that the extent of the distribution is independent 
of the positions of the fixed junctions, i.e., the way in 
which the network is stretched. The position of its center 
turns out to be the point r,, and is related to the position 
of the fixed junctions by 


Toy =CyalatCutet::- (8.16’) 


We now note a useful and suggestive classical analogue 
to the network of Gaussian flexible chains. In this model 
each molecular chain, e.g., that of N;, links connecting 
the rth and vth junctions, is to be replaced by a classical 
elastic thread with the same elastic constant, 


Z/L=2kT /l?Ney. 
The equilibrium positions of the junctions of this clas- 
sical network are just the mean positions of the corre- 


sponding junctions of the molecular network. This becomes 
evident at once when we write (8.14) as 


(8.26) 





C(ri, fo, ---)=exp {—V(r1, 2, ---)/RT}, (8.27) 
where 
kT 1 
V(r, Ta, ** +e, *** ae zz Ne» 


X L(%2—Xv)*+ (Ye — yy)? + (27—2y)*] (8.28) 


is the potential energy of the classical thread model. The 
condition that C shall be a maximum is also the condition 
that V shall be a minimum; (8.15) determines not only the 
mean positions of the free junctions of the molecular 
chains, but also the equilibrium positions of the junctions 
in its classical analogue. 
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It follows at once that the classical analogue and the 
molecular network, being similar in all relevant respects, 
will exert exactly the same mean forces on their correspond- 
ing fixed junctions. Thus if one wishes to treat the elastic 
properties of massive rubber using a model involving a 
network of Gaussian chains, one can go a step farther in 
the discussion of forces and mean positions of junctions, 
completing a cycle, by returning to a classical model in 
which each molecular chain is replaced by a corresponding 
elastic thread. 

Now it is evident from (8.16), or on examination of the 
classical analogue, that if the x coordinates of all fixed 
junctions are increased by a constant factor the same will 
be true of the mean x coordinates of all free junctions, 
while the mean y and z coordinates will be unchanged. 
Generalizing this, it follows easily that if the fixed junctions 
of a network are given the motion of points in a classical 
uniform elastic continuum under uniform stretch in one 
or more directions, with or without superimposed shears, 
then the mean positions of the free junctions will also 
move like points in that elastic continuum. For example, 
if the fixed junctions of the network are constrained to lie 
on the surface of a parallelepiped, which is then distorted 
by stretches and shears, the interior free junctions will 
undergo, as concerns their mean positions, displacements 
corresponding to those same stretches and shears. It should 
be noted, however, that there is no corresponding effect 
on the form or extent of the probability distributions for 
the positions of the free junctions. 

Another observation of importance is that the forces 
exerted by a network are the same whether any given 
junction is treated as free, or as fixed at what would be its 
most probable position if it were free. This is evident from 
(8.12), or from the fact that the classical analogue to the 
chain network is the same in either case. This observation, 
together with the preceding one, forms the basis of our 
discussion in Section 7 of the relation of our model to that 
of Kuhn. If the network is regarded as imbedded in and 
exerting its forces, at the fixed junctions, on a continuum 
which is subject only to uniform distortions, one obtains 
a model which corresponds to that of Kuhn. Since we can 
consider all junctions as fixed, as above, and acting on the 
continuum, it is evident that the connectivity of the 
network is not of importance; we can introduce a parallel 
displacement of the fixed ends of any chains we may select 
without changing the net force on the continuum, and in 
this way we could break up or modify the connectivity of 
the network to any desired degree, arriving in the limit at 
a model which is essentially that of Kuhn. What are of 
prime importance in determining the magnitude of the 
forces exerted by a network are the mean separation and 
relative orientation of the ends of each chain; the con- 
nectivities enter only through their influence on these 
quantities. Indeed, to determine the effective elastic con- 
stants of a Gaussian network one need only specify the 
number of chain segments and the distribution of their 
mean fractional extensions (since these determine the 
mean forces); in the isotropic material all orientations of 
the segments would of course be equally probable. 








ill 



























Now, as in our model for bulk rubber, let us consider 
a network of Gaussian chains stretched through a unit 
cube by having certain junctions constrained to lie on each 
of the six cube faces, x=0, x=1, y=0, y=1, z=0, c=1. 
The free junctions will then all have mean positions within 
or on the surface of the cube. If the separation between the 
x faces is changed from 1 to L, there will be a corresponding 
change in the difference in the x coordinates of any pair 
of junctions, the force between those junctions [cf. (8.10) ], 
and the total normal force on the x faces. The y and z 
coordinates of all junctions will, however, be unchanged, 
as will the total forces between the y planes and thez planes. 
In short, the normal force between opposing forces of the 
model will be proportional to the separation of those faces, 
and independent of all other dimensions. This is exactly 
the behavior of a model ‘in which sets of independent 
chains run between opposite faces of the cube without 
interfering with or affecting each other in the least. If the 
network is isotropic in the mean then this equivalent but 
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simplified network can be chosen as having equal numbers 
M, of chains, with equal numbers of links, N, stretching 
between opposite faces of the cube. This is the model with 
which, for ease of visualization, we have made our cal- 
culations. It should be noted that the forces between 
opposite faces depend only on M/N so long as we deal 
with Gaussian networks. Thus it is not possible to deter- 
mine either M or N separately from any observation on 
rubber under conditions in which this model is adequate. 

It is easily shown that this simplified model will give 
results equivalent to those with the original irregular 
network when shears are present. 

It can be shown by simple examples that reduction of 
an arbitrary irregular network to an equivalent set of 
three independent bundles of chains of equal length, 
acting in the three coordinate directions, is not always 
possible. So long as the departure from Gaussian behavior 
is not extreme it should, however, provide a serviceable 
approximation. 
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At temperatures not too close to its melting point, the molecules in a crystalline solid vibrate 


with simple harmonic motion about fixed equilibrium positions in a space lattice. In order to 
calculate the thermodynamic properties of the solid, one must know the distribution of fre- 
quencies of the normal modes of vibration of the component oscillators. On the basis of the 
Born-Karman model, the frequencies of the normal modes are the characteristic roots of a 
matrix. The moments of the frequency distribution can be found from the traces of powers of 
the matrix. Methods are developed for calculating the frequency spectrum and thermodynamic 
functions from these moments. The frequency spectrum of a simple cubic lattice is studied 
by deriving’ its polynomial approximation from its moments. The spectrum has two maxima, 
one near the middle of its frequency range and the other at the high frequency end. As the 
ratio of the interactions between next nearest neighbors to those of nearest neighbors increases, 
the height of the low frequency maximum increases and that of the high frequency maximum 


INTRODUCTION 


Fl peoeagetinemia NE solid, to a first approxima- 
tion, can be pictured as a set of atoms or 
molecules whose equilibrium positions form a 
periodic lattice in space and which at a given 
time, due to thermal motion or interactions with 
neighbors, may be displaced slightly from their 
equilibrium positions. Neglecting non-harmonic 
displacements (which tend to become increas- 
ingly important at higher temperatures) every 


decreases. Specific heats are calculated from the moments of the frequency distribution. 





possible state of motion of the crystal is a linear 
combination of its so-called normal modes of 
vibration. 

According to classical mechanics and classical 
statistics, to find the thermodynamic properties 
of such a system one needs only to know its 
number of normal modes or degrees of freedom.' 
For example, the specific heat C,, at constant 


‘J. Mayer and M. Mayer, Statistical Mechanics (J. Wiley 
and Sons, New York, 1940), Chapter 11. 
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volume is the product of the number of normal 
modes and Boltzmann’s constant k. Such a 
theory leads to the old Dulong-Petit result of a 
constant C, independent of temperature. To 
eliminate the contradiction between this classical 
result and the experimentally observed decrease 
in C, with temperature, Einstein? suggested that 
the contribution of a normal mode with a fre- 
quency »; is that of a quantized oscillator in the 
Planck sense, so 


Cy=k >; (hv;/2kT)?/sinh? (hv;/2kT), (1) 
where h is Planck’s constant and T is the absolute 
temperature. By assuming independence of mole- 
cules or atoms at individual lattice points and 
equality of frequencies of normal modes Einstein 
obtained his well-known results which agree 
qualitatively with experiments but quantita- 
tively deviate somewhat at low temperatures. 
Debye*—5 improved the situation at low tempera- 
tures by treating a solid as a continuous elastic 
medium in which the fraction of normal modes 
with frequencies in a small interval dv at v are 
proportional to v*dv. Now the Debye theory in- 
volves a parameter @p (related to the largest 
frequency), the Debye temperature of the solid, 
which theoretically should be a constant inde- 
pendent of the temperature. Actually if one plots 
as a function of the temperature the values of Op 
necessary to make the Debye C, agree exactly 
with the experimental ones at the same tempera- 
tures, the resulting curve is far from a constant 
in the low temperature region.® Born and Kar- 
man’ have proposed a molecular approach to the 
frequency spectrum problem that seems to ex- 
plain the variation of Op at low temperatures and 
which, at the same time, is sounder than the 
Debye model in that it provides a molecular 
picture of a molecular situation. 

Born and Karman in a straightforward manner 
deduce the equations of motion for the set of 


2A. Einstein, Ann. d. Physik 22, 180 (1906); 34, 170 
(1911). 

3P. Debye, Ann. d. Physik 39, 789 (1912). 

4E. Schroedinger, Handbuch der Physik, X (1926). 

5R. Fowler, Statistical Mechanics (University Press, 
Cambridge, and Macmillan Company, New York, 1936), 
Chapter IV. 

6M. Blackman, Proc. Roy. Soc. 148, 365, 384 (1935); 
159, 416 (1937). 

7™M. Born and T. von Karman, Physik. Zeits. 13, 297 
(1912); 14, 15 (1913). 
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atoms in the crystal and by substituting periodic 
solutions into these equations they find that the 
squares of the frequencies of the normal modes 
of vibration are characteristic values of a matrix 
which we shall leave undefined at the present and 
merely signify by M*. That is, to find the fre- 
quencies of the normal modes one must find the 
characteristic values of 


My=y = vy. (2) 


(y is an unknown vector which satisfies (2) and 
which shall not interest us here.) For every 

;= v2 which is a characteristic value of 1, there 
are two frequencies +»; and —»;, the positive 
one of which is to be used for finding thermo- 
dynamic properties (i.e., to be substituted into 
equations like (1)). The frequencies {»;} turn out 
to be so densely distributed between — vz and v, 
(A, = vz? is the largest characteristic value of (2)) 
that one can define a continuous function g(v) 
such that g(v)dv/3N is the fraction of frequencies 
in the interval dv at v. From (2) it is apparent 
that 


g(v) =g(—») 
causing the odd moments of the distribution 
function g(v) to vanish: 


vL 
f pntlo(v)dv=0 (3a) 
ae 


VL VL 
f vo(v)dy= 2f g(v)v?"dp. (3b) 
—L 0 


Now it is well known that if {A;} is the set of 
characteristic values of the matrix M? 


Trace M’@=>0;\;=)D; v? 
and 
Trace M*™=> AF =>; v2. 


If the order of M? is 3N, the average value of the 
2kth powers of the frequencies are 


(1/3) Trace M%=(1/3N) ¥; v2 
and since this by definition is the 2kth moment 
of g(v) 


vL 
g(v) v*dp 


—vL 


(1/6N) 


=(1/3N) Trace M*=px. (4) 
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With a knowledge of M? for a given crystal one 
can (at least with an infinite amount of patience) 
calculate all the moments of the function g(v) 
and the problem of determining g(v) is reduced 
to the moment problem of statistics—given a set 
of moments {yx}, to find the distribution func- 
tion which has these moments.® 

The first part of this paper will be devoted to 
outlining some applications of the moments of 
g(v) to the direct calculation of thermodynamic 
quantities and to a discussion of methods of solu- 
tion of the moment problem. In the second part 
a detailed study of the Born-K4rman model of a 
monatomic crystal with a simple cubic lattice 
will be made. 


I. USES OF THE MOMENTS OF THE 
FREQUENCY SPECTRUM 


A. Calculation of Thermodynamic Quantities 


1. General Considerations 


In order to calculate the average value of a 
thermodynamic quantity f (as a function of the 
temperature) of a crystalline solid one must 
evaluate integrals of the form 


f(T) = f fv, T)g(v)dv. (5) 


For certain models of crystals (with 3N degrees 
of freedom) it is not difficult to evaluate the 
moments? 


w= f v*g(v)dv/3N (6) 
0 


of g(v). With this in mind, we shall derive an 
expression for f(T) involving the moments of 
g(v) instead of g(v) itself. Let us define 


F(x, T) =2/9'f fo, T) cos vxdv, (7a) 


G(x) =(2/x)! f iii a 
0 


Then by the well-known Parseval formula® 


f(T) = f F(x, THG(e)de. (8) 


+ Montroll, J. Chem. Phys. 10, 218 (1942). 
Titchmarsh, Introduction to the Theory of Fourier 
literal (Oxford University Press, Oxford, 1937), p. 50. 
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The quantity F(x, 7) can usually be expressed 
in closed form—for example, consider the tem- 
perature-dependent component of the average 
internal energy: 


E(T)= f hvg(v)dv/(e”!*? —1). (9) 


The cosine transform of 
e(v, T) =hv/(—1+ exp (hv/kT)) 


is the integral 


2\! ” @y cos vxdv 
E(x, T) -(-) aT | -- + -—-, 
Tv 0 e”’—1 


when 6=h/kT. It is not difficult to show that!® 
*sinazdz e*+1 
J e—1 A(e*—1) 2a 
Therefore, differentiation with respect to a yields 
af = azdz 1 e 
and 0 e—1 a (e*—1)? 


Qm*h 2\4) 6 e2rz/0 
E(x, T)= (-)| ~ (12) 
T 


92 An 2x? (e2*2/9_1)2) 


(10) 





(11a) 


(11b) 





The cosine transformation of the distribution 
function g(v) is 


2\) p* 
aw)=(=) f a cos vxdy 
= (=n) 
-(- ) fz tdi 
0 (2n)! 
2\3 © (—1)"x2" p* 
={- —_———— v"9(v)dv 
(-) x (2n)! J ew) 


2 bw (-— 1)"x?" yu n 
-3y(-) "pin diet 
vy) = (2m)! 


Therefore, from (8), 


oe) 62 
E(T) = Ey+(12"hN/62) f | 
0 41?x? 





e27=/8 ~ (—1)"x*"yu ndX 
peIE CE 


(e2*2/6 — (2n) ! 


10 E, Whittaker and G. Watson, Modern Analysis (Cam- 
bridge University Press, Cambridge, 1935), p. 122. 


n=0 
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and 


2r fr” WX 
C,/3Nk=— f csch? — 
0/75 0 
wo (—1)"x2"yo,dx 
- » (15) 
n=0 (2n)! 





x TX 
xX 4— coth—-—1 
6 6 


Eo is proportional to the sum >»; which in a 
Born-Karman model is the trace of the positive 
square root of the matrix M?. 


2. High Temperature Series 


Of course at high temperatures a useful formula 
results immediately from the expansion of f(v, T) 
as a power series in 1/7. Again, consider the 
internal energy 


3N ‘ 
E(T) => {hvi/2+hv,/(e/*?—-1)}. (16) 
i=] 
Since 


x/(e#—1) =1—x/2—S. (—1)"Byx®*/(2n)!, 


n=l1 


where the B,’s are Bernoulli numbers and 


|x| <2z, 


an fT hy; hv; 
ar) =a E| ih cieee 
i=1 L2kT 2kT 





sn | 
n (2n)! 


when |hrmax/kT| <2. Now the moments of the 
frequency distribution are 


(1/3) f v'g(v)dv=>> v"/3N= un, 
0 


therefore! 


© (—1)"B,, 
E(T) =3NkT —3NkT > ———— 
(2n)! 


h 2n 
er 2n 17 
x(—) Han (17) 


_"' This type of development was first made by H. 
lhirring, Physik. Zeits. 14, 867 (1913); 15, 127, 180 (1914). 
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and 


«© (—1)"B,(1-—2 
C.=3Nk—3Nk & — : 
(2n)! 


n=1 
h 2n 
x(—) Ken, (18a) 
kT. 


T>hvmax/2nk. 





when 


Remembering that B;= 1/6; B,= 1/30; B;s= 1/42; 
- we have (if 6=h/kT) 
Pus Ou Aus 
C,/3Nk~1———-+—-- +--. 
12 240 6048 


(18b) 


B. Solution of the Moment Problem 


This section will be devoted to introducing 
several methods of finding the distribution func- 
tion g(v) from its moments. The first procedure 
will be based on the Fourier integral theorem and 
the others will be applications of the method of 
least squares. 

If the Fourier integral of g(v) is G(x) (Eq. (13)), 


G(x) = Q/nf £0) cos vxdv 


=3N(2/m)! Yo (—1)"x**pan/(2m)!, (13) 
n=0 
the inversion of (13) by the Fourier integral 
theorem yields 


g() = (7) f ‘Rowe 


- (14) 





6N r® © (—1)"x*"u2, COS vxdx 
J (2n)! 


T n=0 


In general, the power series representation of 
G(x) cannot be expressed in closed form in terms 
of well-known functions, so one must resort to 
numerical integration actually to obtain g(») 
from (14). This is unfortunate, for in order to 
obtain a complete g(v) curve one would have to 
make a numerical integration for every value of 
v; and, furthermore, the resulting g(v) would not 
be in an analytic form, as one might like to have 
it. However, as will be shown later in a discussion 





CRYSTALLINE SOLIDS 


of the method of least squares, some results 
derivable from (13) and (14) are useful in the 
determination of errors involved in the applica- 
tion of the method of least squares. 

Although the properties of Fourier integrals 
do not seem to be suited to our problem, those 
of Fourier series are. Since we are interested only 
in g(v) when |v| <vz (vz is the largest frequency 
of a normal mode) and 


g(r)=0 if |r| >vz, (15) 


we Can express g(v) in the interval <—v ,, vp> 
by the cosine series: 


g(r) = 09/243: a, cos (nrv/vz), 


n=1 


(16) 


where the coefficients {a,} are computed by the 
well-known scheme of multiplying (16) by 
cos (nmv/v,)dv and observing the orthogonality 
of the elements of the set {cos (nrv/v z)}: 


(17a) 


ao =po/V1 


and 


a,= (1/01) f g(v) cos (rnv/v,)dv 


“Te 


=(1/mn) f “0 


—VvL 


XE (—1)™(nav/v1)?"dv/(2m)! 


m=0 


=(6N/»2) z (- 1)"(na/v1)"uom/(2m)!. (17b) 


Thus the coefficients of the (cos mrv/vz)’s are 
linear combinations of the moments of g(y). 
This Fourier development of g(v) is really a 
special example of the method of least squares. 
Suppose we try to approximate g(v) by a function 
G(ao, a1, +++, x; vy) =G(a; v) of a set of para- 
meters {a,} chosen in such a manner that 


I(ao, a1, ++ +, de) =I (a) 


-{ [g(v)—G(a; v)]}*dv (18) 


(where a= (ay, ---, ay)) is a minimum. According 
to the calculus of variations, the necessary condi- 


tion for I(a) to be a minimum is 
aI(a)/da,=0; n=0, 1, -:: 


or 


vL 
f g(v)(0G/da,)dv 


YL 


VE 

-f G(a; v)(0G/da,)dv=0. (19b) 
=e 

The simplest way of finding a set of a’s to satisfy 


(18) is to choose G(a, v) as a linear combination 
of a set of functions 


go(x), ¢i(x), all $n(x) (x= v/v) 
orthonormal in the interval <—1, 1>, i.e., 


1 


f ¢i(x) ¢j(x)dx = 6 jj. 


1 


(20) 


For, in that case, 


G(a; v) => AnGn(v/vr) 


n=0 


and Eq. (19) becomes 


i) * edofonder 


L 


¥L 


= DX om(v/vr) on(v/v1)andv=0 


—vym 


or by (20) 


1 


a,= J g(xvr1) on(x)dx. 


—1 


(22) 


In the special case 
¢n(x) =Cos xn7, 


the coefficients are the linear combinations of the 
moments of {y2,} that were given in (17). And 
for any finite number k, the coefficients deter- 
mined in (17) are just those which give the best 
k-term Fourier approximation in the least square 
sense. This development has the disadvantage 
that each coefficient depends on all the moments 
of g(v). Actually it is possible to find other ap- 
proximations to g(v) such that a, depends only 
on the first 2 ye,’s.* For example, if 


on(x) =[(2n+1)/2 }P,(x), 
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where P,,(x) is the nth Legendre polynomial, 
1 
a,= C(an+1)/2} f g(xv1r)P,(x)dx. (23) 
—1 
When 1 is odd, a,=0, for g(xvz) Pax) is an even 
function and when 1 is even (say 2m) 
n (—1)"(4n—2r) !x?"-2 
Pon(x) _ Z ° 
r=0 22"71(2n —r)!(2n—2r)! 





Therefore 


axn=[(4n-+1)/2} f slees) 


n (—1)"(4n—2r) !x*"-*"dx 
r=0 22"7!(2n—r)!(2n—2r)! 
= 3N[(4n+1)/2]} 
n (—1)"(40—2r) !yon—2r 


r=0 7!(2n —r) !(2n — 2r) !y,2"— 24122" 





(24) 





and a2, depends only on the moments pa, with 
k<n. Finally 


ated => [(4n-+1)/2}arnPon(v/vz). (25) 


Another type of approximating function whose 
k parameters depend only on the first k even 





ao/1+a2/3+a4/5+-- 
ao/3+a2/5+as/7+-- 
Qo/S5+a2/7+a4/9+:- 


ay/(2k+1)+a2/(2k+3) +> >> 


Ww. 
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moments of g(v) is (assuming that x= v/vz) 


2k: 
++) x3 xX) =>, apXx". 


n=0 


G(ao, a1, 


Since g(v)=g(—v) and we would also want 
G(a : x) to be an even function, all the odd co- 
efficients d@2n41 vanish and 

k 


29+, Qxn3X) =D, donr. 
n=0 


(26) 


G(do, Qe, Wa, 


From (19b) 


f " g(o)(v/v1) dv 


—Ph 


k "L 
-> am f (v/v1)?™™dpv=0 (27) 


m=0 —vL 


n=0,1,2, ---,. 


VL 
i) o(»)»*"dy=6N/tan (28a) 


“Fh 


f “pomtmdy = 2v7/L2(m+n)+1]. (28b) 


Therefore the set of k+1 equations (27) become 


’ +ax./(2k+ 1) ai 3Nyo/v1, 
»+a%./(2k+3) =3.Nyu2/v1*, 
. +ax%/(2k+5) = 3Nus/ vz’, 


+a2%/(4k+1) =3Nuoxe/vr*t'. 


Or, if we write the matrix of this set of linear equations as 


1 1/3 


A(k) = 


1/(2k+1) 
and define 
a(k) = (do, 2, aa, 


u(k) =(3N/vz) (uo, oe/ vr, wa/ v4, 


1/(2k+3) 


1/5 1/(2k+1) 


1/(4k+1) 


ee *, Qo), 


wt wor/vt"), 


our set of linear equations becomes the matrix equation 


a(k)-A(k) =k). 


Now suppose A~!(k) is the inverse of A(k). Then 


a(k) =u(k)-A~*(k) 
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and each component of a(k) can be found from yu(k) and A~!(k). By direct multiplication it is possible 


to verify that 
A~*(1) =1, 
( (1-3)(1-3)/1 


(1!)222\ —(3-5)(1-3)/3 


‘( 9 pe 
4\-15 45/° 
(1-3-5)(1-3-5)/1 


—~2(3-5-7)(1-3-5)/3 
(1-3-5)(5-7-9)/5 


A-\(3) =—_—_—_ 
“) (21224 


—2(1-3-5)(3-5-7)/3 
4(3-5+7)(3-5-7)/5 
~2(5-7-9)(3-5+7)/7 


{i oe) 
(3-5)(3-5)/5 


1(1-3-5)(5-7-9)/5 
—2(3-5+7)(5-7-9)/7 
(5-7-9)(5-7-9)/9 


and, in general, the element in the 7th row and jth column of A~'(k), a;.(k), is* 


k—-1\ /k—- 
(—1)‘*#(22 —1) (22 +1) (274+-3)- - -(22-++2k —3) (27 —1)(27+1)- - -(27+2k—3) CG.) 





(k—1)!2%-*[2(¢+ j) —3] 
(—1)4(27-+2k —3)'(2j+2k —3) (k—1)! 





~ 24—4(24 — 2) 1(2j —2) (i-+R—2) (Gtk —2) (kJ —2) (ki —2) [206 +) -3] 


Therefore 
° * 
a;= > f2j0;:(R) 
j=0 


(2i+2k—1)!(k—1)! k 


(34) 


(—1)*7(20-+2k —3) luo; 





~ 24—4(25 — 2) \(¢-+& —2) (k —i—2)! > (27-2) (G-+R—2) Kk —j—2) [2G+7)—3] 


Anyone using the method of least squares soon 
realizes that the accuracy with which an approxi- 
mating function represents an initial function 
after a best set of constants have been chosen 
depends significantly on the form of the approxi- 
mating function. With a given set of moments 
one can determine the best (in the least square 
sense) approximating function to g(v) of all those 
discussed above by seeing which minimizes the 
integral 


1a)= fi “E02 f " eGo; sae 


—vL 


VL 
+f G*(a; v)dv>0. 


—vL 


* The author is indebted to Mr. W. Bower of Cornell 
University for calculating the inverse matrices A~(k) for 
several small values of k. His calculations indicated how 
the general inverse could be obtained. 





(35) 





The component of I(a): /G*(a; v)dv depends 
only on the approximating function under in- 
vestigation, so its calculation causes no trouble. 
And when G(a; v) is taken to be a Fourier series, 
a sum of Legendre functions, or a polynomial, it 
is clear that 


i) " g()G(a; »)d> 


—vL 


can be evaluated in terms of the set of constants 
a and the moments of g(v). Inasmuch as we do 
not know g(v) itself, there is no direct way of 
finding {g*(v)dv. However, if G(x) is the Fourier 
integral of g(v), Parseval’s theorem yields 


[eed f ecodr= f cree. 


—o 


Now from direct calculation in special cases it 
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seems that G(x) has a strong greatest maximum 
at the origin. This indicates that G?(x) would 
have an even stronger greatest maximum there 
and that the integral of G?(x) could be calculated 
numerically. With a knowledge of the three com- 
ponents of I(a) one could find the mean square 
deviation of an approximating function G(a, v) 
from g(»). 

It is also possible to find the maximum error 
in thermodynamic quantities calculated on the 
basis of a frequency spectrum G(a; v) instead of 
the exact one g(v). Suppose we wish to find the 
thermodynamic function 


f(T) = f f(, T)g(»)dv. 


Writing 
g(v) =G(a; v)+€(r), 


where ¢(v) is the error in approximating g(v) by 
G(a; v) at », 


f(T) -{ f(y, T)G(a; drt J f(v, T)e(v)dv. 


Our problem now is to find an upper bound in 
the error resulting from writing 


f(T) = f fl, T)G(a: »)dy. 


Since 
a VL 
f f(y, Te(v)dv= J f(v, T)Le(v) — Gla; v) dv, 
0 0 
the Schwartz inequality” indicates that 


f f(v, T)Le(v) —G(a; v) dv 


< | f i ride] fee) —Ge nae]} 
-|| J “f(y, ry |1(@/2] 


%R. Courant and D. Hilbert, Methoden der Mathe- 
matischen Physik (J. Springer, Berlin, 1931), Vol. 1, p. 40. 
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and 


4(T)= | f(, T)G(a; »)a» 


vL } 
+4f (r0a)/2) f f(y, rar), (36) 


where ¢ is some number such that |¢| <1. 


II. BORN-KARMAN MODEL WITH 
CENTRAL FORCES 


This section will be devoted to developing a 
method of calculating the moments of the fre- 
quency distribution of a Born-Karmén crystal 
in which only central forces exist and in which 
only nearest and next nearest neighbors of a 
given atom directly influence its motion. From the 
moments corresponding to a three-dimensional 
simple cubic lattice, its frequency spectrum will 
be obtained. 

In following the Born-Karman approach one 
first writes the equations of motion for a system 
of N atoms. Assuming central forces, the total 
potential energy is 


N 
Viti, +++, tv) = DL Vixs(ti—r1)). 


i<j=1 


(37) 


Expanding this function in a Taylor series about 
the coordinates of the molecules at equilibrium 
at absolute zero, choosing the zero point of 
energy as the constant term obtained from such 
a development, remembering that, at equilib- 
rium, no forces act on any atom, and assuming 
that displacements from equilibrium are suffi- 
ciently small so that their powers higher than the 
second may be neglected, 


N 
Vitis tw) = DL (&i— Es)? Vij!" (ais)/2, (38) 


i<j=1 


where £; is the displacement vector of the ith 
atom from its absolute zero equilibrium position 
and V;;''(a;;) is the second derivative of the po- 
tential energy function due to interactions be- 
tween the ith and jth atoms evaluated at their 
equilibrium positions. Now if the interatomic 
forces are sufficiently short ranged so that only 
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nearest and next nearest neighbor interactions 
are important, V(rj, ---, ry) becomes a function 
of two parameters, a= V"’(a;;)/2 with i and ja 
pair of nearest neighbors and y= V’’(a,:)/2 with 
k and / a pair of next nearest neighbors. In this 
case 


Vity-sty)=a Do (Ei §,)? 


(nearest) 


+7¥ 2 


(next nearest) 


(E&—E:)*. (39) 


Of course a>y. 


M* =(2aA+4yB)" 
= (2a)"{A"-+(2y/a) © Am "BA 


ni=1 


+(2y/a)? ) > Ar—BAni—n2B A n2-2 


ni=2 no=2 


HQT oy 


ni=3 n2=3 n3=3 


KABA BA2-3BA ng—3 4 oa } 
(41) 


= (2a)" Do (2y/a)*"gs,n- 

Saving details until special cases are introduced — 
in the next section, let us merely recall® that the 
matrix of the secular equation from which one 
obtains the characteristic frequencies of a Born- 
Karman crystal is of the form 


M=2aA +4yB, 


If the order of M is m, the 2nth moment of the 
distribution of its characteristic values is 


Hon = Trace M"/m. (42) 


Since 2y/a is small (of order 0.1), we have, in 
(41) and (42), an expression for we, as a poly- 
nomial in 2y/a, whose first terms are the most 
important. 


(40) 


where A is a diagonal matrix and B a symmet- 
rical one. 





A. Moment Calculations for a Simple Cubic Lattice 


The frequencies of normal modes of vibration of a Born-Karman model of a simple cubic lattice 
of 2N X2N X2N identical atoms are solutions of the secular determinant® 


A11(¢1, ¢2, ¢3) —44*My? 412(¢1, 2) 413(¢1, ¢3) 


22(¢1, $2; $3) —4r* My? A23(¢2, $3) 


| 
| 2i(g2, ¢1) 
| 


a31(¢3, $1) A32(¢3, ¢2) 33(¢1, 2, ¢3) —40? My? 
with g;=2a;/N where the a;,’s are integers satisfying 
—N<d, a2, a3<N 
and 
a11=2a(1—cos g;) +4y(2—cos ¢) cos g2—COS ¢1 COS ¢3), 
ou 2a(1—cos ge) +4y(2—cos g2 cos gi —COS 2 COS ¢3), 
33 = 2a(1—cos ¢3) +4y(2—cos 93 Cos g1—COS 3 COS ¢2), 
4i;(¢i, oi) =4y sin g, sin gy 44)j. 
Incorporating the 4x°M into v*, the kth power of the entire Born-Karm4n matrix is 
(M*(—N, —N, —N) 0 0 --- O 0 
0 M*(—N, —N, -N+1) 0 --- 0 0 


0 0 --» 0 M*(N, N,N) 








(2N+1)* rows and columns J 
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TABLE I. Coefficients Cam of the 2nth moments pon/(2a)"= D Cnm(2y/a)™. 
m=0 


, m 
* 0 2 3 4 
n 











1 * * 
1 * * 


0 
1 
2 3/2 5 ¥ 

3 5/2 15 57/4 ia 

+ 35/8 20 179/4 57 703/16 7 

5 63/8 175/4 990/8 210 3515/16 4549/32 . 

6 231/16 189/2 10,365 /32 690 15,321/16 13,467/16 30,353 /64 2 
7 429/16 1617/8 26,075 /32 133,959/64 58,037/16 272,167 /64 212,471/64 827,031/512 








Since all of the small 14?(a;, a2, a3) matrices are of the same form 
Git Gi2 13 
M?*(a, a2, G3)= | G21 G22 23 
(G31 32 33 
. 


N J N 
Trace M*= } YO > Trace M**(a,, do, a3). 


aiy=—N ag=—N a3=—N 
So, we need only find the traces of one general triple (a1, d2, a3), and sum over all a;’ 


3 
Trace M?(a,aea3) = >. a;i(a1a203) 
i=1 
and 
Trace M?=3 > ay(aiaea3) =3 >> ai. 


a1a2a3 a1a2a3 


By using the fact that the trace of the mth power of the matrix (43) is 


3 3 3 
pw nel Zz QajaoHaga3° * *Aanai, 
aj=1 a2= 


1 an=1 


it is not ‘difficult to show (taking into account all symmetries of matrix elements and summation 
operations) that if we let 
r= LU2X, 


ai a2 a3 


Trace M°=3 = 1, 
Trace M*=3 2 au, 
Trace M4=3 & (a,;?+2a,2”), 
Trace M&=3 & (a4; +6041012? + 2412013023), 
Trace M8=3 & (a114+80117d12? + 4012723? +40110.220 12" + 8011012023031 + 20124), 
Trace M¥=3 = (4135+ 10011912? + 10011023712? + 1001170120 13032 + 100117012722 (46) 
4100130190290 9091-+ 10011012 10012%013052-+ 10011012015"), 
Trace M¥?=3 © (4119+ 120114012? + 2412+ 6011720712? + 120110220124 + 1801170124 
+ 12011920021? + 240417 12013033023 + 1801170127031? + 1201112023031 
+4110 130330320 22021 +4801 10125023031 + 24011012713" 33 + 1201 1012023°d 31 


+ 12021413? + 801270137032" + 6011012723"d33-+ 120117127d23") . 
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Traces of higher powers of M become longer and more tedious to calculate by this method, so it is 
desirable to attempt to develop power series expressions in (27/a) as was suggested in the last section. 
Actually to calculate the moments of our distribution function we proceed as follows: 


pwo= Trace M°/3(2N+1)*=1, 
u2= Trace M?/3(2N+1)*=(2N+1)-3 = an 


=(2N+1)-* & [2a(1—cos 91) +4y(2 —cos ¢1 Cos g2—COS ¥1 COS ¢s) |. 


However, from the expressions 


N 0 if k is odd 
> cos‘ ra/N=} 
a=—N 1+2!-*N(k!)/(R/2)!(k/2)! if k is even. 


N 0 if k is odd 
> sin’ ra/N= ; 
a=—N N(k)!/2*-"(k/2)'(k/2)! if k is even. 





It is clear that 


Ya1=(2N+1) and Yiaicos gi=1 


ws = 2(a+4y)(2N-+1)-*[(2N-+1)?+0(N2) ] 
=2of1+2(2y/a)]+O(1/N). 
Likewise, if c;=cos (wa;/N) and s;=sin (wa;/N), we have 
wa=(2N+1)-* B [(2a)2(1 —c1)?-+2(2a) (4y) (1 —e1) (2— crea — e108) 
+ (4y)*(2—crcs—cx¢s)?-+2(4y) 2512827] 
= (20)*[3/2+4(2y/a)-+5(2y/a)?J+0(1/N). | 


In a similar manner we have derived expressions for ye, us, «++, “is Which are correct to within terms 
of order (1/N). These results are collected in Table I where the coefficients Cnm of 


un/(2a)"= X Cum(2y/a)™ (47) 


are tabulated. 
From Eq. (42), 
Hon = (1/24N*) Trace M"+0(1/N) 


and in the Born-Karman model of a simple cubic lattice, the trace of the first coefficient in (41) 
becomes 


Trace qg1,,= Trace A” 


N 
=3 © (1—c:)"=12N?2 O (1-c1)", 


a\=—N 
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where c;=cos (ai;7/N). But, as can be verified for n>0," 


¥ (1—c1)"= NE (2n—1)!/n!(n—1) !2"-2+0(1/N) ] 


a1 
(1/24.N%) Trace qi,n=(2n—1)!/m!(n—1)!2"—-! 
and, to a first approximation, 
Hon (2a)"[(2m —1)!/n!(m —1) !2"-1+ O(2y/a) ]. 


To obtain the next approximation to ue, we must find 


Trace go,=Trace >> A*™BA"-!, 


nj=1 


If the A matrix is represented by 
Au 0 0 
A=|0 Az O 
0 0 Ass 


and the B matrix by B=(B;;), it is a simple matter to obtain 


ya 
n—1 n—ni ,ni—l n—ni ni—l 
Au Bu Ai Ao Bis Au By3A 33 | 


n—ni ni—l n—ni ni—l 





n—ni ni—l n—ni ni—l 


(Ass BuAu A33 B32A 22 po? ws 


A™" BA" l= 1Aso BuAn i te Ao. B23A33 | 


and 


Trace A""\BA"1=A}, ButAm Bn+Ais Bas 
is independent of ,. Therefore: 
Trace g2,=3n >, An""By 
=3n >> (1—c1)""1(2 —c1€2 — €1€3) 
= (6n)(4N?) Dra (1—c1)"" 


or, for n>1, 
[1/3(2N)*?] Trace go, n=n(2n—3)!/(n—1)!(m—2)!2"-%. 
So 








(2n—1)! n(2n—3)! 27\? 
pa»~(2a)"| +(2y/a) +0( (=) )I (50) 
n!(n—1)!2"-! (n—1)!(n—2)!2"-8 a 


The calculation of Trace q;,, and higher coefficients is more complicated, but can be accomplished in 
essentially the same manner. The results of these calculations are not as simple as those for Trace 
gin and Trace qg2,, so they are not recorded in their most general form here, but are incorporated 
into Table I. 


138 Note that when n=0, we must write (2n—1)!/(n—1)!=(2n)!n/(n!)2n=0!/2(0!) =}. 
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B. Frequency Spectrum for Born-Karman Model 
1. Frequency Spectrum at Low Frequencies 


If one expands all the trigonometric functions 
in the Born-Karman matrix into power series and 
neglects terms of order (a;/N)* and of order 
(2y/a)?, it is clear that when |a;|<N, the normal 
modes » are obtainable from the three equations 


4n?v12M =al_gi27+ (2y/a) (2912+ o2?+ 3?) J 
| +O[(2y/a)*; gi*], 

4m?yo?M = al g2?+ (2y/a)(2¢2?+ 912+ ¢3") ] 
+O[(2y/a)*; gi*], 

4m?ys?M =al gs? +(2y/a) (293? + ¢1?+ ¢2?) ] 


+O[(2y/a)?; gi*], 
Gi =a,7/ N. 


Therefore the number of frequencies less than » 
(if vy is very small) is the total number of lattice 


points of (a1, d2, a3) space which are included in 
the ellipsoid 


v?=a12(a+4y)/ MN? 
+-a22(2y)/MN?2+a3%(2y)/MN?, 


and the corresponding ellipsoids obtained by per- 
muting a; with a2 and a; with a3. This number is 
approximately the sum of the volumes of the 
three ellipsoids, and since, upon rotation of co- 
ordinate axes, all three ellipsoids can be shown 
to be identical, the number of frequencies less 
than » is 


n(v) = (3)(40/3)(N*°Mr*) /2y(a+4y)}. 


Now g(v)dv was defined to, be the number of 
frequencies in the frequency interval dy at »v; 
therefore 


g(v) =dn/dv=12eN*M1y2/2y(a+4y)). . 


They obtain the same g(v) that we have used in 
the last section (that is, the distribution function 
whose moments are given by Eq. (47)), and that 
we shall use in the next section. We incorporate 
47°M into v’, so 


g(v) =3.N* My? /2ry(a+4y)}. 


However, as has been shown by Blackman, the 
square of the largest frequency vz, is given by 
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Fic. 1. Frequency spectrum of simple cubic lattice. 
vig(v)/3(2N)3 plotted as a function of (v/vz) for r=0.075, 
S. 


oi, . 


4n’? Mvz?>=4(a+4y). Using this fact, incorpor- 
ating the 42°M into »’, and defining a new 
parameter 


t= 8y/v_,", (51) 
we finally obtain 


vig(v)/3(2N)* = (v/vz)?/2m*r, (52) 


when » is very small. 


2. Polynomial Approximation of Frequency 
Spectrum 


Using Eq. (35) and Table I one can calculate 
the coefficients a, of the polynomial approxi- 
mation 


7 
G(v/v1) =X dan(v/v1)™ 
n=0 


of g(v). In order to find the a,’s it is necessary to 
know ue2,/vz2" for various n’s. Table I gives this 
ratio as a function of two parameters, but by 
making the substitution (51), r=8y/vz?, it is 
possible to express pen/vz?" as a one-parameter 
function. With this transformation the coeffi- 
cients Dam of 


n 
J 9 
Hen/ VL?" = 2) Damt™ 


m=0 





ELLIOTT W. 


MONTROLL 
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TABLE II. Coefficients bam of the 2nth moments pon/vz2"= LD Bamt™. 
m=0 








1 2 3 


ee 





+ * * 
0 * * 
—1/2 3/4 * 
—3/4 9/8 1/32 
—15/16 119/64 —11/8 
—35/32 355/128 —225/64 
—315/256 7845/2048 —1725/256 
—693/512 20,531/4098 —91,861/8192 


NAUALAONHK SO 


429/2048 


303/256 
1515/512 
7941/1024 
34,419/2048 


25/1024 
—4389/1024 


10,313/4096 
~126,217/8192 


79,947/8192 —58,233/65,536 








TABLE III. Coefficients ¢m; of (53). The large number of significant figures is necessary because of the alternate nature 
of the signs. When sums are calculated in (53) the first four or five significant figures usually cancel out. 








% 


0 1 2 


Coefficients of 7 


5 6 7 





5079 —22.337 297.733 
15.7644 458.440 —11,104.119 

— 295.8016 5,394.480 —15,777.088 
2091.4481 —85,986.577 1,017,799.282 
—6937.2931 374,647.254 —5,196,670.653 
11,680.7148 —722,493.189 10,672,356.522 
—9661.1765 645,337.459 —9,829,994.779 
3115.7088 —217,697.878 3,367,593.382 


—1,751.79 
76,955.12 
—45,293.70 
—5,767,333.61 
31,787,179.83 
—67,064,585.68 
62,531,485.03 
—21,542,684.29 


Coefficient of 2% 
NOurh Wd OO 





— 13,216 
1,387,740 

— 23,591,576 
149,413,313 

— 448,239,940 
687,301,241 
—520,682,759 
154,488,291 


4,999.17 
—221,267.41 

— 80,087.89 
19,496,814.16 
—104,419,659.56 
218,258,344.87 

— 202,398,644.57 
69,437,373.44 


3,556.2 —9,563 
65,651.5 1,257,181 
—6,819,465.1 —26,662.539 
71,690,894.7 172,648,360 
—284,259,563.8 —557,900,611 
523,471,667.7 906,057,942 
—452,304,050.3 —718,614,110 
148,283,675.4 221,355,300 








TABLE IV. Coefficients p, of vig(v/v1)/3(2N)§=Zp,7". 








b3 X10? 


bs X1073 


bs X10-4 De X10 pz X1076 bs X107 





1.8608 
— 1.1879 
— 3.7038 
— 3.3595 
— 1.5638 

1.1881 

3.6712 

4.3579 

2.3076 
— 1.7894 
—5.1835 
— 4.2461 

2.1389 

8.2282 

1.9310 


1.4819 — 14.8749 


—0.9922 
9594 

— 2.3246 
1.9042 
5696 

— 1.2908 
— 2.8309 
— 3.0462 
— 1.3744 
1.6279 
4.0281 
3.3304 
— 1.2523 
— 5.8836 
— 2.3062 
9.1985 


0.0007 
0110 
— .0020 
— .0113 
— .0147 
— .0094 
0045 
0171 
.0190 
.0054 
— .0163 
— .0268 
— .0084 
0277 
0262 
— .0435 


—0.0001 
0012 
.0004 

— .0006 
— 0013 
— .0012 
— .0004 
0008 
0015 
0010 
— .0004 
— .0016 
— .0011 
.0010 
.0017 
— .0019 


0.2797 
—.2978 
—.7017 
— 5616 
—.1351 
4554 
9392 
.9924 
4319 
—.5651 
—1.3715 
—1.1671 
3315 
1.9172 
.9061 

— 2.8059 


— 0.0360 
0082 
.0930 
.0999 
0546 

— .0365 
— .1336 
— .1736 
— .1060 
0595 
2225 
.2255 
— .0165 
— .3187 
— .1935 
4689 








are tabulated in Table II. Introduction of these 
moments into (33) permit one to calculate a 
seven-term polynomial approximation to g(v) 


(x=v/vz): 


vrg(v)/3(2N)?= X (2d Cust™)4™. (53) 


The Cms's are tabulated in Table III. To facilitate 
numerical calculations, >> s¢msx?* has been calcu- 
lated for various x’s in the interval <0,1> and 
tabulated in Table IV. This table is sufficient to 
give vzg(v)/3(2N)* for all expected values of r. 
The largest relative error probably occurs with 
very small x’s, so in this range one should use the 


asymptotic expression (52). In Fig. 1 we have 
plotted vzg(v)/3(2N)* for r=0.075, 0.1, and 0.15, 
on the basis of Table IV and (52). It is to be 
noted that there are two maxima in each curve, 
with the height of the low frequency maximum 
increasing with increasing 7 and the height of the 
high frequency maximum decreasing with increas- 
ing r. If one interpolates between the 7=0.075 
and 7=0.10 excellent agreement is obtained with 
Blackman’s® calculations for 7=0.085. 


C. Specific Heats 


For a simple cubic lattice it is possible to 
calculate the specific heat at high temperatures 
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from Eq. (18). This series converges when Thus, at low temperatures 
kT /hv,>0.159; the first fourteen moments tabu- 7 
lated in Table II being sufficient to determine Co/3hk(2N)'—~(kT/hv1)?/2n? 
C, to within one percent when kT/hv,>0.2. At 
higher temperatures 


( 3b] — 2 - 
Cy/3(2N)®k™1 — (RT /hv1)?/24 and at very low temperatures the upper limit can 
+(kT/hvz)*(3 —47+67?)/1920—. (54) be made infinite. Since 


hv L/kT 
xf x'e7dx/(1—e*)?, 
0 


At very low temperatures we can derive an 


asymptotic T* law by introducing (52) into J xte*dx/(1 —e*)? = 442/15, 


we have 


Cok f g(v)(hv/2kT)*dv/sinh? (hv/2kT). C./3k(2N)*~(2n2/152)(kT/hyz)*. (55) 
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Statistical Mechanics of Binary Mixtures 


T. ALFREY AND H. MARK 
Polytechnic Institute of Brooklyn, Brooklyn, New York 
August 11, 1943 


N a previous brief note! on the statistical mechanics of 
binary mixtures, we started with the usual assumptions 
—that the free energy of mixing is determined by a purely 
configurational partition function, that a given configura- 
tion can be described by specifying the occupants of a set 
of lattice points, and that the energy of interaction is pre- 
dominantly that of nearest neighbors. We then introduced 
a method of describing a configuration in terms of the 
occupants of a reference set of non-overlapping lattice 
point pairs, rather than of the lattice points themselves. 
This was done in order to avoid the usual complexity of the 
relation between the variables used in determining the 
distinctness of a configuration and those used in calculating 
the energy of a configuration. 

In the process of counting configurations, each lattice 
point is ordinarily considered separately; i.e., the lattice 
specifications are the elements of the statistical analysis. 
However, since the configurational energy depends upon 
interaction, it cannot be expressed as a simple function of 
the separate lattice point specifications. This energy is a 
sum of terms, each of which depends upon the simultaneous 
specification of two different lattice points. 

By the use of the above lattice point pair method of 
specification, and the further assumption that the probable 
total number of molecular pairs of a given type was directly 
proportional to the number of such pairs in the reference 
set of lattice point pairs, it was possible to derive directly 
and simply the equilibrium configurational relationships 
for a binary mixture. 

In a recent paper? Dr. Kottler has presented objections 
to our treatment. Dr. Kottler states that we did not 
present a sufficient justification for replacing the configura- 
tional partition function by its maximum term. As Dr. 
Kottler himself points out, this is a step ‘‘which is usually 
taken for granted.” We felt this mathematical step to be 
a well-known and standard procedure in the evaluation of 
this type of sum, and to require no elaborate justification 
in a communication as brief as ours. 

Dr. Kottler further objects to the assumption that the 
probable total number of pairs of a given type is propor- 
tional to the number of such pairs in the reference set. 
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This objection is based upon the argument that “the pairs 
overlap and therefore are interdependent.” We specifically 
stated as a restriction on the reference set of lattice point 
pairs that they were to be chosen in such a fashion that 
each lattice point was the member of one and only one 
reference pair and in this way precluded overlapping in the 
reference set of pairs. It is true that Maz, the total 
number of AB pairs among the total pairs of all kinds, is 
not geometrically identical with g-X, where X is the number 
of AB pairs among the reference pairs. We reaffirm, how- 
ever, that the probable value of Maz is given by 2-X 
and that in the absence of any long range order the dis- 
tribution function for a given value of X is very sharp 
about 2-X. 


1T. Alfrey and H. Mark, J. Chem. Phys. 10, 303 (1942). 
2F. Kottler, J. Chem. Phys. 11, 153 (1943). 





Specific Heat Ratios in CO, 


LEON Katz 
Electronic Department, Westinghouse Research Laboratories, 
East Pittsburgh, Pennsylvania 
August 9, 1943 


ABLE I showing experimentally measured values of 
the ratio of specific heats for a number of gases 
(y° exp.) was recently published by Clark and Katz. 


TABLE I. 








Temp os 
Gas Deg. Cent. 7 exp. * % Diff. 


Argon 24.2 1.6667 F 0.00 
Helium J 1.6669 d —0.01 
Nitrogen F 1.4006 of +0.05 
Hydrogen i 1.404; 
Sulphur dioxide S. 1.2642 
Nitrous oxide 3S. 1.2783 
Carbon dioxide 3 1.294¢ 











We measured gamma as a function of pressure at con- 
stant temperature in the range of 1 to 20 atmos. and then 
extrapolated by the method of least squares to zero pres- 
sure. These experimentally determined zero pressure values 
(7° exp.) were then compared to y° spec., calculated from 
the spectroscopically determined vibrational energy levels. 
It will be seen that excellent agreement was obtained for 
monatomic and diatomic gases, but rather large dis- 
crepancies were found for the triatomic gases N,O and 
CO but not SO». These discrepancies are in the opinion 
of the authors outside the limits of the accuracy of their 
measurements. 

It is the purpose of this letter to point out that similar 
discrepancies exist between all reliable experimentally 
determined values of y for CO found in the literature and 
those calculated by the spectroscopic method. 

Since most published values are given at one atmos- 
phere, it is necessary first to reduce these to zero pressure. 
This was done with the aid of Berthelot’s equation of 
state, which gives 


(1+)? 
1+ (7°—1)0(1—8) © 





(y'—1) = (7°—1) 
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) py BeetHecots Equation 


XExPemmmentar VALUE 
BY Crarmn 4nd KATZ 
VALVE CALCULATED From 
P+V-T DATA OF Mic KEL 


mies. 


TemPerRATuRE IN DeGrees CENTIGRADE 
10 zo 3o 40 


Fic. 1. 


where '!=ratio of specific heats at 1 atmosphere, 
°=ratio of specific heats at zero pressure, 
_27 1 4 (2). 
“32 PANT 
with P,.=critical pressure=72.85 atmos. for COs, 
T.=critical temperature = 304,19°C for COz. 


Using this equation and y° spec., the values of y! spec. 
are calculated, giving (y! spec.—~y® spec.) Berth. 

This difference, for COs, is shown in Fig. 1 as a function 
of temperature. A value experimentally determined by 
Clark and Katz and one calculated from P— V—T data of 
Michel? are also shown. The difference between these two 


TABLE II. 











Source y' exp. 





Partington! 
Partington and Cant? 
Shilling® 

Euchen and Lude+ 
Euchen and Ludet 
King and Partington 
Hubbard and Hodge* 
Sherret and Griffiths? 
Jatkar® 

Kistiakowski and Rice® 
Clark and Katz! 


1 3028 
1.3000 








1J. R. Partington, Proc. Roy. Soc. London A100, 27-49 (1921). 
asda Partington and H. J. Cant, Proc. Roy. Soc. London 369-80 

3W. G. Shilling, Phil. Mag. 3, 273-301 (1927). 

‘A. Eucken and Lude, Zeits. f. physik. Chemie 5B, 413 (1929). 

5 J. R. Partington and F. E. King, Phil. Mag. 9, 1020-26 (1930). 

3s C. Hubbard and A. W. Hodge, J. Chem. Phys. 5, 978-79 (1937). 

7 ore and Griffiths, Proc. Roy. Soc. London A156, 504-17 (1936). 

K. K. Jatkar, Ind. Inst. + 1-11 (1937). 
ra B. Kistiakowski and W. W Rice, J. Chem. Phys. 7, 281 (1939). 
0 See reference 1, in text. 
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Ratio oF Specific Heats ¥° 


—_ TemPERATuRE IN Decrees CENTIGRADE 
~— 10 20 30 40 


Fic. 2. x Partington; © Partington and Cant; ollie Shilling; 


| 
; Eucken and Lude; A\ King and Partington; @ Hubbard and Hodge; 


Kistiakowski and Rice; 


| 
(©) Sherret and Griffiths; ? Jatkar; F 


@® Clark and Katz; 


Straight line through experimental values by 
method of least squares; 


Spectroscopic values. 
points and the above ‘“‘Berthelot” values corresponds to a 
difference of about 0.07 percent in y. 

Table II summarizes the most reliable values of gamma 
at one atmosphere (y' exp.) for CO» as found in the 
literature and their corresponding values at zero pressure 
(y° exp.) obtained from the one atmosphere values with 
the aid of Fig. 1. Where we have assumed that 


(y' exp.—~° exp.) = (7! spec. —y° spec.) 


and that this difference can be calculated with sufficient 
accuracy from Berthelot’s equation of state. 

These values of +° exp. are compared with the spectro- 
scopic values in Fig. 2. It will be seen that all the experi- 
mental values, except the one by Kistiakowski and Rice 
are higher in value than the corresponding calculated ones. 

It would be difficult to see why the values of Kistia- 
kowski and Rice should be considered more accurate than 
those of the other workers. In their paper Kistiakowski 
and Rice point out in a footnote, that some of their earlier 
measurements of C, gave values 0.4 percent lower than 
those published. This difference is in the proper direction 
to throw y° exp. among the other experimental values. 

The straight line shown dotted was fitted to the experi- 
mental values by the method of least squares. 

L. Clark and L. Katz, Can. J. Research A21, 


vy Michels and C. Michels, Proc. Roy. Soc. pe A153, 201-14 
(1935); A160, 348-57 (1937). 

















